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FOREWORD 


A part  of  the  research  work  on  Contract  Nonr  393 (Oi)  has  bean  carried 
out  by  Mr.  Satoshi  Matsushima,  a graduate  student  participating  in  a co- 
operative research  program  in  solar  physics  between  High  Altitude  Observa- 
tory and  the  University  of  Utah.  The  following  document  is  the  Ph.D. 
thesis  summarizing  his  work.  It  is  issued  herewith  as  a contract  technical 
report.  The  thesis  was  jointly  supervised  by  the  astrophysics  group  of  the 
University  of  Utah,  the  research  staff  of  High  Altitude  Observatory,  and 
Dr.  Richard  N.  Thomas  of  Harvard  College  Observatory, 

The  thesis  reports  results  on  the  thermodynamic  structure  of  the  soiar 
atmosphere  derived  principally  from  the  observations  of  the  total  3olar 
eclipse  at  Khartoum,  Sudan,  on  25  February  1952.  The  ob servati onal  materials 
used  for  this  study  <*ero  reported  in  an  earlier  thesis  prepared  under  this 
contract  by  Dr.  Russell  Grant  Athay. 

The  entire  program  of  this  contract,  including  both  the  eclipse  observa- 
tions and  the  post-aclipse  analysis,  has  been  carried  out  in  close  coopera- 
tion with  the  Naval  Research  Laboratory.  This  work  on  the  thermodynamic 
structure  of  the  sun  from  optical  observations  has  great  relevance  to  the 
work  of  the  Naval  Research  Laboratory  on  the  same  problem,  using  solar  radio 
noise  observations.  We  wish  to  express  our  deep  appreciation  to  the  Naval 
Research  Laboratory  group  and  in  particular  to  Dr.  John  P.  Hagen  and  Dr. 

E.  0.  Hulburt,  for  their  invaluable  consultation  and  assistance. 

The  principal,  scientific  results  in  this  thesis  will  later  be  made  the 
subject  of  published  articles  in  the  scientific  literature,  and  reprints 
will  be  available  as  additional  technical  reports  under  this  contract. 


Walter  Orr  Roberts,  Director 
High  Altitude  Observatory 


8 May  195U 
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ABSTRACT 


Recent  observations  have  shown  that  the  classical  theory  of  the 
atelier  atmospheres  is  not  at  all  adequate  for  explaining  the  thermody- 
namic structure  of  the  solar  chromosphere.  In  this  thesis  a special 
effort  was  made  to  develop  the  theory  of  the  structure  of  an  atmosphere 
in  which  departures  from  thermodynamic  equilibrium  are  predominant,  A 
theory  of  the  phenomenon  of  self-absorption  was  also  treated.  These 
theoretical  developments  were  carried  out  with  the  aim  of  solving  various 
discrepancies  encountered  in  the  interpretation  of  the  flash  spectrum  by 
classical  theory.  The  resulting  theoretical  relationships  were  then 
applied  to  the  analysis  of  the  hydrogen  spectra  observed  at  the  1932  and 
1952  eclipses. 

In  general;,  if  the  kinetic  temperature  and  the  radiation  temperature 
in  an  atmosphere  are  not  identical,  the  distribution  of  the  atomic  energy 
states  can  not  be  described  by  the  equations  of  thermodynamic  equilibrium. 
The  departure  from  thermodynamic  equilibrium  may  be  measured  by  a dimen- 
tionless quantity,  bn,  which  is  a ratio  of  atomic  density  in  the  nth  energy 
level  to  that  given  by  the  Boltzmann-Saha  formula.  The  principle  of  the 
theoretical  calculation  of  bn,  developed  by  Menzel  and  Thomas,  was  employed 
in  devising  a method  of  taking  collisional  transitions  into  consideration. 
Numerical  calculations  were  carried  out  for  a pure  hydrogen  atmosphere, 
characterized  by  a kinetic  temperature,  Te  = 35,000°  and  illuminated  by  a 
radiation  field  at  a temperature,  Tr  = 6,000°.  It  appeared  that  the 
largest  differences  between  our  bn-values  and  those  that  Thomas  obtained 
were  in  bp  and  b2«  Thomas  neglected  collisional  excitation  from  levels 
other  than  the  first  level.  We  then  applied  our  numerical  values  of  bn  to 
explain  some  inconsistencies  that  appeared  in  laboratory  measurements  of 
hydrogen  arc  spectra. 


New  theoretical  expressions  for  the  total  Balmer  line  intensities 
observable  in  th9  eclipse  spectra  were  then  developed,  so  as  to  include 
the  effect  of  self-absorption.  The  resulting  formula,  for  total  line 
emission,  E^,  integrated  over  a volume  above  the  moon's  limb  iss 
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where  Nn(h)  denotes  the  total  number  of  hydrogen  atoms  in  the  n^h  level 
within  a unit  column  along  the  line-of-sight  at  a chromospheric  height, 
h,  and  a2n  g is  the  atomic  absorption  coefficient  at  the  line  center,  p 
represents’ the  exponential  height  gradient  of  hydrogen  atom.  Some  approxi- 
mate formulae  were  also  derived  empirically  that  agreed  with  ths  above 
expression  for  the  range  of  physical  parameters  involved  in  the  present 
problem. 


The  following  assumptions  were  made  to  perform  the  integrations 
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1.  The  density  of  hydrogen  atoms  decreases  exponentially  outwards 
along  a solar  radius. 

2.  The  logarithmic  gradient,  p,  is  assumed  to  be  the  same  for  all 
energy  states  of  hydrogen. 

3.  The  kinetic  temperature,  Te,  is  constant  throughout  the 
atmosphere. 

U.  The  central  intensity  of  the  lines  is  solely  determined  by  the 
Doppler  profile. 

The  formula  for  ^elf-absorption  thus  obtained  was  applied  to  the 
analysis  of  Balmer  line  intensities  of  the  1932  and  1952  eclipse  spectra 
to  determine  the  distribution  of  hydrogen  atoms  in  the  2nd  level,  N2,  and 
the  chromospheric  electron  temperature.  It  was  found  that  assumptions  1 
and  3 are  sufficiently  consistent  in  the  lower  region  of  the  chromosphere 
where  our  present  investigations  are  concerned.  Assumptions  2 and  U 
appeared  to  be  most  questionable.  At  the  present  stage  of  this  research 
it  seems  difficult  to  exauiine  quantitatively  the  possible  errors  due  to 
these  simplifications.  More  precise  calculations,  including  pn  on  n and 
fine  structures  of  line  profiles,  would  be  possible  only  by  numerical 
integration. 

Essentially  three  methods  were  developed  in  analyzing  the  1932  and 
1952  eclipse  spectra:  (a)  determination  of  the  population  of  the  Balmer 
ground  state  from  the  relative  amount  of  self -absorption  in  the  Balmer 
lines.  (The  basic  principle  of  this  method  was  first  considered  by  Thomas.) 
(b)  determination  of  the  amount  of  self-absorption  and  the  departure  from 
thermodynamic  equilibrium,  bn,  simultaneously  from  the  line  intensities 
only;  (c)  obtaining  the  amount  of  self-absorption  and  bn  by  combining  the 
line  intensities  and  the  emission  in  the  continuum;  then,  determining 
by  using  the  results  from  the  first  and  second  methods  to  examine  the 
internal  consistencies.  All  of  the  above  results  showed  good  agreement 
within  the  estimated  errors  of  the  observations. 

From  the  1932  data  it  was  found  that  N?,  the  total  number  of  atoms 
in  the  Balmer  ground  state  within  a unit  column  along  the  line- of- sight 
at  a chromospheric  height,  h,  and  the  exponential  height  gradient,  p, 
could  be  well  represented  by  the  following  formulae: 

log  N2  “ 16.13  + slog  -0.36  x 10"8  hcm;  p * 0.82  x 10-^cm-l 

Similarly  from  the  1952  data,  it  was  found  that: 

log  N2  - 16.37  (+0,023)  + i-log  -0.578  (+0.018)  x loA®"1 

p * 1.33  x 10”^cm"^ 

The  most  probable  value  of  N2  at  each  height  was  given  by  using  the 
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Tq  distribution  determined  from  the  continuum  analysis.  It  was  concluded 
that  the  physical  conditions  in  the  chromospheric  regions  observed  at  the 
1932  and  1952  eclinses  were  considerably  different*  Departure  from 
thermodynamic  equilibrium  were  not  predominant  in  the  1932  data  and  only 
the  effect  of  self-absorption  may  be  important  for  the  interpretation  of 
the  Balmer  decrement.  On  the  other  hand,  both  effects  appeared  euqally 
important  in  the  1952  data. 

From  the  1952  data  the  observed  values  of  bn  for  n higher  than  10 
were  found  to  be  almost  constant  in  the  region  lower  than  30,000  km. 
Detailed  calculations  showed  that  the  values  from  b^Q  to  bjjo  could  be 
expressed  by 

log  bn  = -0.13  (5). 022)  + 6.7  (+0.28)  J 


or 


log  bn  = 0.73  (+0.0U9)  - 0, 02U(5)9 0011)  n 

The  value  of  b2  at  each  height  was  also  determined  by  two  independent 
methods.  The  results  agreed  sufficiently  to  prove  the  accuracy  of  the 
methods. 

These  values  of  bn  provide  a measure  of  the  chromospheric  electron 
temperature,  Te,  and  electron  density,  Ne,  T0  was  found  to  increase  up- 
wards but  not  to  reach  30,000°  below  the  height  h - 3,000  km.  More  exact 
values  of  Te  could  not  be  obtained,  because  of  the  lack  of  theoretical 
values  of  bn„  It  was  therefore  suggested  that  the  theoretical  calcula- 
tions of  bn  should  be  carried  out  for  small  intervals  of  the  values  of 
T0  and  N0. 
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I.  INTRODUCTION 


It  has  long  been  knov;n,  even  before  quantitative  measurements  vere 
possible,  that  the  flash  spectrum  of  the  solar  chromosphere  has  various 
anomalous  features  that  could  not  be  explained  by  a simple  analogy  to 
the  classical  theory  of  stellar  atmospheres . One  might  expect,  from 
Kirchhoff's  law,  that  the  flash  spectrum  should  be  a reversal  of  the 
Fraunhofer  spectrum.  That  is,  the  relative  intensities  of  the  chromo- 
spheric emission  lineB  would  be  the  same  as  photospheric.  In  a more 
concrete  definition,  one  may  say  that  the  "classical  theory"  of  the 
chromosphere  is  a model  according  to  which  hydrostatic  equilibrium  holds, 
and  in  which  the  excitation  and  ionization  of  atoms  are  fully  described 
by  the  Boltzmann-Saha  formula  based  on  the  assumption  of  thermodynamic 
equilibrium  at  the  photospheric  radiation  temperature.  Thus,  once  the 
boundary  conditions  are.  set,  the  classical  model  completely  determines 
the  distribution  of  physical  parameters  from  which  one  can  predict  the 
quantitative  character  of  the  chromospheric  spectrum. 

A number  of  eclipse  observations  for  determining  the  relative  in- 
tensities of  lines  and  the  continuum  only  confirmed  the  above  mentioned 
descrepancies  with  the  classical  model.  The  first  quantitative  mea- 
surements were  made  at  the  1932  eclipse  by  CilDie  and  Kenzel  (1935). 

From  the  energy  distribution  in  the  Balmer  continuum  they  obtained  a 
value  for  the  electron  temperature  Te  = 5000°  K,  and,  using  this  tem- 
perature, they  computed  an  electron  density  Ne  - 2.3  x lCP-1  cm- 3 from 
the  Balmer  free-bound  emission. 

The  most  conspicuous  features  of  the  eclipse  spectrum  may  be  the 
anomalous  excitation  of  hydrogen  lines  and  the  presence  of  lines  of  He  I 
and  He  II  that  are  net  observed  in  the  photospheric  spectrum.  Despite 
the  evidence  of  such  superexcitation  of  H and  He,  Wildt  (1947)  obtained 
a value  of  Ne  similar  to  that  of  Gillie  and  Menzel  by  an  application  of 
the  Inglis-Tell er  formula  (1939)  for  the  resolution  of  hydrogen  lines. 
Further,  Goldberg's  (1939)  determination  of  the  excitation  temperature 
from  the  relative  population  of  excited  helium  levels  indicated  similar 
results,  except  for  the  upward  increase  of  temperature.  For  thermo- 
dynamic equilibrium,  the  relative  population  of  atoms  in  excited  levels 
is  given  by  the  Boltzmann-Saha  relation  for  a given  Te  and  Ne.  Hence, 
unlees  we  assume  an  unreasonably  low  electron  density,  the  above  value 
of  Te  is  much  too  low  to  produce  the  lines  of  such  high  excitations. 

The  appearance  of  most  of  the  ionized  metallic  lines  as  more  en- 
hanced than  those  in  the  Fraunhofer  spectrum  was  also  considered  as 
indicating  anomalous  excitation.  Woolley  and  Stibbs  (1953),  for  example, 
insist  that  the  ionization  of  metals  in  the  chromosphere  does  not  require 
a temperature  higher  than  that  of  the  photosphere,  and  that  enhanced 
lines  can  be  adequately  accomited  for  by  a decrease  in  pressure.  But, 
they  state  that,  this  cannot  be  extended  to  the  hydrogen  and  helium 
excitations. 
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Thus,  the  fii’st  inconsistency  with  the  classical  model  appears  in  the 
large  difference  of  intensities  between  the  higher  excitation  lines  and 
the  metallic  lines.  Miyamoto  (1947) , and  the  author  and  Miyamoto  (1947) 
considered  the  possibility  that  the  ultraviolet  radiation  emitted  inwards 
from  the  highly  heated  corona  causes  these  chromospheric  excitations,  so 
that  it  is  not  necessary  to  have  a radiation  temperature  higher  than  the 
photospheric  temperature.  They  computed  the  amount  of  free-free  emission 
due  to  the  free  electrons  of  velocities  corresponding  to  Te  = 1,000,000°  K. 
The  energy  values  thus  obtained  at  the  frequencies  corresponding  to  the 
ionization  potentials  of  He  I and  He  II  were,  respectively,  equal  to 
energies  at  these  same  frequencies  in  8000°  and  16,500°  black  body  radia- 
tion. The  computed  values  were,  however,  inadequate  for  explaining  the 
ionization  of  hydrogen. 

The  next  direct  observational  contradiction  to  the  classical  model 
is  associated  with  the  assumption  of  hydrostatic  equilibrium.  The  observed 
limb  darkening  shows  evidence  that  the  photosphere  is  in  hydx'ostatic  equili- 
brium, under  the  a ssumption  that  the  negative  hydrogen  ion  is  responsible 
for  the  opacity  in  the  continuum.  If  we  a ssume  a layer  in  hydrostatic 
equilibrium  as  an  extension  of  the  photosphere  and  calculate  the  logarith- 
mic density  gradient,  yi  =* * v,e  £.7  ' / o'g  for  T = 5000°  K 

and  the  mean  molecular  weight,^  « I . Contrary  to  this  theoretical  pre- 
diction, the  observed  intensity  decrement  along  the  geometrical  height 
appears  much  too  low  for  the  theoretical  value.  (For  thermodynamic 
equilibrium  and  no  self-absorption,  the  emission  gradient  is  equal  to  the 
density  gradient.)  According  to  Cillie  and  Menzel  (1935),  for  example, 
y3  for  neutral  hydrogen  at  a height  of  about  1000  1cm  above  the  limb  is 
1.54  • 10~°  cm”-*-,  which  agrees  well  with  the  value  1.62  • 10"®  cnr®- 
obtained  by  Pannekoek  and  Minnaert  (1928)  from  the  1927  eclipse  data. 

Again  under  the  assumption  of  hydrostatic  equilibrium,  these  hej.ght  gra- 
dients require  a temperature  of  the  order  20,000°  K,  Wildt  (1947),  using 
an  ingeneous  method  covering  a number  of  previous  observations,  gets 
A ” 0.^  ■ /0~g  v/hich  leads  to  the  temperature  as  high  as  35,000°  K. 

The  disagreement  of  Wildt* s value  with  these  previously  obtained  is  dis- 
cussed by  Thomas  (1950b).  He  argues  that  the  omission  of  the  departure 
from  thermodynamic  equilibrium  would  lower  the  value  of  /3  in  Wildt' s 
method.  Since  Wildt* s determination  is  based  on  the  comparison  of  lines 
of  different  quantum  number,  the  positive  quantity  due  tc  the  departure 
from  thermodynamic  equilibrium  should  enter  directly  into  the  density 
gradient.  Further,  the  eclipse  observations  show  that  the  emission  gra- 
dients of  hydrogen  lines  increase  systematically  towards  higher  members 
of  the  Balmer  series.  As  we  shall  discuss  in  the  later  chapters  of  this 
paper,  this  seems  to  be  one  of  the  pieces  of  evidence  indicating  that  the 
chromosphere  is  not  in  a state  of  thermodynamic  equilibrium.  Thus,  the 
y important  fact  of  these  eclipse  observations  is  that  there  exists  a sharp 

line  of  demarcation  between  the  chromosphere  and  the  photosphere,  where 
the  assumption  of  hydrostatic  equilibrium  must  break  down  unless  we 
assume  a temperature  as  high  as  20,000  - 30,000°. 

In  addition  to  the  above  lover  emission  gradient,  the  observed  low 
decrement  of  balmer  line  intensities  towards  the  series  limit  has  been 


well  known  as  one  of  the  anomalous  features  of  the  flash  spectrum.  As 
we  shall  discuss  in  dei^il  in  Chapter  IV,  Thomas  (I950e ) has  studied 
this  descrepancy  and  shown  that  it  is  mostly  due  to  the  ignoring  of  the 
effect  of  self-absorption,  in  earlier  considerations  of  the  problem. 

In  the  above  we  have  mentioned  various  observational  evidence 
suggesting  that  the  actual  chromosphere  departs  considerably  from  the 
simplest  model  in  which  it  is  considered  as  an  extension  of  the  photo- 
spheric  surface,  and  have  pointed  out  that  spectroscopic  observations 
may  be  better  explained  by  assuming  a higher  temperature  in  the  chromo- 
sphere than  in  the  photosphere. 

Since  the  existence  of  a coronal  electron  temperature  of  the  order 
of  a million  degrees  is  well  established,  it  is  natural  to  assume  a tem- 
pera ture  increase  outward.  That  is,  there  must  exist  a layer  above  the 
photosphere  in  which  the  temperature  gradient  changes  its  sign.  The  only 
question  is  whether  the  minimum  temperature  is  at  the  base  or  above  the 
chromosphere. 

The  determination  of  kinetic  temperature  by  Redman  (1942a,  1942b) 
appeared  most  likely  to  give  a definite  solution  for  this  question. 

Redman  measured  the  line-profiles  of  several  intense  lines  of  the  1940 
eclipse  spectrum,  from  the  relative*  broadening  of  hydrogen  lines  and  the 
lines  of  heavier  elements,  he  concluded  that  .the  line  width  is  mainly  due 
to  Doppler  broadening  by  thermal  motion  of  atoms  corresponding  to  a kinetic 
temperature  of  30,000°  K,  rather  than  turbulent  motion.  These  observations 
apply  to  chromospheric  heights  as  low  as  1300  km.  On  the  other  hand,  tem- 
peratures in  the  photospheric  layers  as  obtained  from  limb  darkening  in 
the  continuum,  excitation  of  molecular  or  metallic  lines,  Balmer  line 
profiles,  etc.  indicated  (following  recent  reviews  by  Minnaert  (1953)  and 
Unsdld  (1952))  a temperature  of  3,800  to  4,500°  in  the  external  photo- 
spheric  layers.  Similar  results  were  obtained  by  Pecker  (1951)  and  others 
from  the  theory  of  radiative  equilibrium  and  the  blanketing  effect.  How- 
ever, contrary  to  earlier  interpretations  of  the  broadening  of  Fraunhofer 
lines  as  due  to  turbulent  motion,  Bell  (1951)  insists  that  it  is  caused 
by  thermal  motion  at  a kinetic  temperature  of  at  least  10,000°  K.  She  has 
examined  the  half-widths  of  line  profiles  made  by  elements  of  vai’ious 
atomic  weight  and  found  that  the  Doppler  width  varies  with  atomic  weight 
in  such  a way  as  to  indicate  a thermal  motion  corresponding  to  Te  = 10,000° 
for  high  excitation,  17,000°  for  low  excitation  neutral  lines  and  20,000° 
for  ionized  lines.  The  value  of  17,000°  is  also  confirmed  by  her  curve- 
of-growth  results.  Further,  Bell  considered  that  the  systematic  change 
of  Te  with  excitation  potential  shows  that  the  kinetic  temperature  in- 
creasing upwards,  since  other  evidence,  such  as  the  heights  to  which 
various  lines  extend  in  the  chromosphere,  agrees  in  indicating  that  the 
lower  excitation  lines  are  formed  at  a greater  e levation  jn  the  atmosphere 
than  the  higher.  Also,  ionization  increases  outward  and  the  ionized  lines 
have  the  highest  kinetic  temperature  of  all.  She  states  that  the  value 
of  17,000°  K agrees  well  with  the  previous  measures  by  Allen  (1937,1949), 
Menzel  and  Dell  (1948 ) , Pierce  and  Goldberg  (1948),  Shane  (1941),  van  de 
Hulst  (1946),  etc.,  except  that  in  the  past  the  interpretation  has  been 
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in  terms  of  turbulent  velocities.  Bell,  however,  finds  no  evidence  for 
velocities  independent  of  atomic  weight.  Bell's  assumption  that  she  may 
ignore  the  detailed  distribution  of  atoms  and  the  ambiguity  of  her  method 
of  measuring  the  half-widths  are  criticized  by  J.  Houtgast  (1952)  and 
others . 

While  Beil's  result  is  thus  open  to  question,  it  may  be  considered 
evidence  to  support  Redman's  hypothesis,  indicating  that  the  reversal  of 
temperature  may  lie  even  as  low  as  the  reversing  layer.  Since  the  tem- 
peratures discussed  above  are  all  kinetic  tempera tures  determined  from  the 
lines  individually,  they  are  not  necessarily  contradictory  with  the  radia- 
tion tempera ture  which  characterizes  the  continuum  and  the  relative  inten- 
sities of  Fraunhofer  lines.  Ultimately  we  must  distinguish  between  radia- 
tion temperature,  kinetic  (electron)  temperature,  excitation  tanperature, 
ionization  temperature,  etc.  These  temperatures  should  be  identical  only 
when  the  atmosphere  is  in  a state  of  strict  thermodynamic  equilibrium. 

If  we  postulate  that  the  physical  state  of  the  chromosphere  departs  from 
thermodynamic  equilibrium,  we  must  specify  the  physical  conditions  in 
order  to  relate  the  temperatures  defined  in  different  ways. 

The  question  of  how  the  energy  necessary  to  raise  the  kinetic  tem- 
perature above  the  photospheric  temperature  is  supplied  now  arises.  This 
problem  has  been  studied  by  several  authors  including  Bierman  (1948), 
Schwarzschild  (1943),  Thomas  (1948a),  Alfven  (1947,  1950),  Hoyle  (1949), 
Miyamoto  (1949),  etc.  We  discuss  only  a theory  by  Thomas  which  starts 
more  directly  from  observed  facts  than  is  the  case  with  other  more  specu- 
lative considerations.  Thomas  interprets  the  observation  of  chromospheric 
spicules  as  a system  of  superthermic  jets.  These  jets  supply  the  energy 
needed  to  give  the  observed  high  kinetic  temperature  in  the  chromosphere. 

His  quantitative  calculation  based  on  hydrodynamics  predicts  theoretically 
a configuration  of  the  spicules  which  agrees  with  the  essential  features 
observed  by  W.  0.  Roberts  (1945)  and  shows  that  the  energy  transferred 
from  the  system  to  the  chromosphere  appears  adequate  to  maintain  a 35,000° 
kinetic  temperature. 

Further,  Thomas  (1943b,  1949a,  b,  c,  1950b,  1952)  has  developed  the 
hypothesis  of  a 35,000°  kinetic  temperature  to  treat  the  radiation  field 
of  the  hydrogen  chromosphere  and  has  explained  the  anomalous  excitation 
of  the  Balmer  lines  mentioned  earlier.  He  assumed  a pure  hydrogen  atmos- 
phere with  a kinetic  temperature  of  35,000°  which  is  illuminated  by  a 
photospheric  radiation  temperature  of  6000°.  In  such  a high  kinetic  tem- 
perature, the  effect  of  collisional  excitation  or  ionization  becomes 
predominent,  causing  a severe  departure  from  thermodynamic  equilibrium, 
so  that  the  distribution  of  populations  of  each  energy  state  does  not 
* obey  the  Maxwell-Boltzmann  distribution.  As  will  be  discussed  in  detail 

in  the  later  chapters  of  this  thesis,  Thomas  has  shown  that  in  general  as 
the  kinetic  temperature  Increases  a larger  number  of  atoms  remain  in 
lower  energy  states  than  there  would  be  in  the  case  of  thermodynamic 
equilibrium.  After  solving  the  cyclic  (steady  state)  eqxiations  for  each 
transition  of  hydrogen  including  the  photoelectric  and  collisional  transi- 
tions, the  numerical  results  by  Thomas  show  a large  excess  of  atoms  in  the 
Balmer  ground  state  for  the  possible  range  of  chromospheric  electron  densities 
assumed. 


The  same  problems  were  treated  independently  by  Giovanelli  (1948a,  b, 
c,  d)  by  a slightly  different  method.  He  computed  the  population  of  the 
Balmer  ground  state  directly,  while  Thomas  accounted  for  its  excess  over 
the  Boltzmann  distribution  by  a parameter,  bn.  Numerical  comparison  of 
the  two  methods  shows  comparable  results. 

Thus  while  Redman's  hypothesis  had  been  supported,  the  recent  develop- 
ment of  radio  observation  provided  another  means  of  measuring  the  electron 
temperature  of  the  solar  atmosphere.  The  electron  temperature  may  be 
determined  from  the  total  amount  of  energy  emitted  at  a certain  wave  length 
from  above  a certain  height  in  the  sun's  atmosphere,  provided  the  distribu- 
tion of  electron  densities  in  the  corresponding  region  is  known.  Hagen 
(1951)  observed  the  quiet  sun's  intensity  at  a wavelength  of  8.5  mm  and 
obtained  the  equivalent  temperature  due  to  the  thermal  radiation  through- 
out the  chromospheric  and  coronal  region.  He  made  an  analysis  of  these 
measurements  in  conjunction  with  many  other  measurements  at  longer  wave 
lengths  to  determine  the  electron  temperature  gradient.  Since  the  observed 
equivalent  temperature  is  an  integrated  radiation  throughout  the  atmosphere 
ebove  a certain  height,  we  have  to  assume  the  electron  density  distribu- 
tion to  reduce  the  observed  temperature  to  tho  kinetic  temperature  at  each 
layer.  Hagen  adopted  the  values  of  Me  obtained  by  Baumbach  (1937),  Allen 
(1947)  and  van  de  Hulst  (1947)  for  the  corona  and  Wildt's  (1947)  values 
for  the  chromosphere.  According  to  Hagen,  Wildt's  Ne  distribution  in  the 
chromosphere  is  somewhat  hypothetical,  except  perhaps  for  the  density  (or 
rather  upper  limit  to  the  density)  of  1.74  • 10^  electrons  per  cubic  cm 
at  a height  of  500  km.  Assuming  such  a model,  Hagen  reached  the  conclu- 
sion that  the  electron  temperature  is  fairly  uniform  at  a value  less  than 
10,000°  K to  a height  of  10,000  Ion  and  then  rises  steeply  to  a temperature 
of  1,000,000°  K at  a height  of  25,000  km.  To  make  radio  noise  data  agree 
with  Redman's  temperature  value  at  1500  km,  it  is  necessary  to  decrease 
the  absorption  coefficients  so  that  the  principal  absorption  for  the 
shorter  wavelengths  will  occur  at  much  lower  heights.  Since  the  absorp- 
tion coefficient  varies , directly  with  Ne^  and  inversely  with  Tfc.3/2,  either 
an  increase  in  the  assigned  temperature  of  the  chromosphere,  a decrease 
in  the  electron  density,  or  some  conbination  of  these  two  is  indicated. 
Hagen's  calculations  show  that,  in  order  to  satisfy  simultaneously  an 
existence  of  a kinetic  temperature  of  30,000°  K at  a height  of  1500  km 
and  the  observed  equivalent  temperature  at  the  \arious  radio  wave  lengths, 
the  electron  density  in  the  lower  chromosphere  must  be  reduced  to  a value 
about  one-tenth  of  the  one  he  first  assumed.  This  would  change  the  posi- 
tion of  maximum  temperature  gradient  from  the  top  of  the  chromosphere  to 
a point  in  the  chromospheric  region. 

Although  the  uncertainty  in  the  Ne  distribution  made  the  interpreta- 
tion of  the  radio  measurements  somewhat  indeterminate,  Hagen's  results  at 
one  time  appeared  to  offer  an  important  contradiction  to  Redman's.  Many 
authors  used  this  radio  result  as  an  argument  for  temperatures  as  low  as 
5000  - 6000°  X for  the  chromosphere.  Since  the  temperature  determined  by 
Redman  is  the  kinetic  temperature,  and  that  determined  by  radio  noise  is 
the  electron  temperature,  we  inquire  whether  or  not  we  should  expect  the 
two  to  agree.  We  may  define  the  electron  tempera ture  as  a measure  of  the 
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thermal  motion  of  electrons  and  the  kinetic  temperature  as  the  same 
function  of  the  motion  of  the  atoms.  If  we  assume  a Maxwellian  distri- 
bution of  velocities  and  equipartition  of  energy,  then  these  two  tem- 
peratures must  be  identical.  Thomas,  Krook,  Menzel  and  Bhatnager  (1953) 
have  concluded  from  theoretical  considerations  that  Te  cannot  differ 
significantly  from  Tjj  for  Te  1C$  . K. 

A number  of  objections  have  been  raised  by  those  who  criticize  the 
hypothesis  of  the  higher  kinetic  temperature  in  the  chromosphere.  The 
most  direct  objection  to  Redman’s  interpretation  was  made  by  S.  Miyamoto 
(1951b).  He  pointed  out  that  the  omission  of  self-absorption  in  Redman’s 
measurement  may  cause  an  erroneous  interpretation  of  the  line-profiles. 
Since  self-absorption  is  more  effective  in  the  central  portion  of  a lino, 
it  would  tend  to  increase  the  half-widths  of  the  observed  line-profiles. 
Miyamoto  calculated  the  half-width  of  H^,  H p,  and  II  y including  self- 
absorption and  has  shown  that  Redman's  measurements  agree  with  his  results, 
if  T is  set  as  5,700°.  This  suggestion  has  been  widely  accepted  by  several 
authors  (van  de  Hulst  1954)  including  Redman  (1952)  himself. 

It  seems  that  most  of  the  criticisms  of  Redman’s  hypothesis  are  based 
on  the  study  of  the  excitations  and  ionizations  of  metallic  lines.  While 
the  high  intensity  of  ionized  metallic  lines  compared  to  the  neutral  had 
been  considered  at  one  time  as  one  of  the  anomalies  of  the  flash  spectrum, 
Menzel  (1931)  and  Woolley  and  Stibfcs  (1953)  have  shown  that  it  can  be 
attributed  to  the  decrease  of  election  pressure  in  the  chromosphere  and 
does  not  require  a higher  temperature.  Their  detailed  calculations, 
using  the  ionization  formula,  have  explained  the  observed  relative  inten- 
sities of  ionized  metallic  lines  to  those  of  neutral  lines.  Miyamoto 
and  Kawaguchi  (1950)  compared  the  metallic  line  intensities  in  the  flash 
spectra  with  those  predicted  by  theories  of  metallic  ionization  and  the 
emission  curve  of  growth  for  the  Ca  II  K line.  They  have  shown  that  for 
Te  = 35,000°  the  metallic  ionization  by  electron  collision  would  be  far 
too  great  and  it  would  resemble  that  of  the  atmospheres  of  the  earliest 
type  stars.  It  is  concluded  that  the  observed  emission  curve  of  growth 
obtained  from  me  1930  eclipse  for  the  Ca  K line  is  best  represented  by 
the  theoretical  curve  for  Te  = 6000°.  Further,  Miyamoto  (1951b)  has 
calculated  the  ultraviolet  emission  from  the  chromosphere  of  Te  = 35,000°, 
in  vrhich  he  argues  that  for  such  a high  temperature,  hydrogen  must  be  so 
completely  ionized  that  the  abundance  ratio  of  hydrogen  to  metals  does 
not  coincide  with  the  observed  value,  10^  : 1. 

VJurm  (1942a,  b)  has  considered  the  presence  or  absence-  of  the  nebular 
forbidden  lines  of  oxygen  and  nitrogen  to  be  a criterion  of  chromospheric 
electron  temperature.  He  has  computed,  as  a function  of  the  electron 
temperature,  the  theoretical  intensity  ratio  of  forbidden  transitions  by 
collision  for  the  ground  states  of  01,  Oil,  0111,  and  Nil  to  the  photo- 
spheric  light  scattered  by  electrons.  He  shows  that  these  ratios  are 
less  than  unity  for  Te  ~ 5000°,  but  are  the  order  of  10^  for  Te  higher 
than  20,000°.  Thus  Wurm  argues  that  none  of  these  forbidden  lines  should 
be  seen  if  the  chromospheric  temperature  is  of  the  order  of  5000°,  but, 
instead,  all  the  lines  should  be  present  if  the  temperature  is  as  high  as 


35,000°.  Since,  in  spite  of  the  non-appearance  of  these  lines,  there  must 
be  some  layer  in  the  chromosphere  where  Te  reaches  to  35,000°,  Vurrn  con- 
cludes that  the  layers  at  such  high  temperature  are  exceedingly  thin;  that 
is,  the  temperature  rises  very  steeply  in  the  upper  region  of  the  chromo- 
sphere as  is  indicated  by  radio  noise  observations. 

This  conclusion  of  Wurm  has  been  supported  by  Woolley  and  Allen  (1950). 
They  considered  a theoretical  model  of  the  corona  and  the  chromosphere, 
trying  to  make  a best  fit  with  observed  data.  Woolley  and  Allen  compute 
the  amount  of  the  ultraviolet  (Lyman  continuum)  radiation  emitted  from 
their  model  as  a function  of  the  electron  temperature  and  conclude  that 
the  Lyman  continuum  emission  from  the  low  temperature  chromosphere  just 
satisfies  the  observed  excitation  of  the  ionosphere,  but  for  temperature 
as  high  as  35,000°  K the  theoretical  level  of  ionization  would  be  much 
higher  than  is  observed.  As  a model,  they  assumed  Baumbach's  electron 
density  distributions,  setting  the  temperature  at  5000°  K up  to  a height 
of  5000  km,  then  rising  to  the  order  of  1,000,000°  within  a layer  of  a 
few  hundred  kilometers.  This  gives  a best  fit  with  the  radio  results. 

As  an  explanation  of  the  appearance  of  He  I and  He  II  linos  that  would 
be  consistent  with  Redman’s  observations,  Woolley  suggests  that  the  actual 
chromosphere  is  not  homogeneous  hut  is  so  irregular  that  any  particular 
line-of-sight  in  the  flash  observation  samples  a variety  of  hot  and  cold 
spots.  That  may  be  the  reason  why  the  high  excitation  lines  like  He  I and 
He  II  appear  at  the  seme  height  as  the  lowest  excitation  metallic  lines. 

On  the  other  hand,  Miyamoto  (1951a)  has  examined  the  excitation  of 
helium  under  the  assumption  of  Te  - 35,000°  and  hydrogen  helium  ratio 
5:1.  He  computed  the  intensities  of  the  resonance  lines  of  He  I and  He  II, 
using  the  optical  depth  in  fhe  chromosphere,  and  shoved  that  the  popu- 
lations of  the  excited  levels  of  He  I and  He  II  should  become  so  large 
that  the  strong  absorption  lines  of  these  elements  would  be  present  in  the 
Fraunhofer  spectrum.  By  such  arguments  he  concludes  that  the  study  of  the 
helium  linc-s  by  L.  Goldberg  (193C.  ) indicates  a kinetic  temperature  not 
exceeding  10,000°.  His  result  appeared  favorable  for  the  lower  tempera- 
ture at  lower  heights  of  the  chromosphere,  although  the  temperature 
appears  to  increase  upwards  already  in  the  lower  layers. 

Finally  we  note  the  temperature  determination  from  the  Balmer  dis- 
continuity. Zanstra  ( 19 50a,  b)  has  attempted  to  separate  the  free-bound 
emission  and  the  scattered  radiation  by  free  electrons  at  the  Balmer  and 
Paschen  series  limit.  Since  the  capture  component  of  Paschen  continuum 
is  sensitive  to  the  temperature  and  decreases  rapidly  toward  shorter  wave 
lengths  for  Te  ~ 5000°,  the  intensity  jump  at  the  Balmer  limit  .is  larger 
for  lower  than  for  higher  temperatures.  After  examining  the  1926  eclipse 
observation  by  Davidson  and  Sti'atton  (1927),  Zanstra  finds  the  disconti- 
nuity is  considerably  larger  than  the  theoretical  upper  limit  for  Te  - 
30,000°.  More  precise  measurements  were  made  by  Koelbloed  and  Veliman 
(1951),  who  obtained  the  upper  limit  of  Te  as  15,000°  - 20,000°.  A 
similar  problem  has  been  also  discussed  by  S.  Ueno  (1950)  based  on  the 
low  chromo s ph er i c temperature. 
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We  have  summarized  above  the  important  studies  related  to  the  chromo- 
spheric temperature  made  before  the  1952  eclipse.  We  have  seen  that  there 
exist  many  difficulties  and  inconsistent  observations  in  the  path  of  a 
self-consistent  theory  of  the  chromosphere.  A large  number  of  expeditions 
vent  to  Khartoum,  Sudan,  for  the  eclipse,  February  25,  1952,  aiming  to 
resolve  these  discrepancies.  While  the  numerical  results  of  these  obser- 
vations have  not  yet  been  published  at  the  time  this  thesis  was  written, 
ve  may  discuss  the  preliminary  results  for  two  important  observations, 

Hagen  (1953)  made  solar  radio  noise  observations  at  8 mm  wave  length  and 
examined  the  intensity  change  with  the  moon's  motion.  He  concluded  that 
his  results  implied  a non-homo gen eous  model  of  the  chromosphere.  He 
suggests  that  the  radio  observations  may  not  be  fully  explained  by  a uni- 
form chromosphere  with  a mono tonic  variation  of  temperature  upwards. 

According  to  his  explanation,  the  chx'omo  spheric  model  may  be  best  described 
as  a mixture  of  spicule-like  radial  columns  of  different  physical  conditions, 
namely,  a mixture  of  columns  with  high  temperature  and  low  density  and 
columns  with  low  temperature  and  high  density.  The  integrated  emission 
observed  at  radio  wave  length  gives  a sort  of  averaged  temperature.  Redman's 
spectra  obtained  at  the  1940  eclipse  may  have  been  representative  of  a 
high  temperature,  region. 

At  the  1952  eclipse,  Redman  (1952.)  repeated  his  1940  observations. 

He  obtained  eight  exposures  for  each  of  the  wave  length  regions,  5550A  - 
6040  A and  3430A  - 4110  A.  Thus  the  Balmer  lines  up  to  H32  are  covered 
by  his  spectrum.  Admitting  the  suggestion  by  Miyamoto  (1951b)  Redman 
states  that  the  strong  self-absorption  for  the-  early  members  of  the  Balmer 
series  may  cause  an  erroneous  interpretation  of  the  width  as  due  to  thermal 
motion  only,  although  the  shape  of  the  line  may  not  reveal  self-reversal 
like  in  metallic  lines.  He  proposed  that  this  trouble  can  be  reduced  by 
going  to  higher  members  of  the  series.  This  may  lead  to  difficulty  from 
the  broadening  by  Stark  effect,  although  Redman  observed  eh at  the  lines 
do  not  grow  broader  and  broader  towards  the  series  limit,  as  has  been 
suggested.  Allowing  for  the  effects  of  both  self-absorption  and  Stark 
effect,  Redman's  preliminary  results  show  some  modification  of  his  former 
value.  That  is,  he  declares  that  the  kinetic  temperature  estimated  from 
the  higher  members  of  Balmer  series  and  the  weak  metallic  lines  cannot 
exceed  17,000°  K.  A higher  value  similar  to  earlier  measurements  is  shown 
for  the  lowest  Balmer  lines.  This  seems  to  imply  the  effect  of  self- 
absorption. 

After  this  thesis  research  had  been  completed,  we  received  a prelimin- 
ary summary  of  the  attempt  by  Woltjor  (1954)  to  construct  a spicule-like 
model  for  the  chromosphere,  as  has  been  suggested  by  Hagen  (1953).  Woltjer 
measured  intensities  of  spicules  and  the  parts  between  spicules  on  a Ha 
monochromatic  plate  taken  by  Lyot.  He  found  the  quotient  of  the  two 
intensities  to  be  1.9  at  5000  km.  He  deduced  a model  in  which  the  larger 
part  of  the  chromosphere  has  a temperature  of  some  5000°,  vhereas  for 
spicule  regions  it  increases  from  8000°  at  a height  of  1000  km  to  20,000° 
at  4000  km,  then  remains  constant  up  to  7000  km.  He  concludes  that  the 
spicules  cover  Z:/j  of  the  solar  surface  and  the  entire  Balmer  emission  is 
coming  from  these  spicules.  For  the  spicule  temperature  value  at  2000  km 
Woltjer  proposes  T «=  10,000°.  This  is  still  far  too  low  compared  to 
Redman's  determination. 
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Lo celly  Intensified  regions  have  been  already  netted  by  Gillie  and 
iSenzel  (1935)  in  the  flash  spectrum  of  1932.  Aih&y,  Evans  and  Roberts 
(1953)  suggest  that  a sln.il ar  hot  region  appears  on.  the  thirl  contact 
spectre  from  the  1952  eclipse.  The  Irregularity  of  solar  surface  evi- 
denced by  such  phenomena  as.  granulation,  suns pets,  spicules,  prominences, 
flares,  radio  noise  observations,  etc.  cay  cause  many  of  the  contradic- 
tions in  the  interpretation  of  chromospheric  observations,  thus,  the 
future  theoretical  model  of  the  chromosphere  ray  have  to  be  based  on  a 
non-uni foru  structure. . 


•O' 


Store  recent  theoretical  studies  are  attacking  these  pro  biers  by 
•ccnsi dering  the  energy  transport  by  a ragneto  hydroctynaralc  wave  or  the 
destruction  of  Hue  state  of  equilibrium  by  a shock  wave,  ill  these 
theories  consider  such  localized  phenomena.  on  the  sun's  surface  as  the 
magnetic  field,  turbulent  notions,  etc.  is  has  been  pointed  cut  by 
T*fco3Jey,  however,  there  is  vide  freedom  in  the  assumptions  if  we  choose 
an  irregular  model.  Unless  more-  direct  observations  like  Hagen's  are 
accuralated  in  the  future,  and  a more  confirmable  choice  of  the  distri- 
bution of  irregularities  fceccr.es  possible,  it  is  very  difficult  to  reach 
an  appropriate  node!. 


Along  with  the  expeditions  by  Hagen  and  fcedteam  to  the  .1952  eclipse, 
that  fcy  Flgb  Altitude  Conservatory  under  the  supervision  of  Roberts,  Ewans 
and  lac. jscs,  v&s  especially  aim’d  at.  the  deterulnotion  of  the  chromospheric 
taupe  rsiure  distribution  as  well  e s the  electron  density  distribution. 

Highly  improved  jumping  film  technique  .trade  It  possible  to  obtain  the 
flash  spectrum  with  much,  greater  height  resolution.  Tne  spectra  thus 
obtained  covered!  the-  wave  length  region  iron  3, GOG  k to  3,400  A,  with  & 
dispersion  of  7^  A ssi~x  in  the  ultraviolet  avd  12  A naT2-  in  the  visible  avsd 
infrared.  lfce  intensities  of  the  Balder  lines  end  the  Balner  continuum 


were  obtained  approximately  every  ISO  Iscl  firm  the  height  500  k®  above 
tie  base  of  the  cfceoBsosffcere.  Ike  major  part  of  the  present  thesis  is 
aimed  at  the  theoretical  reduction  of  - aimer  line  Intensities  frets,  these 
data  together  with  those  of  the  1932  eclipse  by  Gillie  and  Kenzel.  4*e 
are  specifically  interested  In  determining  the  electron  temperature  and 
the  population  of  the  hydrogen  ground  state,  and  the  variation  of  these 
quantities  with  geometrical  height.  The  net.  c-de  developed  in  this  thesis 
are  related  pa'ILicelsrly  to  the  problem  of  depertnree  from  th ersodyosssi c 
equilibrium  and  the  effects  of  sell— absorption  in  the  ehre-rso sphere,  these 
effects  have  been  ignored  in  the  earlier  reductions,  as  is  discussed  above. 

In  the  following  chapter,  the  theoretical  calculations  of  departures  from 
th eansjodyman i a equilibrium  ere  treated.  Chapter  III  consists  cf  the 
development  of  theoretical  forsmlae  for  self-absorption  in  an  optically 
thin  atmosphere,  and  the  naisnericel  results  are  quoted.  ’Chapters  17  and  7 
are  devoted  to  the  reduction  cf  the  1932  srri  1952  eclipse  data,  respectively. 
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II.  TEE  DEPASTURE  IHGK  TKEHIDDIliAKIC  EQUILIBRIA  I!f  TEE  CHROMDSFriERE. 


|1.  Introduction. 


In  the  study  of  the  radiation  field  of  the  chroisDBpfcere,  it  has  been 
shown  by  Keosel  and  Gillie  (1937)  that  the  departures  of  the  populations 
of  the  atonic  energy  levels  from  the  th erne :r rati c equilibrium  values  can 
be  expressed  by  a not— dir lemsi cnal  quantity  bn.  feq  is  the  ratio  of  popula- 
tion, lin,  In  atonic  state  n,  to  the  population  under  conditions  of  thersao— 
dynanic  equilibrium  at  the  sane  electron  density  and  temperature.  lienee 
It  is  defined  by  a swdifi cation  of  the  holtzrann-Safea.  formulas 


N.  = k N:  H, 


< 2KVi?Tt  I 


h 


*1 

A 


e^Te 


(1-1) 


where  IJ,If  ]«j,  and  Be  denote  the  density  respectively  of  hydrogtn  atoms  in 
the  nth  quantum  level,  of  lens,  and  of  electrons.  Xjg  Is  the  ionization 
sciential  frer.  the  level  a.  ‘Jhe  other  notation  is  standard. 


Under  the  assumption  that  within  a small  region  of  the  atmosphere 
there  exists  a unique  kinetic  temperature,  Te,  defined  by  a Maxwellian 
distribution  of  velocities  of  free  electrons,  we  can  note  i-rediately 
that  bn  approaches  unity  when  n ter.es  to-  infinity. 

An  extensive  study  of  the  departure  iron  themso^iuunic  cquilibri-mi 
in  the  solar  chromosphere  in  tarns  of  by,,  has  been  carried  out  and  published 
ir  a,  series  of  papers  by  Thorsas  (2943c,  '1949a,  b,  e,  ISSGb,  1952).  As 
ve  shall,  discuss  later  in  detail,  Thomas  obtained  the  result  that  for  the 
case  in  which  Ie  exceeds  the  radiation  iosperature,  ca  is  greater  than 
unity  and  isonotonical  ly  decrees  to  to  unity  as-  n Increasrs.  Recently , 
however,  a few  theoretical  computations  have  been  in  disagreement  with  this 
conclusion.  lor  example,  Chamberlain  (1953)  computed  the  values  of  bj. 
from,  'the  same  equations  cf  statistical  equilibrium,  but  at  very  scuch 
lower  values  of  He.  For  Te  = 10,000°  K,  he  obtained  bj,  all ght2y  less 
than  unity  for  optically  thick  models  and  considerably  less  for  an  optically 
thin  Kediun.  It  is  not  certain  that  these  contradictory  results  have 
physical  significance.  Since  the  methods  of  deriving  bjt-cxpresoJoms  were 
practically  the  same  for  both  casts,  the  disagreement  eeens  due  tc  the 
uncertainty  of  the  factors  mploycd  for  the  ccnputailon,  such  as  tee 
cross-section  for  collisions!  excitation.  Thus  It  Is  important  to  determine 
the  bjj‘s  ofceervat tonally.  Ve  return  to  this  point  later,  in  Chapters  IV 
and  V. 

In  the  following  section,  vs  formulate  the  method  of  computing  b-j 
and  discuss  the  collisional  crosr-secticr.3  employed  in  the  equations.  The 
numerical  results  are  discussed  In  section  3.  In  section  4,  vs  shall  con- 
eider  a method  of  determining  the  departure  from  thermodynaralc  equilibrium 
empirically  Iron  laboratory  measurements  of  the  hydrogen  arc  spectrum.  The 
present  theoretical  results  are  compared  to  the  results  of  this  experiment. 
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§2.  Formulation  of  the  equations  si*  analysis. 


The  relative  distribution  of  tie  a taros  over  n energy  levels  Is  de- 
fined 'ey  the  e quaiion  (1.1)  which  includes  bjj,  ss  a parameter.  Hence, 
the  expressions  Tor  fc*,  as  a function  of  1r>  Tej  and  Se  can  be  obtained 
free  a set  of  cyclic  equations  representing  the  statistically  steady 
state.  The  cyclic  equations  for  a very  dilute  radiation  field  end  negli- 
gible coilisicn&l  excitation  have  been  discussed  by  Kernel,  and.  his  asso- 
ciates (1937-1945).  Thoms  has  later  extended  this  seethed  to  the  case 
of  a pure  hydrogen  a troc sphere  illuminated  'ey  a radiation  field  of  inten- 
sity I tv)  mid  characterized  by  a kinetic  temperature  Te  > Tr,  where  Tr  is 
the  temperature  representing  the  continuum  of  the  radiation  field.  Com- 
putations were  carried  out  for  the  cases  Te  = 35 ,QCCP  and  7r  = 6COO° . 

He  included  the  collisions!  ionization  iron  the  discrete  levels  and 
collisiaral  excitation  from  the  ground  level,  but  ignored  the  excitations 
by  collision  from  and  between  the  other  levels.  Here  we  follow  the  sane 
procedure  to  compute  bjj  for  the  case  including  the  collisional  excitations 
from  all  the  levels. 


In  general,  we  employ  the  same  notations  adopted  try  Kensel  and 
Thosa&s,  except  for  a few  cases,  mentioned  specifically.  In  particular 
FV;_,  denotes  the  number  of  optical  transitions  from  state  A to  B per 
cubic  centimeter  per  second,  whereas  C»y  refers  to  collisional  excitations, 
or  ionisations.  'The  condition  for  a statistically  steady -state  is  then 
expressed  ass 


(2.1) 


where  ve  use  the  convention  for  principal  quantum  numbers  ar  ^ n 
X denotes  a certain  state  in  the  continuum  for  which  we  can  set  X = -in 
to1  extend  the  fbjpaal  expressions  for  discrete  transitions  to  the  contimnm. 
Then  ^ , and  for  hydrogen  2 = 1.  end  Cmi  represent 

the  collisional  de-excitations.  Tie  neglect  the  three-body  collisional 
recombination,  end  oaii  the  fine-structure  of  the  energy  levels. 


For  writing  the  explicit  expressions  for  the  terms  in  equation  (2.1), 
ve  use  the  following  symbols: 


-ft*  S-'e'k"  2 * _ 

3 it  tc 


3,2xsc 


12 


X, 

Yv 


_ i»k. 


*Tt 

£-  v'«. 


*cTe  # 

F,(-X)  = 


Yw  - r*  - v 


*°  -x 
e 


o(  x 


c = 


at  V 

Ai ry*,3-  k 


-2-6 


= /.  U y to 


•n 


Hi®  formulae  for  optical  transitions  are  given  by  Kernel  and  Ms 
associates  (1937  -1945).  For  the  spontaneous  transitions  frun  the  dis- 
crete level,  n",  or  the  continwss,  K , to  the  level  n,  we  have 


F,.„  - K 


Nc  It  i _ t 

T,! 


k „ 


-e 


3 h*5  i 


^♦1 
n / 


/ F„, 


For  the  radiation  Induced  transitions,  ve  have 


! eK'zhJ  1-'  t 


I w'*1 


- j j > 


; ! ft3  # 


>l  * 


I 


f^-K^i'sukU^r.)) 

. 

i 


(2.2) 

m 

(2.3) 


(2.4) 


(2.5) 


where  dW  de-notes  the  rr*unt  factor,  which  is  tabulated  by  Wenzel  (1937- 
1945).  Dr.e  quantity  is  a factor  intended  to  account  fer  the  presence 
of  absorption  lines  or  any  ether  affects,  such  as  dilution,  which  cause 
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the  radiation  to  depart  from  a true  Planck  distribution.  In  order  to 
perform  the  integration  over  the  range  of  frequencies  higher  than  , 
for  the  transitions  from  or  to  the  continuum,  we  use  the  averaged  quan- 
tities ^ In  the  actual  computation,  we  set  = 1,  further,  we 

note  that  equations  (2.4)  and  (2.5)  include  the  radiation  induced  emissions. 


fc  ve  proceed  to  consider  the  collisiocai  tern,  As  discussed 

in  detail  by  I homes  (1943  s),  ve  stay  generally  neglect  coliisional  excitation 
by  stems  relative  to  that-  by  electrons.  I her.  ve  can  write  the  number  cf 
collisions  per  second  per  or.-5'  in  which  the  a hex  is  raised  from  level  n 
to  n": 


C. 


jrs 


v?Q 

T'Yi 


_'#*V 

3 Zt 


d'J , 

(2.6) 


where  V is  the  electron  velocity  which  is  assumed  to  follow  a Kaxwellian 
distribution  at  the  temperature  ie.  Qw, denotes  the  cross-section  to 
excite  the  atom  from  level  n to  r.".  iT»  corresponds  to*  the  lowest  energy 
at  which  Q is  not  aero.  Recalling  the  definition  of  the  'fc>iv  » ’**’6  see 
that  the  factor  in  the  bracket  enters  because  of  coilisional  de- 

excitations.  After  completely  eliminating  and  representing  by 

jC  , ve  have: 


c 


Vitn” 


„ i 4k 

7C0..  - — r 

8 


M he  '-4- 


k('-fe)4e 


EQ^tAEM) 


Me  proceed  to  consider  the  value  of  Q 


Hln 


At  the  present  stage,  the  theoretical  calculations  of  Qnn'&re  rather 
obscure,  the  lorn  first-order  approximation  for  excitation  from  the 
ground  stats  holds  only  for  the  collision  by  electrons  of  energies  higher 
than  ICG  eT.  Thomas  has  tried  to  esuinate  a more  reliable  cross-section 
for  low  energy  excitation  from  the  .ground  state  only.  Giovanelli  (1943a,  b} 
has  suggested  that  era  assumption  setting  the  value  of  proportional 

to*  the  corresponding  transition  probabilities,  vs-l^-  for  colli- 

sions by  electrons  whose  energies  are  slightly  above  the  corresponding 
excitation  potential.  Ke  further  assumed  that,  for  high-energy  electrons, 
the  value  of  On*-  varies  as  £"  and  derived  the  general  expression  for 

n-  es  a function  of  E.  Ifce  integration  in  equation  (2.6)  was  approxi- 
mated by'  separating  the  problem  into  tvo  cases,  J <s;  E /%'J^  end 


I £/<Te 


Is  the  present  case  the  values  cf 


X 


Jit* 


are 


very  close  to*  unity  for  the  lysian  and  Rainer  series,  therefore,  we  assume 
that  Q >„*-  is  constant,  and  set  it  proportional  to  B „„•*>  throughout 
the  whole  range  of  possible  energies  of  electrons.  Sheri  ve  write. 
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(2.8) 


The  vaxue  of  the  proportionality  constant,^,  may  be  evaluated  by  intro- 
ducing the  appropriate  value  for  certain  definite  transitions,  say  Qlx. 
Thus 


1,01  * to1"'  . 


(2.9) 


As  a first  approximation  we  assume  that  Q (measured  in  units  of  *rr£l,A'  ) 
is  equal  to  1.  Substituting  the  numerical  values,  equation  (2.8)  can  -be 
expressed  as 


-io 


Q ~2,Z5*\0  — ■ 


1 (tr-k*y  • 


(2.10) 


where  the  spontaneous  transition  probabilities  fi\  are  given  by  Menzel 
and  Pekeris  (1935). 

Integrating  equation  (2.7)  and  arranging  the  constant  factors,  we 
obtain, 


— (2.11) 


'■yf  V’V  ) 

or  substituting  the  constant  C&  , which  represents  the  terms  shov/n  on 
page  12,  we  write 


AL  .IV 

' V 7\*t  » 


(2.12) 


Now  we  consider  the  number  of  collisional  ionizations  by  electrons, 
C-nl  , where  we  use  the  suffix  i instead  of  X-  to  denots  the  integration 
of  over  transitions  to  the  whole  range  of  continuum  energies  rather 

than  to  a certain  state,  x.  , in  the  continuum.  We  have  an  expression  of 
0 • similar  to  equation  (2.10)  as  follows: 

= k(.-tKM  Q„;  . 


(2.13) 
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i'or  the  value  of  Q.-y*f  , Giovanelli  (1948a)  has  also  adopted  a theoretical 
formula.  On  the  other  hand,  Thomas  gives  an  empirical  expression  for 
evaluated  from  laboratory  experimental  data.  We  may  compare  the 
values  from  these  expressions  v/ith  experimental  data  in  the  following 
table. 


Table  I. 


20  eV 

25  eV 

30  eV 

50  eV 

100  eV 

Giovanelli 

1.79 

2.03 

2.03 

1.61 

0.96 

Thomas 

0.18 

0.32 

0.45 

1.00 

0.95 

Experimental 

0.4 

0.6 

0.8 

1.2 

1.2 

* Mott  and  Massey.  Theory  of  Atomic  Collision  (2nd  Ed:  Oxford:  Clarendon 
Press  (1949).)  p.  244. 


Thus  we  see  that  for  the  more  important,  low  energy  values,  Thomas^  empiri- 
cal formula  agrees  better  with  experiments.  Therefore,  we  adopt  this 
empirical  cross-section  for  ionization  from  the  ground  level,  and  for  the 
ionizations  from  the  higher  levels  we  again  assume  proportionality  to  the 
optical  transition  probabilities,  B*;  . Thus, 


B.i 


- 0,o<H 


B,, 


(2.14) 


after  substituting  the  integrated  value  of 
R may  be  expressed  as, 


Q\/l  of  Thomas'  formula  (1948b). 


Q « *1  24  Q 

Sv  £ 


n &1- 


(2.15) 


where  the  oscillator  strength, 


is  given  by  the  formula. 


I *m  C.  3 a m a 6 x 


in  which  is  the  Gaunt  factor,  tabulated  by  Menzel  and  Pekeris 

For  our  present  calculations  we  set  9^  equal  to  unity. 


(2.16) 

(1935), 


After  integrating  equation  (2.15)  and  substituting  it  in  equation 
(2.14),  ve  have, 
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B, 


3itzx<~(>  y\.  __  ~ , -// 

,3  o3  — 3.6^7  x | o n . 


*l  3f3t,5  3'R3 

Thus,  we  finally  obtain  from  equation  (2.13) 

Qvv  = °-09l?  n ■ 

Substituting  the  above  in  the  equation  (2.12),  we  obtain 

Niti 


C-ril  C0  rp. 


L(o  ~n 

or,  after  evaluating  the  constant  terms,  it  becomes 


Cya  ~ 2.13 


3 ^ 'oiT  -^r1  K 0 - jjn'  <X+  0 . 


(2.17) 


(2.18) 


kO'-b^Cxn+D  Q* , iza9) 

^ * 


(2.20) 


Now  we  substitute  equations  (2.2),  (2.3),  (2.U),  (2.5),  (2.12)  and  (2,20) 
into  the  cyclic  equations  (2,1)  and  have,  finally. 


(2.22) 
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N^CX^i)  j. ^ , (2#23) 

\ y(x  ri*-J 

and 

7 . - 3.5  * lO^TeW*  eX"  (X,  -H  ) • <*•*> 

ht 

For  computing  purposes,  wo  may  rewrite  these  equations  in  the  following 
forms 


(2.25) 


VJo  obtain  the  values  of  bn  by  solving  the  above  system  of  equations  by 
successive  approximations. 


3«  Numerical  results  of  the  values , bn,  and  discussion. 

We  3ee  in  the  equation  (2.25)  that  the  values, bri,  are  given  for 
four  parameters,  rIe  Ne  , and  5»>.  Since  the  factor  % drops  oat 

in  the  final  equation,  the  usual  assumption  N-i  “ N0  is  not  necessary. 
The  computations  were  carried  out  under  the  same  assumption  as  Thomas', 
That  is,  we  assume  a constant  radiation  temperature  of  6000°,  a kinetic 
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temperature  of  35,000°  and  an  electron  density  N»  = 2 x 1C^,  (Thomas 
put  the  Boltzmann  constant  k equal  to  1.4  x 10~T°}  however,  we  take  the 
more  exact  value  1.37  x 10”l6>  for  the  present  computation.) 

i).  Results  neglecting  the  chromospheric  radiation  field. 

Following  Thomas,  we  consider  the  cases  » / and  « 0 • 

Our  numerical  results,  along  with  those  obtained  by  Thomas,  are  tabu- 
lated in  Table  II. 


Table  II. 


bn  for  Bo  Chromospheric  Radiation  Field. 


s 

bl 

b2 

b3 

b4 

b5 

b6 

b7 

h8 

b9 

b10 

(A)  1. 

1.6lxl06 

2.1x10°* 

9. 03 

3.32 

2.80 

2.31 

2.01 

1.75 

1.54 

1.38 

1.20 

10.9* 

4.67* 

3.45* 

2.9? 

2.62* 

2.42* 

2.27* 

2.17* 

1.06* 

(B)  0. 

1.85xl06 

10.3 

4.41 

3.24 

2.72 

2.43 

2.17 

1.90 

1.61 

1.43 

2.3x10°* 

12.3* 

5.33* 

3.92* 

3.31* 

2.96* 

2.75* 

2.60* 

2,47* 

2.36* 

* Denotes  Thomas  solution,  Paper  II. 


We  see  from  Table  I that  there  are  no  great  differences  in  the 
numerical  results  of  bn  between  the  two  calculations.  The  higher- 
level  terms  are  altei'ed  the  most,  relatively,  as  is  to  be  expected. 
Any  great  differences,  if  they  exist,  must  arise  in  the  next  case 
where  the  effect  of  the  chromospheric  radiation  field  is  included. 

ii).  Consideration  cf  the  chromospheric  radiation  field. 

If  we  add  the  chromospheric  absorption  and  emission  to  the 
radiation  field  from  the  photosphere,  £■„  may  or  may  not  be  greater 
than  unity.  For  the  lyman  lines,  whose  chromospheric  opacity  is 
high,  this  consideration  is  particularly  important.  Thomas  has  dis- 
cussed the  problem  and  has  shown  that  the  essential  perturbation  to 
the  radiation  field  by  the  chromosphere  is  in  the  Lyman  region,  with 
a small  effect  also  in  the  center  of  the  lower  order  Balmer  lines. 

He  obtained  in  his  paper,  for  in  the  lyman  lines, 


In  a private  communication,  however,  Thomas  has  mentioned  that  equa- 
tion (3.1)  is  in  error  because  of  the  faulty  inclusion  of  the  induced 
emission.  The  correct  expression,  which  follows  from  the  fifth 
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paper  of  Thomas'  series  (1949b)  under  conditions  of  large  values  for 
j where  y3  is  the  atmospheric  logarithmic  density  gradient, is 


(3.2) 


With  this  expression  for  , the  bracket  in  equation  (2.25)  involv- 
ing the  Lyman  line  radiative  transitions  vanishes.  In  the  following 
solutions  v/e  neglect  the  chromospheric  emission  outside  the  Lyman 
lines.  As  in  the  earlier  results  by  Thomas,  the  significant  para- 
meter is  ^5“  , the  residual  intensity  in  the  Balmer  lines.  Therefore 

we  have  shown  the  solution  for  bn  in  Table  III  for  two  extreme  values 
of  iS  . 


Table  III. 

tv,  including  the  Chromospheric  Emission. 

by  b2  b3  b4  bg  b7 

(G)  nS=l(l(n).  4.24x10+  1.18x103  57.4  19.0  11.2  8. 14  6.54 

(D)  2.43x105  2.36x105  71.9  26.6  15.9  11.34  3.60 

nS=l(2\n) 


t’3  b9  b1Q 
5.57  4.87  4.41 
6.75  5.53  4.59 


From  the  results  in  Table  III  ve  see  that  bn  has  a larger  value  for 
high  n than  in  the  case  without  chromo spheric  emission.  The  differ- 
ence of  b2  in  the  two  cases  is  especially  great,  bg  has  almost  the 
same  value  as  by.  We  also  find  considerable  difference  betv/een  the 
results  for  the  case  n8=l,  n>l  and  the  case  23=0,  nS=l,  ipH. 

It  may  be  interesting  to  check  these  results  by  comparing  the  value 
of  bo  with  that  which  v/e  can  compute  from  the  population  of  the  Balmer 
ground  state,  N2,  which  Thomas  obtained!  from  the  Balmer  decrement  in  the 
1932  flash  ppectruip.of  Menzel.  Thomas  obtained  N2  = 2.4  x 10^6  at  670  km 
and  N2  » 0.52  x lG-1*3  at  1500  km  where  II2  is  the  total  number  of  atomo  in 
the  second  energy  level  in  a 1 cm2  column  of  chromosphere  along  the  line 
of'  sight.  (These  numbers  are  corrected  by  a factor  of  2 for  a reason 
discussed  in  the  following  section.)  We  may  evaluate  the  values  of 
from  these  results,  assuming  an  isothermal  chromosphere  at  Te  = 35,000° 
and  Ne  = 2 x 10H.  We  obtain; 

at  h - 670  kmj  b2  = 1.15  x 10 5 

at  h = 1500  kmj  b2  = 2.5  x 10^. 

These  results  agrer  well  with  those  in  Table  HI  if  we  note  that  the  most 
likely  valu  / of  2^  is  somev/hat  less  than  1.  While  ve  must  note  that 
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Thomas'  values  of  N2  are  open  to  question,  as  will  be  discussed  in  the 
following  chapter,  the  above  agreement  seems  rather  interesting. 


§4.  Reduction  of  bfo  from  the  laboratory  experiments. 

The  theoretical  value  of  bn  discussed  in  the  preceding  sections  may 
be  compared  indirectly  with  laboratory  experiments  on  the  hydrogen  emission 
spectra  taken  under  conditions  in  which  thermodynamic  equilibrium  does  not 
hold.  Unfortunately,  only  a small  amount  of  work  has  been  done  for  differ- 
ent excitation  conditions. 


H.  Edels  and  J.  D.  Craggs  (1951)  have  investigated  the  hydrogen  arc 
spectrum  excited  by  the  electric  discharge  to  determine  the  excitation 
temperature  from  the  relative  intensities  of  H*,  and  Hv  lines.  They 
measured  the  ratios  of  the  total  intensities,  !</,  and  from 

the  line  profiles  of  H^,  and  H^.  The  experiments  were  performed  under 
various  conditions  of  gas  pressure  and  electric  current,  Th9  effects  of 
self-absorption  on  the  line-profiles  were  removed* 


The  emission  line  intensities  due  to  the  spontaneous  transition, 
X — •,  ^ , can  be  expressed  by, 


t K 


* a 


(4.1) 


where  i?,*.  is  the  total  intensity  corrected  for  the  self-absorption,  and 
Nn  is  the  total  number  of  atoms  in  the  nth  level  along  the  line-of-Bight. 
If  thermodynamic  equilibrium  is  assumed,  Nn  may  be  written  in  terms  of 
1J2  as, 


(4.2) 


where  % represents  the  excitation  temperature.  From  (4.1)  and  (4.2), 
we  have 
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(4.3) 


From  the  above  relation,  Mels  and  Craggs  determined  T«  using 
their  measurements  on  Balmer  line  arcs  taken  at  1 to  2 atmospheres 
pressure  and  carrying  currents  varying  from  about  3 to  10  amperes.  They 
found  that  % determined  from  differs  considerably  from,  and 

is  greater  than,  the  values  obtained  from  ,^/j  . Apart  from  the 

experimental  errors,  the  excitation  temperature  should  be  identical  for 
every,  pair  of  the  lines  in  the  state  of  thermodynamical  equilibrium. 
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Therefore,  it  is  natural  to  consider  that  finding  non-identical  tempera- 
tures nay  be  due  to  the  effect  of  the  departure  of  the  experimental  con- 
dition from  thermodynamical  equilibrium.. 

We  will  consider  this  effect  in  terms  of  t^.  from  our  theory,  the 
alternative  expression  of  equation  (4.3)  may  be  given  as 


I, 
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Since  bj^  is  defined  by  equation  (l.l),  the  temperature  appearing  in  the 
above  equation  refers  to  the  electi'on  temperature.  (In  thermodynamical 
equilibrium,  however,  only  a unique  temperature  should  exist,  and  Te  =■  Te^ . ) 
Taking  the  logarithm  of  the  above  equation  v;e  have 
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Substituting  the  numerical  values,  we  write, 
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(4.7) 


hj- 

TVji  and  T^jjr  are  conventional  notations  for  temperatures  derived  from 
and  fP/jj  i respectively.  The  factors  bn  are  absent  in  the 
equaxion3  for  tnermo dynamic  equilibrium,  so  that  bn  may  be  considered  as 
the  correction  factor  to  give  the  equality  T«a  = Ta?  • Thus,  equating 
(4.6)  to  (4.7),  we  have  the  relation  ' ‘ 
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Table  IV  shown  the  results  of  computation  of  the  right  hand  side  of  the 
above  equation,  using  the  values  of  Ia/Ja  and  Ineasur'£'^  by  Edels 

and  Craggs.  It  is  interesting  to  seethdx  there  are  no  great  systematic 
changes  in  the  last  column  for  different  conditions  of  pressure  and  elec- 
tric current.  The  values  .in  the  last  column  of  Table  IV  may  be  compared 
with  corresponding  quantities  computed  from  the  theoretical  data  in  Tables 
II  and  III.  We  have  the  values  corresponding  to  the  left  hand  side  of 
equation  (4.8)  as; 

(A)  -0.10  (B)  -0.08  (C)  -0.04  fl))  -0.12 

where  (A),  (B),  (C),  and  (D)  correspond  to  each  row  designated  by  the 
same  letter  in  Table  II  and  Table  II.  We  see  that  the  value  varies 
slightly  for  the  different  conditions  assumed  for  computations.  If  we 
consider  that  the  physical  situation  postulated  in  the  theoretical  com- 
putation ic  considerably  different  from  that  under  which  each  experiment 
ip  performed,  the  theoretical  values  may  be  considered  harmonious  with  the 
experimental  results.  We  may  be  able  to  conclude  at  least  that  the  orda'a 
of  the  relative  values  of  bjj  are  reasonable  as  a first  approximation. 


Table  IV. 


Arc.  No. 

Current 

(amp) 

V 

T 

Dif  br 

& 

1 

5.7 

1.84 

1.32 

-0.40 

2 

6.0 

1.76 

1.37 

-0.59 

s 

3 

7.3 

1.65 

1.27 

-0.48 

4 

9.4 

1.50 

1.26 

-0.61 

CO 

5 

9.9 

1.52 

1.19 

-0.44 

6 

3.35 

9.95 

4.54 

-0.35 

& 

7. 

3.53 

9.92 

4.58 

-0.38 

8 

5.2 

9.97 

3.73 

+0.06 

g 

9 

5.4 

5.85 

3.47 

-0.30 

o 

10 

8.6 

4.22 

3.50 

-0.65 

3 

11 

8.68 

4.14 

3.46 

-0.65 

12 

8.80 

— 

3.26 

X 

13 

3.70 

5.89 

3,73 

-0.47 

6 o 

14 

4.05 

5.33 

3.58 

-0.51 

m 

15 

5.1 

4.94 

3.67 

-0.60 

0 

16 

5.9 

4.52 

3.44 

-0.56 

17 

6.0 

4.51 

3.31 

-0.47 

a 

18 

8.3 

3.74 

3.46 

-0.75 

19 

8.9 

3.53 

3.66 

-0.94 

Since  we  are  comparing  the  relative  values  of  bn  rather  than  the  absolute, 
the  above  agreement  is  significant  even  though  the  experimental  conditions 
are  very  different  from  those  assumed  in  the  computation. 


III.  THEORETICAL  CONSIDERATION 


OF  THE  SELF-ABSORPTION  IN  THE  STEI-1 A a ATECSP'HFEE 


§1.  Introduction  and  assumptions; . 

Despite  the  importance  of  the  possible  effect  of  self-absorption 
in  the  chromosphere,  little  attempt  has  been  mude  to  consider  this  effect 
in  the  Interpretation  of  eclipse  spectra.  The  chief  difficulty  for  a 
theoretical  study  of  self-absorption  seems -to  be  the  lack  of  information- 
on  the  distribution  of  the  density  of  atoms  and  electrons,  as  well,  as  on 
the  temperature  gradient.  Further,  for  stronger  lines  where  the  self- 
reversal appears  predominant,  one  meets  the  more  fundamental  problems 
related  to  the  details  of  line-formation. 

A quantitative  analysis  of  the  effect  of  self-aosorption  in  the 
Balmer  lines  was  first  carried  out  by  Thomas  (1950b)  to  interpret  the 
Balmer  decrement  for  the  higher  series  members,  as  will  be  discussed  in 
more  detail  in  the  following  chapter.  Thomas  derived  an  expression  for 
the  total  amount  of  self-absorption  in  the  Balmer  line  observed  with  a 
slitless  spectrograph.  He  assumed  an  exponential  decrease  of  hydrogen 
density  with  height  and  a constant  temperature  throughout  the  atmosphere. 
In  the  analysis  of  the  1932  eclipse,  Cillie  and  Kenzcl  (1935)  gave  a 
similar  s eries  formula,  obtained  from  the  same  assumptions;  but  as  they 
considered  the  emission  from  unit  cross-section  of  the  chromosphere,  the 
convergence  of  uhe  series  was  too  slov>  to  perform  a numerical  computation 
Hence  this  method  was  not  practical  for  eclipse  analysis.  The  series 
treated  by  Thomas  converges  considerably  more  rapidly,  but  he  uses  a 
cumbersome  calculation  involving  double  integrals. 

In  this  chapter,  we  sha!3  reconsider  the  mathematical  formulation,  of 
the  problem,  and  a considerably  simpler  expression  than  Thomas'  will  be 
given.  Further,  Thomas  took  half  of  the  t' tal  population  of  the  Balmer 
ground  state  in  place  of  the  total  number  in  a unit  column  along  the  line 
of-sight,  and  the  convergence  of  his  series  is  much  slower  than  in  the 
series  using  the  total  number  as  shown  in  the  following  derivation. 

(The  values  of  N2  given  by  Thomas  in  Table  I]  are,  as  he  stated,  half  of 
the  total  number  instead  of  the  total  number,) 

As  we  shall  see  later,  the  self-absorption  for  the  Balmer  emission 
depends  only  on  the  population  of  the  hydrogen  second  state,  H2>-  and  the 
atomic  absorption  coefficient,  C-( . The  temperature  enters  in  the 
latter  only  if  we  use  the  absorption  coefficient  at  the  line  center. 

Thus,  it  is  important  to  note  that  the  effect  of  temperature  and  density 
may  be  treated  independently,  particularly  if  we  want  to  repeat  the 
method  of  numerical  integration  by  successive  approximations.  Hence  in 
the  following  discussion,  we  first  assume  an  exponential  decrease  in 
density  along  the  radial  direction,  then  set  the  ebserption  coefficient 
constant.  Cur  results,  for  the  weaker  lines , are  not  toe  sensitive  to 
the  latter  assumption  because  most  of  the  emission  and .bsorption  origi- 
nates at  the  small  central  region  of  the  moon's  limb  so  that  the 
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temperature  may  be  essentially  constant,  fciince  for  the  stronger  lines, 
we  have  to  take  the  failure  of  these  assumptions  into  account,  the  direct 
expansion  of  the  absorption  formula  may  not  be  applicable  for  these  cases. 
Further,  we  note  that  the  parameters,  N2  and  enter  as  a product,  so 

that  ve  can  always  adjust  the  actual  value  of  the  temperature  in  the  final 
computation,  as  will  be  shown  later. 

In  the  following  section,  we  will  develop  the  formula  for  the  observed 
intensity  and  the  self-absorption  at  a certain  point  within  a line,  then 
integrate  over  a line  contour  to  get  the  total  amount  of  energy  and  self- 
absorption (section  3).  The  resulting  formula  is  much  simpler  than  the 
one  obtained  by  Thomas.  The  numerical  values  of  the  self-absorption  are 
discussed  in  the  last  section  and  are  shown  graphically  as  a function  of 
Jtoq.  . Then,  from  this  curve,  \ie  derive  some  analytical  expressions 

which  agree  almost  exactly  with  the  numerical  computation. 


FIGURE  I Geometrical  Diagram 


§2.  Formulation  of  the  self-absorption  formula  for  monochromatic  light. 


Figure  I shows  a diagram  of  the  sun's  limb,  indicating  that  we  take 
the  y-coordinate  along  the  linc-of-sight  and  the  x-  and  z-axis  perpendicular 


to  the  y-axis.  Let  us  consider  i.,o  energy  per  unit  solid  angle  emitted 
vithin  a frequency  int  erval  to  V t «ti>  by  transitions  from  the  nth  to 
the  second  level  of  the  hydrogen  atoms  in  a small  volume  element 
at  point  (x,  y,  z).  Thus, 
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(2.1) 


where  nn  denotes  the  manber  of  hydrogen  atoms  in  the  nth  level  within  a 
unit  volume,  and  p\nzis  the.  Einstein  coefficient  for  the  transition  prob- 
ability. (Throughout  the  remainder  of  this  paper  we  shall  use  small  latter 
for  the  number  of  particles  per  unit  volume,  and  capital  letters  for 
quantities  integrated  over  a certain  volume. ) $ refers  to  the  atomic 

absorption  coefficient  within  the  frequency  interval  \>  and  ptiP  . i/e 
use  the  suffix  2X;)  for  E and  oC  to  designate  these  quantities  for  a 
small  region  at  the  frequency  p within  a line  and  the  suffix  jj  r.  for  the 
total  quantities  integrated  over  the  line.  Thus  we  write, 

to 
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?be  other  notations  are  conventional. 


It  ”ill  be  easiest  if  we  consider  the  radiation  from  a slab  of  width 
1 'i  in  the  z~ direction,  so  that  we  can  eliminate  the  z-component  by 
operating  with  J - on  equation  (2.1).  Thus,  from  equation  (2.1 ) 

and  (2.2)  we  get  - 
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Vie  will  nov;  consider  the  effect  of  self-absorption  on  the  above 
radiation.  Since  all  hydrogen  atoms  in  the  second  state  from  the  point 
(x,  y)  to  the  observer  contribute  -bo  the  absorption,  the  amount  of  energy 
that  finally  emerges  from  the  solar  atmosphere  will  be 
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where  rg>  (x,  y)  is  the  population  of  hydrogen  atoms  in  the  second  level 
at  the  point  (x,  y).  V/c  introduce  the  notation 
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Then  equation 


(2.4)  becomes 
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Prom  the  photometric  measurements  of  the  flash  spectrum  taken  through 
a slitless  spectrograph,  ue  obtain  the  total  energy  emitted  from  the 
volume  of  the  chromosphere  in  a 1 cm  slab  above  a certain  height,  h.  There- 
fore, we  integrate  equation  (2.6)  along  the  x~  and  y-axis, 
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here  (x,  y)  may  be  expressed  as  a function  of  bn,  T0  and  Ne  according 
to  the  Boltzmann-Saha  formula  which  includes  the  departure  from  thermodyna- 
mic equilibrium: 
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Thus,  we  see  that  the  integration  of  equation  (2.7)  is  difficult, 
for  it  depends  on  the  distributions  of  !he  pirameters,  bn,  Te,  nj,  ana  ne 
throughout  the  atmosphere.  If  \;c  know  these  quantities  as  functions  of 
x,  y,  and  z,  we  may  at  least  perform  the  integration  numerically.  Since, 
however,  these  parameters  enter  implicitly  in  equation  (2.7),  one  method 
is  to  assume  some  analytical  form  for  the  distribution  of  the  atoms  instead 
of  assuming  each  parameter.  7'or  the  present  case,  we  introduce  the  assump- 
tion of  spherical  symmetry  fox'  the  distribution  of  hydrogen  atoms,  and 
assume  an  exponential  decrease  in  density  along  the  radial  direction,  r. 
Thus,  we  write 


(2.9) 


where  the  constant  nn  (0,  0)  represents  the  number  of  atoms  in  the  nth  level 
at  the  base  of  the  chroma son ere.  In  general,  the  density  gradient,  £ also 
should  be  considered  to  vary  for  different  r egiens  in  the  atmosphere.  While 
we  set  £ constant  in  the  following  development,  wc  may  have  to  consider 
its  variation  with  height  in  the  numerical  integration.  According  to  the 
usual  geometric  approximation,  we  set 
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Substituting  equation  (2.10)  into  (2.1l),  we  have 
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Or,  recalling  our  notation  given  in  (?.!)),  we  may  write  (2.1l)  a,s 
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Let  U3  now  consider  the  double  integral 
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We  will  first  consider  the  integration  over  /,  then  over  x.  Expanding 
the  last  exponential  term;,  wc  write 
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where  N2  (x,  y)  may  be  written  as 
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and 
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Hence,  for  y 0,  wc  have 
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Introducing  a notation 
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we  may  write  equation  (2.16)  as 
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Here  we  see  that  the  quantity,  71*. C0,0)  e ^ ,j^p  represents  the  total  number 
of  hydrogen  atoms  in  the  second  level  within  a column  of  unit  cross- 
section  along  the  line.-of-sight  at  a height  x above  the  base  of  the 
chromosphere.  We  denote  this  number  by  Nq(x)  or  simply  N2,  so  that 


Thus,  substituting  (2.19)  into  (£.18),  we  have 
fory>0, 

By  a similar  method,  we  readily  obtain 
for  y <,  0, 
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For  substitution  in  the  power  series  (2,14),  we  calculate, 
for  y)>  0, 
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and,  for  y 0, 
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Now  returning  to  the  integral  (2.13),  we  first  write  the  y-dependent 
part  of  the  integral  as  follows: 
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From  equation  (2.14)  we  may  write, 
for  y } 0, 
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and,  for  y<(  0, 
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Hence,  equation  (2.24)  may  be  written 
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Again  integrating  over  y,  we  have 
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where  ve  introduce  a notation  3^: 


J*  " -k  f*F  rcr']  ap  • 
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Now  equation  (2.28)  involves  a series  of  Im,  which  requires  a tedious 
computation  if  tne  convergence  of  the  series  is  slow.  We  can,  however, 


simplify  the  equation  by  methods  fully  developed  in  Appendix  IIIA,  (see 
page  37  ).  Namely,  we  can  write; 
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Substituting  the  above  into  equation  (2.28),  v;e  have 
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Returning  to  the  integral, , given  in  (2.13),  we  multiply  both 
sides  of  equation  (2.31)  by  a factor  q~  P**,  and  integrate  over  x from  h 
to  infinity.  Thus,  recalling  our  definition  (2.10), 

Hx  00'  = N (o)  e (2.32) 


we  finally  obtain 
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We  again  note  that  N2(h)  is  the  total  number  of  hydrogen  atoms  in  the 
second  level  in  a 1 cm2  column  along  the  line-of-sight  at  a height,  h. 

wow,  from  equation  (2.10),  we  finally  obtain  the  equation  for  the 
observed  Balmer  line  emission  within  the  frequency  interval  i)  and  •p  + U'p 
in  the  slitless  flash  spectrograms  in  the  form: 
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or,  using  the  relation  (2.19),  >r\ 
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The  sex'ies  representing  the  last  factor  of  the  above  equation  corres- 
ponds to  the  fraction  of  energy  absorbed  in  the  atmosphere,  which  we  here- 
after denote  by  a symbol,  |Aksnj  • Namely, 


(Na^n.,)?)  , ( ^2  tj) 
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(2.36) 


Returning  to  the  integration  (2.13),  we  recognize  that  0^2’n;>>  has  been 
implicitly  assumed  to  be  constant  to  perform  the  integration  over  the 
atmosphere.  According  to  our  definition,  oC^  p refers  to  the  absorption 
coefficient  at  a frequency,  -p  , so  that  it  depends  on  the  atmospheric 
conditions  determining  the  line-profile.  The  line-profile  may  be  deter- 
mined by  various  effects,  and  p becomes  a function  of  various  para- 

meters. Therefore,  in  a rigorous  sense,  a numerical  integration  only  may 
be  possible  for  a given  model  of  atmosphere.  In  the  following  section, 
we  shall  integrate  equation  (2.35)  over  frequency  to  derive  the  total 
emission  of  the  line.  As  we  shall  discuss  later,  we  may  assume  Doppler 
broadening  to  be  predeminent  for  chromospheric  conditions,  for  which 
cX-zri  is  a function  of  Te  only.  Hence,  only  the  isothermal  atmos- 
phere Wist  be  assumed  in  performing  the  integration  (2.13).  We  shall  re- 
turn to  this  point  later. 


§3 . Consideration  of  total  line  emission  and  self-absorption. 

In  the  previous  section,  we  have  considered  the  monochromatic  emission 
corrected  for  self-absorption  in  the  solar  atmosphere.  For  the  analysis 
of  the  eclipse  spectrum,  it  may  be  desirable  to  have  the  corresponding 
expression  for  the  total  amount  of  energy  involved  in  a line.  We,  there- 
fore, integrate  equation  (2.35)  over  the  line.  While  the  line  contour 
may  be  determined  by  various  effects,  we  recognize  that  Doppler  broaden- 
ing is  predomlnent  for  the  chromospheric  conditions  we  assumed,  and  neglect 
the  terms  due  to  collision  damping  and  Stark-effect . Further,  if  we  con- 
sider the  formation  of  the  detailed  line  structure,  we  have  to  include 
the  dependence  of  the  emission  at  the  different  frequencies  on  the  distri- 
butions of  various  parameters.  However,  we  are  considering  only  the 
higher  members  of  the  Balmer  series  for  which  optical  depths  are  small 
and  we  may  consider  that  most  of  the  emission  comes  from  a small  region 
of  the  chromosphere,  as  we  have  discussed  earlier.  For  stronger  lines, 
this  simplification  may  not  be  justified. 


Thus,  we  may  approximate  the  a bsoi’ptlon  coefficient  as  follows: 
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where 
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Then,  from  equation  (2.35),  we  have 
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The  integration  of  the  nth  term  of  the  series  will  be 
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represents  the  absorption  coefficient  at  the  center  of  the  line.  Thus, 
we  finally  obtain 

iTt 
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for  the  total  emission  in  the  nth  Balrner  line  in  the  volume  of  chromo- 
sphere considered.  Or  we  may  write 


The  above  series  is  absolutely  convergent,  but  for  N^a(2n^0  larger 
than  10  the  convergence  becomes  very  slow.  For  example,  to  obtain  three 
decimal  places  for  CfiW}  > we  have  to  carry  out  six  figure  computations 
of  27  terms.  The  amount  of  computation  increases  very  rapidly  for 


we 


exceeding  10.  Comparing  (3.8)  to  Thomas'  equation  (I950u),  however, 
sec  that  the  form  of  the  series  is  considerably  simplified,  and  the 
coefficient  of  t)^  is  fl  factor  °f  2.*^  as  large.  Thu3  the  range 

of  calculations  is  extended  widely. 

As  ve  have  mentioned,  following  equation  (3.5),  the  absorption 
coefficient  at  th  line-center,  o^-no  i£  proportional  to  t;'1  , hence  is 
a slowly  varying  function  of  the  kinetic  temperature , TeJ  — • 

Since  we  have -to  assume  an  isothermal  chromosphere  to  perform'' the  inte- 
gration developed  in  the  preceding  section,  equation  (3.7)  may  not  be 
adequate  for  stronger  lines  (lower  members  of  Calmer  series)  where  the 
optical  depths  are-  so  high  that  the  outer  layers  of  the  chromosphere  con- 
tribute considerably  to  the  absorption.  For  the  higher  series  linos,  how- 
ever, most  of  the  emission  originates  in  the  small  central  r egion  at  the 
lowest  layer  observed.  Hence,  we  inay  use  this  method  as  a first  approxi- 
mation and  a numerical  integration  must  be  perform- d for  a detailed  analysis. 


§4 . Kuroerical  results  and  the  approximate  formula  for  self-absorption. 

In  the  preceding  section  we  have  devt loped  the  theoretical  expression 
for  self-absorption  as  a function  of  0 . If  the  observed  amount 

of  ael f-absorption  in  the  Calmer  flash  spectrum  line  Hn  is  determined, 
then  we  can  compute  the  value  of  from  equation  (3.8).  Hence, 

the  population  of  hydrogen  atoms  in  the  "second  level,  M2,  may  be  obtained 
if  we  have  the  value  of  the  chromospheric  kinetic  temperature  for  calculat- 
ing the  absorption  coefficient,  oOj>n  0 . Equation  (3.8)  is  absolutely 

convergent,  but  it  requires  cumbersome  computations  for  higher  values  of 
Nxo<xhjO  • After  repeating  the  computation  for  a very  snail  interval 
of  Ni*m,o  > however,  the-  author  has  found  that  the  logarithms  of  fAy$-A 
can  be  represented  as  a function  of  log  o '°y  & smooth  curve,  shown  fn 

Figure  II.  The  exact  calculation  of  the  series  has  been  Carried  to 

o ~ / 0 • As  we  have  discussed  in  the  previous  section,  we  have 
to  compute  27  terms  of  the  series,  carrying  six  figures,  to  get  three 
decimal  places  for  f At>s-n_|  . It  is  practically  impossible  to  extend  the 
computation  beyond  N^in  o“  /0  • However,  ve  recall  that  formula  (3-8) 
is  derived  from  certain  simplifying  assumptions  discussed  prior  to  its 
derivation,  so  it  may  not  be  applicable  tor  stronger  lines  (higher  values 
of  Ni  c4aTVjo  )•  Log  Os  3.;  iOj  o is  found  to  be  -14.59  and  -14.97 

for  Te  equal  to  6,000°  and  35,000°,  respectively;  and  log  o<£;)j-,  o > Ls 
-15.16  and  -15.54  for  the  same  temperatures.  Therefore,  wc  sec  that 
Figure  II  can  cover  the  lines  as  high  as  Iiio>  if  log  is  smaller  than 
16,  and  to  H15  for  log  N2  of  the  order  16.5.  From  the  results  obtained 
in  the  following  two  chapters  (cf.  equations  (6.1),  Chapter  IV,  (2.6), 

Chapter  V,  and  Figure  IV),  we  estimate  the  value  of  log  N2  to  be  of  the 
order  16  and  16.5  corres ponding  to  the  chromospheric  heights  of  1400  Ion 
and  500  km,  respectively.  That  is,  wo  can  obtain  sAbs^  for  the  lines 
of  higher  order  than  HpQ  for  the  upper  levels  and  H15  for  Ihe  lower 
levels  in  the  chromosphere. 

Vhile  we  may  thus  either  use  the  curve  of  Figure  II  to  obtain  the 
amount  of  self-absorption  for  given  N2  and  Te  or  to  obtain  the  value  cf 


3: 


N*  atfr  0 inversely  from  some  otherwise  determined  self-absorption,  it 
will  be  v-.ry  useful  if  we  can  find  some  analytic;:;]  '-xpression  representing 
it.  For  example,  if  we  calculate  N2  from  the  relative  amount  of  self- 
absorption  as  will  be  described  In  the  following  chapter,  we  nc.ed  to  have 
such  a formula,  further,  if  we  rely  on  our  observation  that  the  curve 
chang-3  its  slope  very  slowly  for  log  KJ  a 0^ xr\ , 0 higher  than  about  0.3, 

(that  is,  the  curve  appears  to  be  nearly  a straight  line)  the  formula 
thus  obtained  may  be  used  to  extrapolate  the  function  into  the  region 

/V  > I Q 

By  repeated  use  of  the  method  of  trial  and  error,  we  finally  found 
that  the  expression 

- -0,0 <1  e ■ ’ 


gave  the  beat  representation  of  the  curve,  Equation  (',.!)  may  also  be 
written 


In  Table  V the  values  computed  from  equation  (4.2)  are  compared  to  those 
read  from  various  points  in  Figure  LI. 


Table  V. 

l>gM^  -1.0  -0.7  -0.4  -0.2  0.0  0.2  0.4  0.6  0.8  0.9  1.0 

7-^VU0*008  c'016  0*030  0.045  0.070  0.106  0.157  0.224  0.303  0.353  0.414 


^^Vp.OOf  0.010  0.025  0.043  0.069  0.107  0.158 


0.225  0,303  0.345  0.389 


Thus,  we  sec-  that  equation  (4.1)  gives  the  exact  value  with  an  accuracy 
of  three  figures  except  for  the  first  three  and  the  last  two  columns. 

Since  the  lest  figure  is  insignificant  in  the  practical  analysis  of  obser- 
vations, the  expression  (4,1)  may  be  considered  to  give  the  exact  solution 
for  x^2n,o  < l0  . However,  it  is  obviously  incorrect  to  extrapolate 
to  larger  vulues  of  N^x^o  because  we  see  in  equation  (4.1)  that  £*,«  f At- 
tends to  zero  as  n goes  infinity.  That  is,  the  gradient  of 

evidently  changes 'its  sign  at  some-  value  of  *3 

obviously  without  physical  significange.  We  muy  see  the^  direction  of  this 
change  appearing  in  the  slight  deviation  of  the  values  in  the  later  part 
of  the  Table  IV,  It  vdll,  therefore,  be  best  tc  limit  the  use  of  equation 
(4.1)  to  the  analysis  of  smaller  Nx^in^o  • 


We  have  also  introduced  a second  approximate  formula  as  follows: 


a 

Aks„j  = -0.01k  (Ni ct^y  (4. 3 ) 

The  comparison  of  this  formula,  with  the  exact  value  i s shown  in  Table  VI. 


Table  VI. 


-0.B  -0.6  -c.4  -0.2  0.0  0.2  0.4  0.6  0.8  0.9  1.0 

10.012  0.C19  0.030  0.045  0.07O  0.106  0.157  0.227  0.303  0.352  0.414 
Fig. II  


10.019  0.029  0.038  0.054  0.076  0.107  0,152  0.214  0,302  0,355  0.42 

Equ . (4 .3} 


In  the  above  Table,  we  see  that  equation  (4.3)  gives  ar.  even  better 
approximation  than  equation  (4.1)  for  higher  values  of  NzP< 

Further,  for  the  whole  range  of  o the  Bimpleiv  expression  (4.3) 

has  two  decimal  place  accuracy.  Since  tint-  accuracy  of  the  present 
measurements  of  the  eclipse  spectrum  could  not  exceed  two  decimal  places 
In  the  logarithm,  equation  (4.3)  may  well  be  taken  in  the  actual  analysis. 
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IIIA.  APPENDIX 


Ills 


Consider  the  following  series  involved  in  equation  (2,23), 


5U 


+ — ^ --  L + 


i 

,2  10-2)1  r 4 I (n-f)  I ^ 


where  IK  represents  the  integral 
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From  equation  (A. 3),  we  can  write 
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Hence,  frora  equation  (A. 2),  we  have 
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Substituting  the  above  equation  into  (A.l),  we  get 
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Now  consider  the  polynomial  expansion 
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Comparing  the  above  expansions  to  series  (A.l),  v<=- readily  See 
may  be  replaced  by  an  integral 
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Hence,  integrating  (A. 9)  directly,  we  have 
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IV.  APPLICATION  TO  THE  FLASH  SPECTRUM 
OF  1932  ECLIPSE 


§1.  Introduction. 

In  Chapter  I,  it  was  pointed  out  that  the  low  decrement  of  the 
Palmer  series  lines  is  one  of  the  anomalous  features  of  the  flash  spec- 
trum that  cannot  be  explained  by  a model  of  the  chromosphere  that  is  a 
simple  extension  of  the  photosphere.  Thomas  (1950a)  has  first  considered 
the  problem  quantitatively,  and  he  suggested  that  self-absorption  is 
responsible  for  this  characteristic  of  the  spectrum.  By  comparing  the 
ratio  of  the  observed  intensities  of  a pair  of  higher  order  lines  of  the 
E aimer  series  with  the  theoretical  intensity  ratio  computed  from  the 
assumption  of  thermodynamic  equilibrium,  Thomas  found  that  the  differences 
of  these  ratios  show  a systematic  decrease  both  with  order  of  the  lines 
and  height  in  the  chromosphere.  Since  this  systematic  difference  in- 
creases with  increasing  self-absorption,  Thomas  assumed  that  the  differ- 
ence between  the  observed  and  the  computed  Balmer  decrement  was  due  to 
the  effect  of  the  self-absorption  only.  By  introducing  this  difference 
into  his  absorption  formula  (of.  equation  U),  Thomas  (1950a)),  Thomas 
computed  the  population  of  hydrogen  atoms  in  the  second  state  at  various 
heights.  In  his  computation,  however,  he  has  taken  into  account,  at  each 
height,  only  one  pair  of  very  high  order  members  of  the  Balmer  series 
because  of  the  difficulty  of  computation  due  to  the  slow  convergence  of 
the  absorption  formula  he  derived.  Thus,  the  numerical  results  for  N2 
are  open  to  question  if  there  is  evidence  of  considerable  fluctuation 
in  the  observations.  In  fact,  the  fluctuation  of  N2  values  computed 
for  different  pairs  of  lines  were  quite  large.  We  will  discuss  this 
point  in  the  next  and  last  sections. 

Further,  Thomas  has  explicitly  neglected  the  factor  bn  v;hich  should  j 

be  included  under  the  assumption  of  high  chromospheric  temperature.  Hence, 
as  he  pointed  out,  the  amount  of  absorption  that  he  estimated  is  to  be 
considered  as  a lower  limit.  In  section  2,  we  will  discuss  the  effect  of 
these  simplifications,  and  point  out  some  additional  inconsistencies  in- 
volved in  Thomas'  assumptions.  Then,  in  the  later  sections,  we  will 
consider  several  possible  methods  for  determining  a self-consistent  model 
of  the  chromosphere. 


§2.  Comment  on  Thomas'  determination  of  N?  from  self-absorption. 


As  discussed  in  the  previous  section,  we  will  first  follow  Thomas' 
method  to  check  his  results.  We  start  from  the  equation  for  Balmer  line 
intensities.  Chapter  111’,  equation  (3.7): 
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To  make  the  factor  fc-a  explicit  in  the  stove  equation,  we  replace 
the  spontaneous  transition  probability,  by  the  total  atomic  absorp- 
tion coefficient,  , according  to  the  relation 
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from  the  definition  of  bn,  Chapter  II,  equation  (1.1). 

Hence  we  can 

write 
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From  equation  (2.2)  we  note  that  the  second  term  in  the  bracket  of  equation 
(2.4)  gives  the  contribution  due  to  the  radiation  induced  emission.  Sub- 
stituting equation  (2.4)  into  (2.1),  ve  have 
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If  we  work  with  the  logarithmic  gradient  of  , we  eliminate  all  the 

factors  that  do  not  depend  on  the  quantum  number,  n.  Thus  we  have 
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where  * "ft  "T  **  ■f  C»t  k.  ) • . The  above  e quation  corresponds  to 

equation  (3a)  in  Thomas 1 paper.  Note  that  Thomas,  in  equation  (l), 
implicitly  omitted  the  radiation  induced  emission  term  by  placing 
outside  the  integration.  This  is  not  generally  justified,  even  though  as 
ve  shall  see,  it  leads  to  the  same  equation. ■ 

We  first  consider  tho  term  "jffeV  The  numerical  values 

of  -e-Xj'a  ('RJx.yi/ 1?')  are  shovn  Tn  Table  VII* 


a 

Table  VII. 
Value  of  6.  K'tZ  . 


Te  -■ 

6000 

10,000 

20,000 

35,000 

^>1° 

8.7 

• 109 

2.3  • 105 

8.7  • 10 

3.0 

^d,30 

9.1 

• 109 

2.5  • 105 

9.1  • 10 

3.1 

Thus  we  see  that  for  temperatures  lower  than  20,000°,  we  can  imme- 
diately neglect  the  induced  emission  term,  since  b*.  ^ b ^ .•  Prom  the 
results  of  the  theoretical  computations  in  Chapter  II,  we  know  that 
k/b„^\05  > lea3'i:  f°r  n higher  than  10  and  Te  = 30,000°.  Therefore 
we  can  state  that  v*  ivxp  (fv^/lcTe)^  ^f°r  the  range  of  conditions 

in  tho  nrestnt  study,  and  ve  can  ignore  the  radiation  induced  emissions . 
Accordingly,  we  vrrite 
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Thus,  we  see  that  the  log  62  term  subtracts  out  of  equation  (2.6),  leaving 

E- =*  3 (~2~  ~ •**  ^ tnr, 

To  show  the  effect  of  temperature  variation,  we  have  tabulated  in 
^abie  VIII  the  differences  of  fourth  terms  in  equation  (2.8). 


Table  VIII. 

"ft  ^2115  . yAC<H, 
kTe 


n - n 

i + k 

Te=6000° 

10,000° 

20,000° 

35,000° 

25 

- 30 

0.01 

0.00 

0.00 

0.00 

20 

- 25 

0.01 

0.01 

0.00 

0.00 

15 

- 20 

0.02 

0.01 

0.01 

0.00 

10 

- 15 

0.06 

0.04 

0.02 

0,01 

10 

- 13 

C.05 

0.03 

0.01 

0.01 

Since  the  observed  data  used  for  the  present  computation  are  given 
to  two  decimal  places  only,  ve  can  say  that,  at  least  for  n higher  than 
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15,  logarithmic  decrements  are  not  "t  all  sensitive  to  temperature  varia- 
tions. for  n between  10  and  15,  there  is  also  no  temperature  effect  for 
Te  higher  than  20,000°.  Slight  differences  due  to  the  temperature  term 
may  have  to  be  considered  for  n :£  |S"  and  TbSsIOjOOo*  • 

Thomas  has  ignored  the  bn  factor  in  equation  (2.8).  Unfortunately, 
at  the  present  time,  there  is  no  exact  computation  of  bn  for  n higher 
than  10.  However,  we  may  use  the  asymptotic  values  computed  by  Thomas 
to  estimate  the  order  of  magnitude  of  Av,  6*^.  b * . (It  should  be 

noted  that  these  values  are  obtained  from  computations  in  which  colli- 
sional  excitations  are  neglected,  Cf.  Chapter  II.).  In'  Table  IX  we 
compare  these  values  to  the  values,  at  each  height  for  which 

Thomas  computed  the  values  of  No-  The  values’*  inside  the  brackets  are 
computed  using  the  small  corrections  in  Table  VIII  for  Te  = 35,000°. 

Table  IX. 


and 

Height 

670  km 

1500  km 

2330  km 

3170  km 

n - n + k 

25  - 30 

15  - 20 

10  - 15 

10  - 13 

0.02 

0.07 

0.13 

0.09 

0.04 

0.05 

0.10(0.09) 

0.05(0,04) 

Thus,  from  the  above  table,  we  see  that  the  factors  23^ 4m are 
even  larger  than  ^ . Even  if  we  consider  that  the  above 

b-Y-1  are  upper  limits,  we  may  at  least  say  that  both  factors 
are  about  the  same  order  of  magnitude,  so  that  the  neglect  of  the  bn 
term  is  by  no  means  justified.  On  the  other  hand,  if  we  assume  the 
electron  temperature  to  be  as  low  as  6,000°,  where  bn  terms  vanish  theore- 
tically, the  exponential  term  becomes  significant  (Table  VIII).  Ihere- 
fore,  the  omission  of  both  these  terms  leads  to  much  too  low  a value  for 
N2.  We  will  return  to  this  question  later.  As  long  as  we  assume  the 
electron  temperature  to  be  of  the  order  of  35,000°  the  largest  factor  in 
determining  the  temperature  as  veil  as  its  variation  seems  to  be  the  bn 
factor  rather  than  and  the  contribution  of  the  exponential  term  ia 

negligible.  In  the  later  sections  of  his  paper  (1950a),  Thomas  introduced 
the  logarithmic  height  gradient  of  N2  obtained  from  his  computation  into 
the  equation  of  hydrostatic  equilibrium  at  Te  = 35,000°,  and  determined 
the  temperature  gradient,  ignoring  a gradient  in  bn.  It  should  therefore 
be  noted  that  the  computations  thus  developed  are  not  self-consistent,  as 
Thomas  himself  has  pointed  out . 

Ignoring  fob  the  moment  questions  involved  in  the  simplified  compu- 
tations of  Thomas,  we  will  follow  his  assumptions  and  compute  Hp  from  the 
other  pairs  of  lines  to  check  the  consistency  of  the  results.  Since  the 
series  for  converges  more  rapidly  in  our  revised  formula  discussed 

in  Chapter  III,  we  can  now  perform  the  computal ion  for  pairs  of  lower  order 
lines. 


A3 


Table  X. 


The  population  of  the  Balmer  ground  state 
along  the  line  of  sight. 


n-n+lNj 

670  km 

1500km 

2330km  3170km 

25  - 30 

2.53 

20  - 25 

1.77 

15  - 20 

0.82 

0.46 

10  - 15 

0.68 

0.52 

0.22 

10  - 13 

0.87 

0.67 

0.33  0.13 

The  results  of  computations  are  shown  in  Table  X.  The  first  entry 
in  each  column  corresponds  to  the  value  Thomas  computed  at  each  height", 
is  va  have  discussed  in  Chapter  III,  p.23,  Thomas  used  half  the  popu- 
lation of  the  Balmer  ground  state  in  place  of  the  total,  so  that  the 
results  shown  in  Table  2 in  his  paper  should  be  multiplied  by  two. 

We  notice  in  Table  X that  the  values  are  considerably  scattered  for 
the  different  combinations  of  lines.  However,  the  fluctuations  are  not 
random.  The  values  decrease  rather  systematically  towards  the  lower 
order  line  combinations  for  the  lowest  height,  while  the  effect  is  in 
the  opposite  direction  for  the  next  lowest  height.  Considering  this 
change  in  the  direction  of  the  effect  for  two  different  heights,  as  well 
as  the  dependence  of  the  results  on  a particular  choice  of  combination 
of  two  lines,  one  might  s^ply  attribute  the  effect  to  the  random 
scattering  of  the  observed  data.  . Before  doing  so,  however,  we  should 
consider  how  the  simplifying  assumptions  made  in  the  computations  might 
affect  the  results.  Let  us  consider  the  direction  of  the  contribution 
of  in  equation  (2.8).  Since  b,%>  ! and  it  approaches 

unity  mcnotonically  as  n increases,  we  see, 

( ^ C 'by*'  ~ *V+  fcT)  dd  < !ry  , (2.9) 

That  is,  the  inclusion  or  b ^ adds  larger  numbers  for  smaller 

n than  for  larger  n to  the  left  hand  aide  of  equation  (2.8),  which  means 
the  relative  absorption  should  be  larger  for  the  lower  lines 

than  the  amount  we  have  used  in  the* computation.  Hence  the  contribution 
due  to  the  bn  factor  would  increase  the  values  of  N2  more  for  combinations 
of  lover  than  for  higher  order  lines.  The  direction  of  the  correction 
expected  by  the  inclusion  of  the  departure  from  thermodynamic  equilibrium 


thus  seems  promising  (at  the  lowest  height).  If  we  assume  the  chromo- 
spheric temperature  to  be  as  low  as  6000° , then  the  effect  of  the  inclu- 
sion of  the  Boltzmann  factor  A*.  e*7'-  will  also  be  to  increase 

No  in  the  same  direction  although  the  amount  of  this  effect  is  smaller 
than  the  inclusion  of^^i*^.  bu« 

As  far  as  th6  values  shown  in  Table  IX  are  concerned,  the  above  argu- 
ment holds  only  for  the  lowest  height.  For  the  height  of  1500  km  the 
first  two  entries  show  the  opposite  effect.  The  increase  in  the  last 
row  for  all  three  heights  is  not  easily  justified  without  knowing  the 
most  probable  values  of  bn  for  n higher  than  10,  because  the  interval  of 
n is  not  the  same  as  for  the  others.  We  also  cannot  depend  too  greatly 
on  the  disagreement  at  1500  km  compared  to  the  lowest,  since  the  calcu- 
lations are  made  for  only  two  combinations.  We  realize  that  the  relative 
intensity  of  to  H 2 r at  1500  km  implies  a negative  absorption,  which 
of  course  is  an  indication  of  the  observational  error. 

Thu3,  we  again  note  the  possibility  of  a considt3rable  dependence  of 
the  results  on  the  choice  of  two  lines.  We  may  not  be  able  to  draw  any 
conclusions  from  data  at  greater  heights  unless  we  make  computations  for 
a much  larger  number  of  line  comcinations . The  importance  of  the  sys- 
tematic variations  of  H2  values  for  670  km,  therefore,  is  not  affected 
by  the  failure  of  data  at  greater  heights  to  follow  that  system.  Thus 
we  conclude  that  it  is  worth-while  to  consider  methods  of  including  the 
bn  factor  in  the  calculation. 

§3 • Consideration  of  the  effect  of  self-absorption,  including  departure 

from  thermodynamic  equilibrium. 

There  may  be  two  alternative  methods  to  include  departures  from 
thermodynamic  equilibrium.  In  order  to  determine  the  amount  of  self- 
absorption we  have  essentially  two  parameters,  bn  and  N2,  both  of  which 
are  functions  of  Te  and  Ne.  From  the  slitless  observations  of  the  flash 
spectra,  only  the  continuum  measurements  have  provided  a direct  detei'- 
mination  of  Te  and  Ne  as  a function  of  chromospheric  height.  That  is, 
from  the  free-bound  emission  we  can  determine  Te  and  Ne  and  the  electron- 
scattered  photospheric  continuum  gives  us  Ne  independently.  If  we  thus 
have  the  values  Te  and  Ne  at  each  height,  we  can  compute  first  bn,  then 
Nn,  theoretically.  Hence  one  method  is  to  substitute  the  theoretical 
values  of  bn  thus  computed  into  equation  (2.8)  to  estimate  the  real  amount 
of  self-absorption  from  the  measured  line  intensities.  Then  from  this 
amount  of  self-absorption  we  can  determine  the  value  of  N2  at  each  height. 
Since  N2  is  a function  of  b2,  Te  and  Ne,  we  may  recompute  the  values  of 
Te  and  Ne  from  the  line  intensities.  The  second  method  is  to  assume  M2 
at  each  height,  and  from  it,  compute  the  amount  of  self-absorption.  We 
then  use  the  self-absorption  to  make  corrections  for  the  observed  line 
intensities,  and  from  the  corrected  line  intensities  and  equation  (2.8) 
we  determine  bn  at  each  height.  Then,  from  these  values  of  bn,  we  can 
estimate  Te  and  Ne.  Thus,  from  either  one  of  these  methods,  we  can 
obtain  a preliminary  distribution  of  Te  and  Ne  that  we  cun  use  to  start 
the  second  approximation  in  the  analysis.  Repeating  the  same  procedure 
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we  may  finally  determine  a self-consistent  model  of  the  chromosphere  from 
both  line  and  continuum  measurements. 

At  the  present  time,  unfortunately,  tje  have  only  Wildes  determination 
(1947)  of  the  K2-gra.dis.nt,  which  was  obtained  by  setting  the  emission  gra- 
dient equal  to  the  density  gradient;  but  ye  have  no  absolute  values  of  N2 
except  the  one  by  Thomas  discussed  in  the  preceding  section.  On  the  other 
hand,  the  application  of  continuum  measurements  mentioned  above  involves 
various  uncertainties.  Separation  of  the  free-bound  emission  and  the 
scattered  continuum  is  difficult.  Second,  since  Te  appears  only  in  the 
Boltzmann  exponential  factor  of  the  freo-bcund  emission,  ve  have  to  depend 
on  small  differences  in  the  observed  data  to  determine  Te,  so  that  small 
errors  in  the  measurements  are  of  importance. 

further,  the  theoretical  computation  of  bn  involves  a great  amount 
of  computation.  We  have  available  only  the  asymptotic  values  of  for 
n higher  than  10,  computed  for  the  case  Te  = 35,000°,  to  which  ye  have 
referred  in  section  2.  We  must  also  none  that  the  final  determination 
of  N2  in  the  first  method  requires  another  laborious  computation,  if  we 
follow  the  same  procedure  described  in  section  2 involving  every  pair  of 
lines  in  the  analysis,  because  of  the  difficulty  of  solving  the  exponen- 
tial formula  to  obtain  N2. 

It  is  therefore  desirable  to  determine  a T0  and  Ne  distribution  irom 
the  line  intensities  only,  and  then  to  compare  these  results  with  the 
results  obtained  from  the  continuum.  In  general,  the  temperature  deter- 
minations discussed  in  Chapter  I that  were  based  on  line  intensities 
yielded  higher  values  than  those  based  on  continuum  intensities.  It  may 
be  that  line  emission  and  continuum  radiation  originate  under  different 
atmospheric  conditions. 

The  following  method,  based  on  least  square  analysis,  provides  a 
somewhat  independent  determination  of  bn  and  N2  and  minimizes  as  well 
the  amount  of  computation  necessary  to  include  all  the  line  intensities 
at  a certain  height.  The  least  square  calculation  assumes  that  all  the 
small  fluctuations  from  the  theoretical  Gainer  decrement  may,  after  in- 
cluding both  effects  of  the  departure  from  thermodynamic  equilibrium  and 
the  self-absorption,  be  attributed  to  observational  errors.  If  we  assume 
some  expression  for  bn  as  an  analytical  function  of  n,  we  may  express  the 
Balmer  line  intensity  as  a function  of  n.  Ihus  we  write  equation  (2,5) 
in  the  following  logarithmic  form,  retaining  only  terms  depending  on  the 
quantum  number,  n: 


where 


(3.1) 


and 
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^e’KTe* 


We  node  that  the  Boltzmann  factor  is  also  included,  and  vs  have  ig- 
nored the  radiation  induced  emission.  In  the  above  expression  A and  B 
are  independent  of  n but  involve  Te  and  b2. 

Now  we  wish  to  express  bn  ae  a function  of  n.  Using  Thomas'  asymp- 
totic values  for  bn  mentioned  in  section  2,  we  find  the  following  to  be 
the  most  simple  formula; 

$ T 

+ ~ (3.2) 


Actually,  the  first  two  terms  are  sufficient  to  fit  Thomas'  values. 
Table  XI  shows  the  comparison  of  the  computed  values  and  the  values 
obtained  from  the  formula: 


(3.3) 


Table  XI. 

Theoretical  values  of  b^ 


n 

10 

13  15 

20 

25 

30 

Thomas 

0.68 

0.59  0.55 

0.48 

0.44 

0.42 

Equation  (3.3) 

0.68 

0.59  0.55 

0.49 

0.45 

0.42 

Thus,  we  aesume  that  the  expression  (3.2)  with  the  third  term  intro- 
duced is  adequate  for  the. whole  range  of  electron  temperatures.  Substi- 
tuting (3.2)  into  equation  (3.1),  we  have 

(3.1*) 

where 

c,=  A + p , Ci=  % > c,=  B +r. 

Since,  the  coefficients,  f and  T of  the  bn  formula  are  combined 
with  the  other  constants  Including  Te,  we  cannot  determine  the  absolute 
value  of  bn*  Since,  however,  as  we  will  show  later,  the  coefficient  C3 
is  negligibly  small,  and  we  can  obtain  the  ratio  bn  / fc>n>  . 


47 


The  principle  of  the  present  method  is  to  assume  a certain  value 
for  K2  from  which  we  compute  [ A bsnl  to  substitute  into  equation  (3.4) 
From  this,  we  compute  a correction  factor  and  the  coefficients  for 
bn  by  assuming  that  the  observed  Calmer  decrement  can  be  fixed  by  (3-4) 


with  a small  differential  correction 


A I" 


lhus  we  obtain  both  a 


corrected  value  Ng  and  the  ratio  ^ ^ , which  we  use  in  repeating 

the  analysis  successively.  After  adding  a differential  of  log 
we  can  arrange  equation  (3.4)  in  the  following  form  for  convenience  of 
computation: 


rT'  +■  c3  ti3*  + 


where 


B N. 


<p 


(3.5) 


oo 

2.  c-iy 

■m=  i 


4-  Y)*- 

2*1  C N a,  S< 


O-Vl 


} 


>- 1 


C ^41)1  C T*  t I ^3/jl  ^ ^ 0 


B N 


a. 


X,  c-b 
‘0 


’"V'  C hi  2-  -J-Vi 


-^,^0  + (Wv-t-l)  /2- 

j represents  the  observed  quantity,  and  we  can  determine  C^,  C2, 
C3,  and  A Ni<  One  advantage  of  this  method  is  that  we  do  not  have  to 
assign  a definite  value  for  Te  and  Ke.  Since  C|  is  determined  by 
least  square  calculation,  values  of  electron  temperature  and  electron 
density  above  a certain  height  may  be  fixed  observationally  through  the 
coefficient  C7. 


§4.  Numerical  results  and  discussion. 

The  calculation  was  carried  out  first  neglecting,  then  including, 
the  term  C /ni- In  equation  (3.5)  and  comparing  the  two  results.  It  was 
concluded  that  this  term  was  not  appreciable  at  all,  which  implies  that 
the  exponential  Boltzmann  factor  was  negligible,  and,  also,  that  the 
last  term  of  equation  (3.2)  is  not  significant  for  a bn  approximation. 

In  general,  we  have  included  every  line  for  which  the  convergence  of  the 
absorption  terra  is  sufficient  to  carry  out  the  computation  at  each  height, 
excluding  the  few  lines  recognized  in  the  literature  as  blended  with  other 
lines.  According  to  Menzel  (1931),  H3Q,  H23,  H25,  Hi  , and  Hg  are  blended 
with  strong  lines  of  Ti+,  Y+,  V* , Ti+,  and  He,  respectively.  H22  is  also 
blended  with  a weak  line  of  Fe,  but  is  included.  For  the  lowest  height, 
however,  we  first  included  all  these  blended  lines  and  then,  for  comparison 
purposes,  we  repeated  the  same  analysis  excluding  them.  We  found  that  the 
successive  computations  did  xiot  show  a monotonic  convergence,  probably 
indicating  that  the  observational  data  were  too  inaccurate  to  determine 
the  numerical  figures  to  the  high  precision  we  expected.  Therefore,  in 
the  following,  we  will  show  the  number  obtained  at  each  step  of  the 
successive  computations. 

i).  Result  for  the  lowest  height,  h = 670  km. 

The  lines  higher  than  H15  were  included.  We  took  0,8  6 
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and  obtained  A 6,13>  and  Cx  * 3.7  . To  see  the  effect  of 

small  computational  fluctuations,  we  interchanged  the  order  of  the 
C 2 and  the  Terms  in  equation  (3.5).  This  gave  C'a.^'S.G 

and  A C Hj.^t30>6,|a..  It  is  interesting  to  note  that  A C « o, 

was  also  obtained  from  the  analysis  covering  only  the  lines  from  H25 
to  H^o* 

All  the  blended  lines  were  included  in  the  above  computation. 
Repeating  the  analysis  with  the  same  initial  value  Na.^^  —0.86  t 
but  excluding  H22»  H25,  ^28  an^  ^30»  ve  obtained  A N»y-j6  = 0,  |o  and 
£ =5  3 , £,  . Then,  as  the  second  approximation,  we  started 

from  the  value  30  =*0,95"  • We  obtained  and 

Cz  = 3,3  from  both  equation  (3.5)  and  the  same  equation 

with  the  order  of  two  terms  interchanged.  Thus,  both  approximations 
give  Na&ss-  j and  indicate  that  this  value  is  reliable. 

However,  the  accuracy  of  the  present  data  is  not  sufficient  to  give 
the  value  to  two  decimal  places  for  a reason  that 

will  be  pointed  out  later. 

The  final  normal  equations  derived  from  the  analysis,  from 
which  the  above  values  were  computed,  are  as  follows: 

for  NA°<i0=  0.86,  ^ + Cay  + - <f 

13  6.3438  -4.3590  0.5440 

3.2168  -2.2449  0.2970 

1.6007  -0.2100 


for 


°-95> 


-0.1210  0.1177  0.0315 

-0.1391  0.0282 

0.0246  0.0024 

= 0.10, 

-4.0868  0.7350 

-2.0902  0.3947 

1.3805  '-0.2588 

0.0958  -0.0360 

-0.957  0.0277 

0.0199  0.0008 


8ince  the  observed  data  is  given  to  only  two  decimal  places, 
only  the  first  two  figures  after  the  decimal  point  are  significant 
in  the  above  calculations.  Therefore,  for  Uxo(y 0 0.9  S'  the 

value  A ( =>  has  no  meaning  and  even  for  * 0,26  t 

is  not  actually  determinate  within  the  indicated  accuracy. 
This  means  that  we  cannot  determine  more  than  one  decimal  place  for 
M H • Thus  from  this  analysis,  ve  can  only  state  that  at 
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Ritnofc*-  (4.1) 

3-'0‘  , 

It  may  be  interesting  to  note  that  the  value  C2  - 4 agrees  well 
with  Thomas ' numerical  solution  of  equat3.on  (3.3). 

The  same  computations  for  this  height  were,  repeated,  including 
the  term,  Cx/r^  in  (3.5).  It  was  found  that  C3  was  zero  for  all 
cases,  so  we  concluded  that  this  term  was  not  appreciable  for  the 
present  computation.  We  will  discuss  in  a later  section  tne  accuracy 
of  the  results  shown  in  equation  (4*1). 


ii).  Results  for  h = 1500  km. 


We  will  next  show  the  normal  equations  computed  for  the  height 
h = 1500  km  for  various  initial  values  of  N 2.  • Again,  alD 

lines  from  Hi 5 +0  H30  except  the  blended  linds  were  included  in  the 
computations.  Actual  computations  were  carried  to  four  decimal 
places,  but  only  three  decimals  are  shown. 

c.  + c.-n-1  “ 'f 

K^a=  0.18  13  6.3438x10 

X 3.2168 

-7.140 

-3.844 

5.005 

-1.249 

-0.597 

0.654 

-0.1210 

0.359 

-1.084 

-0.013 

0.032 

0.017 

0.006 

«=  0.  tf 

K*>„=  0.49 

-4.085 

-2.089 

1.380 

-0.391 

-0.139 

0.075 

0.096 

-0.096 

-0.052 

0.048 

0.020 

-0.007 

* -0,^ 

°*26 

-6.673 
-3.573 
4. 261 

-0.996 

-0.461 

0.440 

0.317 

-0.836 

-0.025 

0.062 

0.007  0.00/ 


^ K)  c4  •«=  0.£ 
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W o<,  = 0.86  -4.359  0.394 

-2.245  0.274 

1.601  -0.218 


0.118  -0.082 

-0.139  0.086 

0.025  -0.007 

^x**  ^ ~0'3 

We  see  that  for  all  four  initial  values  ranging  from 
to  0.86,  the  final  values  of  (f>  are  less  than  the  last  significant 
value,  0.01.  That  Is,  we  cannot  obtain  physically  significant 
values  for  even  to  one  decimal  place.  If  we  calculate 

using  the  third  figure  for  , v.’e  have  results  that  are 
inconsistent  with  what  we  would  expect  from  the  method  of  successive 
approximations,  for  N1^36«dii(6 > we  n higher  positive  value  of 
than  for  N3«<i()-=  ogsr  For  as  high  as  0.86,  4 CM**  a*) 

is  less  negative  than  for  0.49.  Thus,  we  cannot  obtain  even  the 
first  decimal  place  for  at  this  height.  Ihis  means  that  the 

effect  of  self-absorption  and  the  departure  from  thermodynamic 
equilibrium  in  the  higher  members  of  the  Balmer  series  are  not 
sufficient  to  lead  to  self-consistent  values  of  No  within  the  accuracy 
of  the  observation. 

The  same  computations  were  carried  out  for  the  height,  h = 2330km 
and  3170  km,  covering  the  lines,  Hqo  - H17  and  Hg  - Ifjj,  respectively. 
As  we  would  expect,  the  results  showed  that  the  situation  becomes 
even  worse  for  higher  regions,  and  the  conclusion  stated  above  is 
confirmed. 

Now  returning  to  the  value  shown  in  (4.1 ),  and  comparing  the 
value  at  this  height  with  those  at  greater  heights,  it  may  not  be 
possible  to  conclude  immediately  that  at  least  the  first  figure  of 
NiVio  is  determinable  even  at  the  lowest  height.  Unless  we  repeat 
the  calculations  taking  various  initial  values  at  large  intervals, 
such  as  N Wso  = 0.5,  1.0,  1.5,  we  cannot,  prove  that  the  result 
(4.1)  has  meaning.  However,  it  seems  unwise  to  repeat  such  laborious 
computations;  as  a matter  of  fact  it  is  impossible  for  Nq3>4,4  greater 
than  1.0,  because  of  the  slow  convergence  of  the  {Ats*]  aeries, 
"therefore,  in  the  next  section,  we  consider  some  simple  calculations 
designed  to  estimate  a rough  order  of  N2,  following  the  first  method 
discussed  in  section  3. 


§5.  Determination  of  No  and  No-gradlent. 

In  the  preceding  section,  we  found  that  successive  calculations  by 
ths  leaet  square  method  did  not  converge  to  determine  a solution  of 
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sufficient  accuracy.  The  fundamental  equation  (3.5)  was  chosen  for  the 
analysis  under  the  assumption  that  the  bn  factor  was  solely  responsible 
for  the  systematic  errors  of  the  earlier  method.  1'hat  is,  we  assumed 
the  non-id6ntical  values  of  N2  shown  in  Table  X were  due  to  the  omission 
of  bn.  It  is  not  very  clear  whether  this  basic  interpretation  was  faulty 
or  simply  that  the  accuracy  of  the  observation  was  not  sufficient  for 
obtaining  quantitatively  consistent  results.  In  the  discussion  of  Table 
X,  we  mentioned  that  the  systematic  decrease  of  N2-value  may  not  have 
appeared  for  a series  of  different  combinations  of  the  lines  and  suggested 
that  the  non-identical  r esults  may  be  attributed  to  observational  errors. 
The  only  way  to  prove  this  possibility  seems  to  be  to  repeat  the  same 
calculation  for  a sufficiently  large  number  of  different  line-combinations 
and  draw  a statistical  interpretation.  Since  we  have  found  some  approxi- 
mate formulae  for  the  self-absorption  series,  the  computation  described 
in  section  2 is  not  too  difficult.  Thus,  we  made  the  following  calcula- 
tions using  equation  (A. 3)  in  Chapter  III. 

1/e  sot  ^ 

j" ~ - 0.016  C 


Her.ce,  the  logarithmic  ratio  of  the  self-absorption  of  two  different  lines 
is  given  by 


(5.D 


Substituting  the  above  into  equation  (2.8),  we  have 


-A.4j  k-0.076  N;(a„<4{.(5.2) 


The  value,  N2,  can  be  determined  directly  from  the  a.jove  equation. 
The  first  and  second  term  of  the  left  hand  side  is  the  logarithm  of  the 
intensities  observed  to  the  theoretical  intensities  in  thermodynamic  ^ 
equilibrium  neglecting  the  self-absorption,  lhe  Boltzmann  factor,  <2.  K , 
is  included  in  the  second  term,  and  is  not  negligible  for  lower  tempera- 
tures. Apart  from  the  bn  term,  the  temperature  factor  appears  only  in 
the  absorption  coefficient,  and  in  the  smal]  1 quantity  e*r'‘ . In 

the  following  computation,  we  consider  an  isothermal  atmosphere  at  two 
extreme  temperatures,  Te  = 6000°  and  Te  = 35,000°.  For  Te  - 6000°  we  can 
ignore  the  term  b^  since  the  assumption  of  the  thermodynamic  equilibrium 
may  hold.  If  we  set  Te  - 35,000°,  the  greatest  amount  of  uncertainty 
comes  in  through  the  term  b-r.  . Vfe  have  made  computations  for 

two  cases,  one  neglecting  b which  corresponds  to  the  compu- 

tation shown  in  Table  X,  and  the  other  including  the  asymptotic  values 
of  discussed  in  section  3 end  computed  from  equation  (3.3). 

Prom  recent  results,  as  we  have  discussed  in  Chapter  I,  the  most  probable 


Table  XII 
No  x lCf16 


3.14  *35.7 
2.72  * 9.15 


1.65 

4.U 

1.23 

2.49 

1.03 

3.91 

0.72 

2.26 

0.85 

1.13 

*0.26 

0.93 

1.12 

1.02 


0.08 


5.22 
U .li3 

l.uu 

4.33 

tt 

ii.ue 

*13.86 

5.U1 

0.76 

1.23 
0.63 

*0.53 

*0.27 

0.65 

1.07 

1.28 


0.95 

tt 

# 

tt 

tt 

0.36 

0.39 

2.30 

0.6? 

0.79 

0.67 

0,75  *1.00 
0.75  *0.58 
0.08  0,37 

tt 

* 4.60 
*11.03 
* 6.11 
§ 

U 

tt 

2.14 

2.56 

0.97 

# 

1.97 

1.38 

0.71  0.53 
# 0.81 
0.30  0,22 
0.94  1.07 


Te  - 
670 

35000*  "°  bn  _ 
inc.  e"1 

1500  2330  3170 

2.53  *33.1 

2.5l  * 8.40 

1.49 

0.86 

3.74 

tt 

1.03 

tt 

1.79 

tt 

0.94 

tt 

3.41 

0.31 

0.63 

0.33 
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0.72 
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0.88 

0.63 

*0.22 

0.51 

0.69 

0.54  *0.75 

0.76 
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I 0.71 

0.52  0.25 

# * 4.47 
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1.36 
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// 

4.23 

*12.60 

4.92 

0.69 

1.13 

0.56 
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0.54 
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0.07 
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tt 

*0.45 

*1.18 
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0.34  0.13 

0.79  0.35 

0.67  0.29 

0.14  0.11 

0.37  0.30 

0.30  0.24 

0.852  0.422  0.221  0.107 

2.223  1.086  0.537  0.280 

1.966  0.9120.4370.223 

0.85  0.78  0.86 

0.86  0.65  0.78 

0.92  0.91  0.80 

53 


value  of  Te  seems  likely  to  be  lower  Limn  35,000°,  so  that  the  values, 

A*  b-n.  we  used  may  be  considered  a3  the  maximum  values,  further, 
since "wo  will  see  later  that  the  final  results  are  not  too  sensitive  to 
the  bn  factor,  we  can  determine  M'2. 

The  values  of  NpxlC"^  thus  computed  are  tabulated  in  Table  XII,  For 
each  height,  we  included  every  line  higher  than  H^o  in  our  computations. 

For  iines  lower  than  Hpo,  the  self-absorption  becomes  too  large  for  appli- 
cation of  our  approximate  formula.  The  first  column  of  the  table  gives 
the  numbers  designating  the  pair  of  lines  taken  for  the  computation  of  the 
N2  values  in  the  corresponding  rows.  While  there  is  no  reason  for  selec- 
tion of  particular  pairs  of  lines,  we  computed  with  every  possible  combina- 
tion of  line  numbers  having  intervals  of  five  and  three,  for  the  two 
highest  regions,  we  introduced  also  the  combinations  oi  lines  having 
intervals  of  two  in  order  to  have  more-  data.  Thus  we  have  repeated  the 
computation  a sufficient  number  of  times  to  obtain  for  the  three  lower 
heights  enough  results  to  carry  out  a statistical  analysis.  For  the 
highest  region,  the  number  of  lines  observed  is  not  large  enough  to  get 
a reliable  result. 

In  Table  XII,  the  blank  spaces  occur  where  observational  data  are 
not  available.  The  entires  marked  by  # indicate  that  the  relative  inten- 
sities corresponding  to  these  lines  show  the  wrong  sign  for  self-absorption. 
That  is,  the  basic  interpretation  for  carrying  through  these  computations 
leads  to  inconsistencies.  As  we  see  in  the  Table,  # appears  primarily  for 
the  height,  1500  km,  and  in  the  table  as  a whole,  the  great  majority  of 
the  results  are  consistent.  We  see  that  the  I^-values  are  scattered  at 
random  over  n rather  than  systematically  decreasing  towards  the  lower  n 
as  in  lable  X.  'that  is,  the  non-identity  of  N2  seems  more  likely  to  be 
caused  by  the  observational  error.  If  we  compare  the  results  for  two 
cases  of  Te  = 35,000°,  i.e.,  the  cases  with  and  without  the  bn  factor,  we 
see  the  differences  are  very  small.  This  implies  that  the  bn  term  is  not 
appreciable  and  not  large  enough  to  correct  the  decrease  of  N2  in  Table 
X.  It  is  also  interesting  that  values  of  11?  in  Table  XII  show  consistent 
decrease  with  height  for  each  combination  of  lines.  Thus  we  conclude  that 
variation  of  M2  at  a given  height  is  not  systematic  but  random,  and  one 
may  take  an  arithmetic  mean  to  determine  a most  probable  value  of  N2  at 
each  height.  In  averaging  all  the  values  shown  in  Table  XI,  we  exclude 
the  entries  designated  by  #,  because  of  the  unreasonably  large  variation 
of  these  values  from  the  mean. 

The  mean  values  thus  computed  and  the  logarithmic  height  gradient  of 
N2  are  shown  in  the  last  two  row3  of  Table  XI.  It  may  be  interesting  to 
note  that  the  height  gradient  of  N2  is  very  insensitive  to  a change  in 
temperature  from  one  to  the  other  extremes,  in  spite  of  a difference  of 
more  than  a factor  of  2 in  the  absolute  value  of  N2.  This  may  be  a 
necessary  consequence  of  the  assumption  of  a constant  temperature.  We 
note  that  in  equation  (5.2)  all  the  temperature  dependent  terms 
L , cCiVl  0 were  taken  as  constant  throughout  the  atmosphere.  The 
No  height  gradient  for  Te  = 6000°,  where  bn  is  equal  to  unity  and  tho  tem- 
perature is  constant,  is  subject  tc  less  uncertainties  than  for  the  higher 


temperature.  It  should  ho  equal  to  half  of  the  alec  Iron  density  gradient 
and  also  to  the  true  emission  height  gradient  (i.e.  the  observed  emission 
gradient  corrected  fax-  sel  f-  ahsorplion)  • Iho  r- salts  agree  very  closely 
with  Wildt' s (1947)  determination  of  the  emission  gradient.,  0lll'iQSO'-1 
and  a r<.  smaller  than  the  value  of  Gillie  and  I-'enzel.  1 . \ o s ^ which 

is  not.  corrected  for  self-r. bsor y. li on . Thomas  (1950)  and  Athay  et  a~(.  (3,954) 
have  argued  that  the  density  gradient  roust  be  higher  than  the  Wildt  emission 
gradient  because  of  the  effect  of  the  bn's,  which  Wildt  ignored.  Regard- 
less of  this  point,  however,  the  present  results  show  practically  no 
difference  after  including  bn.  further,  the.  small  difference  between  the 
last  two  cases  in  7 able  XII  indicate  a direction  of  increase  opposite  to 
their  suggestion.  It  is  also  important  to  note  that  the  present  Ng- gra- 
dients are  almost  constant  throughout  the  different  heights,  whereas 
Thomas'  x-esults  showed  a large  decrease  towards  ihe  upper  r egion  of  the 
chromosphere.  Despite  the  fact  that  the  values  obtained  for  the  highest 
height,  h = 3170  km,  in  Table  XII  are  most  uncertain,  because  only  three 
computations  were  averaged,  the  final  results  of  N?- gradient  seem  consis- 
tent with  the  other  heights. 


Since  the  ^-gradient  seems  constant  through  all  heights,  we  computed 
the  most  probable  N2-gradient  by  the  least  square  analysis.  The  following 
results  were  obtained: 

For  Te  = 6000°,  p.  « 0.^3  • \o  ^ . 

ft  ~ t#', 


For  Te  = 35,000°,  including  bn, 

For  Te  = 35,000°,  neglecting  bn,  ^ « o.  87  • \S* 


(5.3) 


Or,  we  can  write  the  following  expression  for  log  N2  as  a function 
of  the  chromospheric  height,  h: 

for  Te  = 6000°,  ^ U>.  I n - 0. 35-9  **6®  TL  ^ . (5.4) 

For  Te  = 35,000°  including  b , JU  |\)  = 16,5"?  - 0,34  3L  MO* . 

* (5.5) 

For  Te  = 35,000°  neglecting  bn,  ^ M _ (£,5^-0,  i')')  X |0* 

> • (5.6) 


§6.  Discussion. 

In  the  previous  section,  we  derived  expressions  for  N2  as  & function 
of  height  for  cases  of  Te  = 6000°  and  35,000°.  The  effects  of  departures 
from  thermodynamic  equilibrium  were  examined  by  including  a bn  factor  in 
the  case  of  Tfe  = 35,000°.  The  important  r esult  is  that  the  density-height 
gradients,  p , are  practically  the  same  tor  both  temperatures.  The  differ- 
ence between  the  values  of  N2  for  the  two  temperature  cases,  at  each  height, 
is  due  to  the  temperature  dependence  of  cfx-n,  o • '1‘hat  is,  since  we 
evaluated  the  absorption  coefficient  at  the  line  center  using  the  approxi- 
mation  that  the  line  profiles  are  determined  by  Doppler  broaden ings , 
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0 varies  as  (Te)~a.  Hence,  it  is  necessary  to  know  the  value  of  Te 
to  determine  N2.  However,  we  should  re-emphasize  that  the  empirically 
determined  quantity  Je  the  product  td  ^ o<  a.y,  g , and  that  the  only  assump- 
tion concerning  temperature  necessary  for  its  de  termination  is  that  Te 
is  constant  throughout  the  chromosphere,  the  absolute  value  assigned  to 
Te  does  net  affect  this  determination.  Except  for  calculating  the  values 
of  £**  and  , we  introduced  values  of  Te  - 6000°  and  35,000°  only  in 
the  final  computation  to  deduce  N2  from  0 • If  we  carried  out 

the  analysis  ignoring  e^and  b^.hen  Q would  necessarily  be  t he  same  for 
Doth  cases.  The  insignificant  difference  we  found  for  p>  in  two  cases, 
implies  that  the  numerical  results  are  not  sensitive  to  departures  from 
th enno dynamic  equilibrium.  This  conclusion  is  confirmed  further  by  the 
following  calculations. 

Ignoring  the  factors  bn  and  e11,  we  consider  the  variation  with 
temperature  of  the  constant  term  in  equation  (5.4)  and  (5.5)  to  be  due 
solely  to  the  variation  of  O^xn^O  • Then  we  express  the  constant  t erm 
as  a_(  function  of  Te.  Since  N2  varies  as  (Te)a  for  varying  as 

(Te)'c,  ye  can  write  the  following  expression  from  equation  (5.4), 

^ Nx“'6-i3 

From  the  above  equation,  we  get  the  value  of  16.55  as  the  constant  term 
for  Te  = 35,000°.  This  differs  from  that  in  equation  (5-5)  by  only  0.03. 
Thus  we  may  conclude  that  the  factor  of  bn  does  not  appreciably  affect 
the  present  N2  results.  We  may  therefore  consider  that  equation  (6.1) 
glve3  the  most  probable  N2  distribution,  at  least  for  a chromospheric 
region  in  which  Te  is  lower  than  35.000°.  It  should  be  noted  that  the 
N2-height  gradient  in  equation  (6.l)  is  determined  empirically  without 
any  a priori  assumption  on  the  value  of  Te  or  . Hence,  if  the  electron 
temperature  at  each  height  is  given,  v/e  can  determine  N2  for  the  corres- 
ponding height  directly  from  equation  (6.1 ). 

As  was  discussed  in  the  above,  it  appears  that  the  1932  data  show 
very  small  departures  from  thermodynamic  equilibrium.  This  implies  that 
the  chromospheric  electron  temperature  was  much  lower  than  35,000°.  How- 
ever, the  small  value  v/e  obtained  for  , compared  to  the  emission-height 
gradient  determined  by  Cillie  and  Menzel  (1935)  for  the  higher  order  lines, 
gives  sojne  evidence  that  the  effects  of  b2  may  be  present  and  that  the 
temperature  increases  with  height  in  the  lower  chromosphere.  In  this  case, 
the  value  of  may  vary  for  different  lines.  We  will  return  to  xhese 
points  in  later  sections  in  the  following  chapter. 
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V.  APPLICATION  TO  THE  FLASH  3PECTRTTK  OF  1952  ECLIPSE. 


§1 . Introduction. 

In  the  preceding  chapter,  we  have  considered  several  methods  for 
determining  the  absolute  value  of  the  population  of  hydrogen  atoms  in 
the  second  state,  N2,  from  the  numerical  values  of  self-absorption  in 
the  higher  order  lines  of  the  Calmer  series  and  have  applied  these 
methods  to  the  spectra  taken  at  the  1932  eclipse.  We  discussed  critically 
the  uncertainties  appearing  in  the  previous  analysis  by  Thomas  and  com- 
puted the  most  probable  values  of  No  at  each  height  within  the  accuracy 
of  the  observations. 

In  his  earlier  work,  Thomas  (1950a)  suggested  that  the  above  re- 
sults may  be  extended  to  a determination  of  the  chromospheric  tempera- 
ture distribution  as  well  as  an  electron  density  distribution,  thereby 
deducing  a model  of  the  chromosphere  from  flash  spectra.  The  1952  eclipse 
program  of  the  High  Altitude  Observatory'  under  the  supervision  of  Roberts, 
Evans  and  Thomas,  was  primarily  to  obtain  data  for  this  purpose.  A highly 
improved  jumping  film  technique  made  it  possible  to  obtain  the  flash 
spectrum  with  much  greater  resolution  of  chromospheric  heights.  The 
photometric  measurements  and  a preliminary  discussion  of  the  results  has 
been  given  by  Athay,  Dillings,  Evans,  and  Roberts  (1954).  The  intensities 
of  the  higher  Balmer  lines  and  the  Balmer  continuum  were  obtained  at 
approximately  every  100  km,  starting  from  the  height  of  500  km  above  the 
base  of  the  chromosphere.  A method  has  been  developed  to  determine  both 
the  electron  temperature  and  the  electron  density  simultaneously  by 
combining  the  free-bound  emission  and  the  electron  scattered  photospheric 
continuum.  However,  the  method  itself  Involves  various  uncertainties  as 
discussed  in  the  above  mentioned  paper.  In  fact,  the  results  show  a 
considerable  range  of  values,  depending  upon  the  different  assumptions 
made  in  the  calculation.  It  is  therefore  desirable  to  consider  the 
possibility  of  analyzing  the  line  intensity  data  independently,  especially 
since  the  eclipse  program  was  planned  primarily  for  the  line  analysis 
following  the  basic  principles  described  In  the  preceding  chapter. 

The  photometric  measurements  of  the  1953  spectra  were  carried  out  for 
2 points  on  the  arc,  and  the  heights  for  each  point  were  determined  inde- 
pendently from  moon  profile  pictures.  It  turned  out  that  the  heights  for 
the  2 points  differed  by  the  same  amount  that  occurred  between  six  con- 
secutive exposures,  and,  as  a result,  at  many  of  the  heights  intensities 
were  measured  at  both  points  on  the  limb.  In  the  following  computation 
we  took  the  mean  values  of  the  intensities  for  these  heights  in  all  cases. 
In  section  2,  we  shall  reconsider  the  method  of  determining  bp  and  N2 
simultaneous 1.  and  in  the  later  sections  the  general  characteristics  of 

the  Balmer  decrement  will  be  examined. 


§2 , Method  of  simultaneous  deterr.ii nation  of  Np  and  bn . 

In  sections  3 and  4 of  Chapter  IV,  we  developed  a method  to  uetermine 
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bv,  and  H2  simultaneously  from  the  numerical  amount  of  self-absorption. 

It  was  a method  of  successive  trial  and  error,  based  on  the  least  square 
principle.  We  found  that,  using;  the  1932  data,  repeated  calculations 
did  not  give  consistent  solutions.  As  ve  discussed  in  detail,  it  was 
not  certain  whether  the  cause  of  the  failure  was  due  to  insufficient 
accuracy  of  the  observations  for  this  detailed  calculation,  or  if  the 
theory  itself  was  not  adequate  to  explain  the  general  character  of  the 
phenomena.  The  method  was  first  developed  from  the  consideration  that 
from  the  1932  data  the  relative  absorption  at  the  lowest  height  decreased 
with  increasing  quantum  number.  This  decrease  was  in  such  a direction  as 
would  be  explained  by  a monotonic  decrease  of  bn  with  quantum  number.  We 
have  mentioned  that  the  systematic  decrease  at  the  lowest  height  was  con- 
cluded from  only  a part  of  Ihe  data,  however,  it  was  found  that  when  all 
lines  \jere  considered  the  systematic  change  with  quantum  number  was  not 
so  evident.  The  next  possible  fault  in  the  theory  may  lie  in  the  assumed 
form  of  the  bn,  because  the  theoretical  form  for  a kinetic  temperature  of 
35,000°  was  used.  We  shall  therefore  examine  the  general  character  of 
the  1952  data,  with  regard  to  these  two  points  before  going  into  the 
calculations.  Since  the  number  of  heights  involved  in  the  1952  data  was 
much  greater  than  in  1932,  we  may  draw  more  affirmative  conclusions  by 
statistical  methods. 

We  define  the  quantity  byj 

and,  in  Figure  III,  we  plot  Rn  at  each  height  against  the  quantum  number 
n.  denotes  the  observed  intensity  of  line  n and  En^dT. E.)  the 

theoretical  value  computed  for  thermodynamic  equilibrium  without  self- 
absorption. Hence  R^  represents  the  combined  effects  of  bn  and  self- 
absorption. Since  E"Wi(T.t.)  is  a function  of  M2,  and  N2  depends  on  Te, 

Ne  and  b2,  we  cannot  determine  the  absolute  value  of  Rn  without  knowing 
these  parameters  at  each  height.  However,  we  are  interested  only  in  the 
relative  amount  of  Rn,  which  can  be  determined  as  a function  of  n at  each 
height.  (As  long  as  we  assume  a spherically  symmetric  distribution,  Te, 

Ne  and  bq  are  constant  at  the  same  height,  00  that  ve  can  determine  the 
relative  values  of  JzmC T, E,)  by  computing  n-dependent  terms  only.)  The 
shapes  of  the  curves  in  Figure  III  demonstrate  the  general  characteristic 
of  bn  and  . If  the  chromosphere  is  in  thermodynamic  equilibrium 

and  self-absorption  effects  are  negligible,  the  curves  should  be  hori- 
zontal straight  lines.  'Ihe  effect  of  self-absorption  tends  to  lover  the 
line  below  the  horizontal,  and  the  departure  from  thenno dynamic  equili- 
brium affects  it  in  the  opposite  direction.  The  broken  lines  were  drawn 
horizontally  near  the  maximum  point  at  each  height  to  indicate  the  com- 
bination of  these  effects.  Ihe  solid  horizontal  lines  were  drawn  through 
the  average  value  of  Rn  at  each  height.  The  vertical  scale,  indicated 
in  the  top  curves,  is  the  same  for  all  neights. 

For  lower  heights,  the  curves  lie  lover  at  lower  quantum  numbers,  but 
this  effect  appears  to  decrease  systematically  with  height.  If  we  ignore 
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bj-j,  we  see  that  this  agrees  well  with  the  absorption  effect,  but  for 
greater  heights  the  effect  reverses,  i.e.,  ^ decreases  for  higher  lines 
and  the  amount  of  decrease  becomes  greater  for  the  greater  heights.  The 
latti-r  effect  may  only  be  explained  by  a changing  bn.  Since  decreases 
with  n,  we  may  get  a consistent  picture  for  this  observed  feature,  if  we 
assume  simply  that  bn  increases  with  height,  From  the  theoretical  con- 
sideration discussed  in  Chapter  II,  this 'implies  that  the  electron  tem- 
perature increases  upwards  in  the  chromosphere. 

The  last  two  curves  in  figure  III  were  drawn  for  the  1932  data  for 
comparison.  It  is  interesting  to  note  that  the  1932  data  do  not  show 
the  same  tendency  for  greater  heights,  an  interpretation  being  that  the 
bn  effect  was  smaller  in  1932.  Thus,  our  conclusion  fi'om  the  1932  data 
discussed  in  Chapter  IV  is  predicted  by  this  diagram.  For  the  1952  data, 
however,  in  contrast  to  these  1932  results,  we  see  that  the  bn  term 
should  not  be  excluded  if  we  are  to  get  a consistent  picture  over  the 
whole  region. 


Thus  we  shall  repeat  ihe  calculations  based  on  the  least  square 
principle  described  in  section  3,  Chapter  IV.  In  section  3,  we  used  an 
approximate  expression  for  log  brj  deduced  from  the  theoretical  values 
computed  for  Te  = 35,000°.  The  equation  is  hyperbolic  in  form,  with 
small  curvature.  Since  the  continuum  data  seem  to  indicate  a temperature 
lower  than  35,000°,  we  assume 

log  bn  = a + bn  (2.2) 


for  the  present  case.  We  also  see  in  Figure  III  that  the  above  form 
seems  sufficient  within  the  scattering  of  points  plotted.  Substituting 
equation  (2.2)  into  (3*1)  of  Chapter  IV,  we  have  the  following  formula 
corresponding  to  equation  (3.5). 


f\  t-  8 n 


-r 


°^,3D,0 


^ C^aj^sc 


(2.3) 


T = Em.-  3 S ^ - 4)  ‘ " S f dfcs»V 


For  convenience  of  computation,  we  have  written  the  differentiation  by 
instead  of  Ni.  . 

Hie  successive  computations  were  repeated  until  the  probable  erx'orB 
of  ACUi_«<i0')  and  B became  the  same  order  of  magnitude  as  their  solu- 
tions. Contrary  to  the*  former  results,  we  found  that  v/e  could  determine 
two  decimal  places  within  the  observed  accuracy  at  all  heights.  To  illus- 
trate this  point  both  coefficients  of  the  final  normal  equations  from 
which  we  obtained  the  solutions,  ^ 6^2*36^)  anc^  ® are  shown  below.  The 
probable  errors  are  also  added. 
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In  the  following,  the  lowest  line  Included  in  the  computation  is 
indicated  in  the  bracket  following  the  height.  Ihe  first  value, 
is  the  initial  value  of  Nr«2j>6/0  taken  for  the  computation,  and  the  last 
is  the  final  value  obtained.  lor  the  lowest  height, 
h = ‘>30  km,  the  last  two  computations  are  shown  in  order  to  illustrate 
the  convergence  of  the  computation,  and  for  all  the  other  cases  we  tabu- 
lated only  the  final  results,  -‘he  individual  computations  were  carried 
out  at  e very  other  height  for  the  lower  region,  and  for  h higher  than 
1820  km  we  found  solutions  at  each  height  in  order  to  minimize  the 
effects  of  random  errors.  Ihe  results  are  as  follows: 
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Since  the  observed  data,  as  we  have  discussed  in  Chapter  IV  , give 
only  two  decimal  places,  the  last  figure  in  the  above  results  for 
has  not  much  significance.  However,  from  the-  probable  errors,  we  may 
say  that  at  least  two  figures  are  reliable.  The  values  of  B represent 
the  logarighmic  ratios  of  two  successive  b^'s,  log  ^ ^ ( , and  are 

also  reliable  to  two  figures. 

The  directly  determined  solution  is  ^'a<^2>o/o  > instead  of  N2 
itself.  Again  we  recall  that  thQ  absorption  coefficient  introduced  here 
refex-s  to  the  absorption  at  the  line  Renter  so  that  it  13  a function  of 
temperature.  Since  U varies  as  (Te)Tf,  N2  varies  by  a factor  of  2.6  as 
Te  varies  from  5000°  to  35,000°.  The  values  of  N2  x 10~-'-°  are  shown  in 
Table  XIII  for  various  temperatures.  If  the  chromosphere  is  assumed  to 
be  isothermal,  then  the  logarithmic  h'2-gradient,  ft  , is  the  same  for  all 


Table  XIII. 
N2  x 10-16 


*r 

h km\sQ 

5000 

10,000 

20,000 

35,000 

530 

1.28 

1.81 

2.56 

3.39 

750 

0.92 

1.30 

1.84 

2.43 

960 

0.61 

0.86 

1.22 

1.61 

1180 

0.46 

0.65 

0.91 

1.21 

1390 

0.50 

0.71 

1.00 

1.32 

1610 

0.24 

0.34 

0.4S 

0.64 

1820 

0.21 

0.30 

0.42 

0.55 

1930 

0.19 

0.27 

0.37 

0.49 

2040 

0.13 

0.18 

0.26 

0.35 

2150 

0.064 

0.090 

0.13 

0.17 

2260 

0.060 

0.085 

0.12 

0.16 

2400 

0.102 

0.144 

0.20 

0.27 

3060 

0.042 

0.059 

0.83 

0.11 

3180 

0.015 

0.021 

0.030 

0.040 

temperatures.  In  Figure  IV,  log  N2  is  plotted  against  h2  for  Te  = 
35,000°.  Four  points  marked  by  / corresponding  to  h = 1390,  2150,  2260, 


FIGURE  ISC  Log  N2  vs.  Height 
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3180  km,  are  scattered  considerably,  but  all  the  rest  of  the  points  seem 
to  define  a straight  line.  Hence,  v:e  fitted  these  points  to  a straight 
line  by  ‘the  least  square  method.  The  point  for  h = 3060  km  is  apparently 
in  a good  position,  but  it  is  700  km  from  the  next  point,  and  we  have  to 
consider  the  relative  error  to  be  greater  for  the  weaker  lines.  Thus  we 
computed  two  cases  as  follows:  1)  including  every  point  except  four 

marked  by  X , and  2)  excluding  the  highest  point  from  l).  We  obtained 
the  following  formulae: 


1). 

log  Hp  = 16.79  (±0.023)  - 0.578  (*0.018) 

• 10"8  hcm  . 

(2.4) 

p = 1.33  • 10-8  cm"1 

2). 

log  N2  = 16.81  (±0.025)  - 0.597  (±0.021) 

. 10"8  h°m  . 

(2.5) 

£ = 1.37  • HT8  car1 


data.  (cf.  Chapter  IV,  equation  (5-3)).  We  notice  that  the  emission  gra- 
dient measured  from  1952  data  is  more  than  twenty  per  cent  greater  than 
from  1932  data.  This  implies  that  either  chromospheric  or  technical 
conditions  for  the  measurements  wore  considerably  different  at  the  two 
eclipses.  If  we  consider  this  point  together  v/ith  the  various  assump- 
tions made  in  the  analysis,  the  agreement  between  the  two  cases  is 
satisfactory.  Under  the  assumption  of  hydrostatic  equilibrium,  the  ob- 
served density  gradients  for  1952  and  1932  leads  to  the  temperatures  of 
25,000°  and  40,000°,  respectively.  We  will  return  to  this  comparison 
of  p>  in  the  two  eclipses  later. 


Since  the  Np  gradient  does  not  depend  on  the  assumed  temperature, 
we  obtain  the  expressions  for  log  Np  for  various  lower  temperatures 
simply  by  changing  the  first  term  in  equation  (2.4)  or  (2.5).  In  fact, 
we  can  write  equation  (2.4)  in  either  the  form 


or 


log  N2  = 16.37  (±0.023)  + £log( 


- 0.578  (±0.018)  x 10 h™ 

(2.6) 


log  Np  = 14.52  (±0,023)  4 ilog  Te  - 0.578  (±0.018)  x 10-8  h®1 

(2.7) 

pince  varies  aB  (Te)“^.  Three  broken  lines  in  figure  IV  repre- 

sent log  N2  given  for  T©  - 20,0000,  10,000°,  and  5000°.  If  we  know  the 
temperature  at  come  height  and  the  chromospheric  temperature  gradient  is 
not  large,  we  can  determine  the  N2  value  at  each  height  from  equation 
(2.6).  Later  in  section  5,  we  will  calculate  the  value  of  Np  at  every 
500  km  in  the  lower  region  of  the  chromosphere,  using  the  temperature 
distribution  obtained  from  the  continuum. 


The  above  argument  is  especially  important  when  w e r ecall  that  our 
final  aim  is  to  find  the  Te  and  Ne  distributions  by  successive  approximations. 
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If  tho  temperature  gradient  is  small,  the  form  of  equation  (2.6)  tail  not 
change  significantly  in  higher  order  approximations,  and  we  can  use  it  to 
calculate  N2  immediately.  Of  course,  in  a rigorous  sense,  the  empirical 
value  of  p may  be  altered  from  one  model  to  the  next,  because  the  self- 
absorption will  change  as  H2  is  changed.  However,  the  amount  of  correction 
may  be  negligibly  small  in  the  higher  regions  where  N2  values  are  small. 
Thus,  the  actual  amount  of  work  in  successive  calculations  can  be 
considerably  reduced. 
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400  1200  2000  2800  3400 

JL 

FIGURE  3E  B vs.  Height 

Now,  returning  to  the  results  summarized  in  pages  61  and  62,  we  plot 
the  coefficient  B against  height  in  Figure  V,  The  points  are  widely- 
scattered  and  it  seems  difficult  to  presume  a correlation  with  height. 
However,  if  we  again  exclude  the  four  points  corresponding  to  those  marked 
by  X in  Figure  VI,  we  see  that  for  the  region  lower  than  h =>  2000  Ion,  the 
values  of  B are  almost  constant.  Therefore,  wo  take  an  arithmetic  mean 
of  these  values  and  write  t 

log  bn  - log  bn  + k =»  0.02l(  k (2.8) 

for  the  region  lower  than  h “ 2000,  and 

log  bn  - log  bn  + k « 0.022  k (2,9) 


for  the  whole  region. 
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FIGURE  HI  Theoretical  Value  of  Log  bn  vs.  n 


As  we  have  discussed  in  the  previous  chapter,  the  only  theoretical 
values  of  bn  available  for  n higher  than  10  ore  those  that  Thonus  calcu- 
lated approximately  for  Te  = 35,000°.  These  are  plotted  in  Figure  VI. 

The  solid  line  represents  a best  fit  to  these  theoretical  values.  The 
broken  line  was  drawn  from  equation  (2.9)  with  the  constant  adjusted  so 
as  to  fit  the  first  four  points.  The  gradient  of  the  solid  line  is  0.012. 


§3. 
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In  Figure  III  of  the  previous  section  wo  plotted  the  relative 
effects  of  departure  from  thermodynamic  equilibrium  and  self-absorption 
for  each  height  and  observed  that  the  diagrams  indicated  a systematic 
change  of  these  relative  effects  with  height.  However,  the  quantity  thus 
plotted  is  relative  and  the  consistency  of  the  results  is  examined  only 
qualitatively.  In  section  3 a method  ms  discussed  for  determining  both 
the  self-absorption  and  the  departure  from  the  equilibrium  simultaneously 
under  some  simplified  assumptions.  We  now  wish  to  determine  the  absolute 
scale  of  Figure  III,  so  that  ve  can  proceed  to  calculate  more  rigorously. 


b^  must  approach  unity  as  n goes  to  infinity,  since  the  free  elec- 
trons can  be  assumed  to  obey  the  Maxwellian  distribution  at  a uniqiia 
kinetic  temperature.  Hence,  either  the  continuum  (free-bound  emission) 
intensity  or  the  emissions  in  the  lines  near  the  series  limit  may  be 
regarded  as  free  from  the  bn-effect.  At  the  some  time,  self-absorption 
should  be  negligibly  small  in  these  legions . Therefore,  using  the 
continuum  or  the  high-order  line  intensities,  we  may  calculate  the  ab- 
solute amount  of  bn  and  the  self-absorption  observational ly  so  that  we 
can  fix  the  zero  point  of  Figure  III. 


6a 


The  free-bound,  emission  in  the  Balmer  continuum  in  a frequency 
range  eltf  integrated  along  the  line  of  sight  and  above  the  moon's  limb 
may  be  w ritten  as : 


p — 2,6  3 • 

*-  r-  a. 
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£- 


'AD 
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where 
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and  where  E,  and  Nt  refer  to  the  integrated  quantities  of  Tx^and  ft 
along  llne-of-sight,  as  was  stated  on  pagers’.  A factor  Renters  through 
the  integration  over  x.  . The  integration  is  carried  under  the  same 
conditions  and  assumptions  we  made  in  Chapter  III. 


The  total  Balmer  line  intensity  is  given  by  equation  (2.1)  of 
Chapter  III.  We  replace  Nn  in  this  equation  by  the  Boltzmann-Soha 
formula  to  make  the  terms  Nj.,  Ne  and  Te  explicit.  Then  we  have  the 
following  expressions 
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E 


NiNe 


%•»> 


(3.2) 


where 


4X 


( 


xic  ™ k. ) 


7 Z '10 


28 


Taking  the  logarithmic  ratio  of  (3.2)  to  (3.1),  we  write 

A-j  EooX  = 5_.4-3'?  dp  --£z) 

B 2El  . r ~ 

E ^ e v*.  ■ bv  \ (3.3) 

where  is  the  intensity  at  the  series  limit.  Hence,  from  equation 

(3.3),  we  can  obtain  the  observed  amount  of  • Since 

the  continuum  measurements  were  considerably  more  scattered  than  the 30  of 
the  lines,  we  used  the  values  smoothed  by  least-square  fitting  for  Eo oi. 
and  the  raw  values  for  lines,  although  the  results  were  not  seriously 
different. 

Figure  VII  shows  fiw  €Y"  fc„  [««"}  at  about  every  500  kilometers, 
plotted  against  the  quantum  number.  The  shape  of  the  line  at  each  height 
1b  essentially  the  same  as  the  corresponding  one  in  Figure  III.  While,  in 
Figure  VII,  the  value  lt-&  e * " b (fri *]  is  plotted  for  every  line,  it 
should  be  noted  that  the  lines,  , H22,  H25.  H28  and  H3Q  have  been  found 
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to  be  blended  vdth  other  lines,  as  is  discussed  in  Chapter  IV.  Hence, 
these  lines  should  be  excluded  in  considering  the  general  characteristics 
of  the  diagram  that  are  discussed  below. 

The  essential  facts  of  physical  significance  in  Figure  VII  are  sys- 
tematic changes  of  e.*  " b * fflbs  ^ with  chromospheric  height  as  well  as 
vdth  line  number.  The  quantity  , knj  t*  ’s  b .*•  f A bs  » '( , hereafter  denoted 
by  Rn,  increases  upwards  by  a greater  amount  for  the  lower  members  of  the 
series  than  for  the  higher  members.  Towards  the  series  limit  it  approaces 
a constant  value  monotonically , and  for  n greater  than  26,  becomes 
slightly  greater  than  zero  (except  for  the  height,  2400  km).  ihis  verifies 
observationally  that  both  departure  from  thermodynamic  equilibrium  and 
self-absorption  are  negligible  for  these  higher  members.  If  we  draw  a 
horizontal  line  at  the  apparent  asymptote,  we  see  that  Rn  for  h lower 
than  1000  km  comes  below  this  line,  and  at  higher  regions  it  comes  above 
it.  ^ince  jLo  ^ Bbs*}  is  a negative  quantity,  decreasing  with  height,  and 
Hox>  b-r,  is  positive,  we  may  say  that  the  s elf-absorption  exceeds  bn  in 
the  dower  region,  and  the  latter  becomes  predeminent  in  the  higher 
chromosphere.  We  may  not,  however,  conclude  immediately  that  the 
increase  outwards  in  the  chromosphere  unt.il  we  attempt  by  estimation 
to  separate  and  independently,  as  will  be  done  in  the  follow- 

ing section. 


t Before  proceeding  to  the  bn  determination,  it  will  be  interesting 

to  calculate  the  Rn-height  relation  using  the  intensity  of  higher  order 
lines  instead  of  the  continuum.  As  we  have  discussed,  the  continuum 
measurements  are  considerably  scattered  and  we  have  used  the  smoothed 
values  for  the  present  analysis.  We  have  shown  that  bj^  and  terms 

become  negligible  for  the  highest  lines  of  the  series 'both  theoretically 
and  observationally.  Hence  we  can  simply  calculate  the  difference  be- 
tween the  line  intensity  of  a given  line  and  one  of  the  highest  order  lines 
in  order  to  determine  Rn.  From  equation  (3.3)>  vie  may  write  immediately 

h eY'-  ^ ^ - h - h 

(3.4) 

in  which  n'  refers  to  the  line  whose  self-absorption  and  departure  from 
thermodynamic  equilibrium  are  assumed  to  be  negligible.  We  take  the 
mean  value  of  the  intensity  of  H^i,  Hgo,  and  H27  (H30,  H28  are  blended) 
as  E’rr'i.  and  set  n*  = 29.  Figure  VIII  shows  the  result  thus  obtained 
for  heights  corresponding  to  those  taken  in  Figure  VII,  Ihe  increase  of 
Rn  with  the  chromospheric  height  is  not  as  great  as  when  vie  used  the 
. continuum,  but  the  systematic  change  with  height  remains. 

Vfe  have  computed  Rn  for  all  heights  and  obtained  a gradual  increase 
with  height.  For  each  individual  n,  variation  of  the  line  intensity  by 
• a fraction  of  the  observational  error  would  shift  Rn  either  positively 

or  negatively.  In  the  latter  case  where  we  used  high  order  lines  for  a 
reference,  most  anomalies  in  Rn  appear  for  heights  higher  lhan  1930  km. 

Rn  at  these  higher  regions  is  smeller  than  for  the  next  lower  height, 


Log(bne*nAbs.) 

FIGURE  301  Log(bn-eXn-Abs.)vs.  n (a) 
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FIGURE  IVULL  Log  (bn-eXn-Abs.)  vs.  n (b) 


which  was  not  shown  in  the  results  when  we  used  the  continuum  as  a 
reference  point.  We  may  interpret  this  discrepancy  as  due  primarily  to 
the  inaccuracy  of  the  measurements  of  ’Ey/*.  values  for  greater  heights, 
where  these  lines  are  very  much  weakened.  Since  the  smoothed  values  of 
the  continuum  are  used  in  the  former  case,  such  discrepancies  were 
corrected. 

Except  at  the  greatest  heights,  the  only  anomalies  in  Rn  appear  at 
h = 960  km.  Here  Rn  seems  abnormally  high,  and  the  Rn-curve  has  a rather 
peculiar  shape.  In  order  to  minimise  such  effects  we  constructed  Rn 
diagrams  hi'  taking  mean  values  of  Rn  at  several  succeeding  heights,  except 
for  960  km  and  heights  higher  than  1930  km,  in  the  foli owing  way.  We  took 
the  mean  of  the  Rn  at  530,  640,  750,  and  850  lan  and  plotted  the  value  for 
the  mean  height,  h = 685  km.  Similarly,  we  chose  combinations  of  1070, 
1180,  1280,  1390  lan  and  1500,  1600,  1820  lan  to  determine  the  mean  values 
for  h = 1230  and  1640  km,  respectively.  Thus  we  obtained  an  Rp  curve  for 
three  heights  shown  in  Figure  IX.  'ihe  deviate  ns  of  points  for  individual 
heights  are  smoothed  by  this  averaging  procedure,  and  all  the  points  lie 
well  within  a clearly  defined  Rn  curve.  Further,  the  height  at  which  the 
slope  of  the  Rn  curve  changes  sign  for  n given  n,  more  clearly  indicated 
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by  Figure  IX.  Ihe  general  characteristics  of  the  effects  of  bn  and  self- 
absorption,  concluded  from  the  systematic  change  of  the  Rn  curves  in 
Figure  VII,  may  be  explained  more  conclusively  in  Figure  IX. 

In  the  preceding  discussions  about  R-n  a Srj.  e * * • b „ 1 ^ , we 

implicitly  ignored  the  term  g**  , which  depends  on  the  temperature.  If 
the  chromosphere  is  assumed  to  be  isothermal,  the  qualitative  interpreta- 
tions should  not  be  altered  at  all.  On  Ihe  other  hand,  if  there  is  a 
temperature  gradient  the  e*"*  term  may  affect  slightly  the  shape  of  the 
Rjj  curves.  However,  e**  has  a very  small  value,  even  for  low  temperatures, 
and  decrease  exponentially  with  temperature.  We  shall  come  back  to  this 
question  again  in  the  next  section,  but  in  lable  XIV,  we  show  the  values 
lrj.-e><V|  computed  for  various  Te  to  compare  the  order  of  magnitude  with 
Rn  in  Figures  VII,  VIII,  and  IX. 

Table  XIV. 


J Eofr  ex» 


Tb\n 

10 

11 

12 

13 

15 

17 

20 

25 

30 

5,000 

0.14 

0.11 

0.09 

0.08 

0.06 

0.05 

0.03 

0.02 

0.02 

10,000 

.07 

.06 

.05 

.04 

.03 

.02 

.02 

.01 

.01 

20,000 

.03 

.03 

.02 

.02 

.02 

..01 

.01 

.01 

.01 

30,000 

0.02 

.02 

.02 

.01 

.01 

.01 

.01 

.00 

.00 

We  see  in  the  above  table  ■e,*v'  is  negligible  for  Te  higher  than 
20,0000  and  is  appreciable  only  for  lower  order  Balmer  lines  at  lower 
temperatures.  In  the  next  section  we  will  include  this  factor  for  an 
appropriate  temperature  to  determing  the  observed  value  of  bn. 


§4.  .fPA t^g  q tfre.obperye&Jfo. 

In  the  previous  section  we  calculated  the  observed  values  of  Rjj  that 
are  shown  in  Figures  VII,  VIII,  and  IX.  To  eliminate  -e.**  from  Rn,  we 
must  know  the  temperature  distribution  in  the  chromosphere.  But,  the 
contribution  of  this  term  is  very  small  and  we  can  treat  it  as  a secondary 
correction.  If  we  once  obtain  bn  or  , wc  may  estimate  the  tem- 

perature d istribution,  after  which  wc  can  eliminate  e.*Y>  and  repeat  the 
process  successively.  A larger  contribution  to  Rn  comes  from  b.*- fAbj-J. 
Either  bn  or  ^ AV>s,,\  can  be  determined  if  we  know  the  other.  If  we  have  a ' 
knowledge  of  the  chromospheric  temperature  and  electron  density  (the  latter 
is  much  less  sensitive  to  bn  than  the  former),  we  can  calculate  bn  theore- 
tically, then  determine  ; or  we  can  obtain  N2  from  b2,  Te  and  Ne, 

then  compute  the  theoretical  value  of  fftbs^to  compare  with  the  observed 
» At  present  we  have  only  a value  lor  Te  derived  from  ihe 
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continuum.  Vie  wish,  however,  to  yet  a mode]  from  the  line  intensity 
IndejAinctcnt  of  the  continuum.  In  section  2 we  developed  a Method  to  got 
< observations lly  at  each  height.  To  derive  I?2  i'rom  this  result,  we 
have  to  know  the  temperature  at  each  height,  since  «<  in  this  case  is  a 
function  ol  T'e.  however, was  obtained  try  a purely  empirical  method 
and  {AbsjJ  is  a function  of  NjfC  , so  we  can  calculate  without  assuming 

the  temperature,  Then  we  can  obtain  the  observed  values  of  bjj,  choosing 
the  parameter  e to  give  a best  fit  with  the  theoretical  bn  computed  for 
various  temperatures.  Thus,  vie  can  determine  the  temperature  empirically 
and  repeat  Ihc-  procedure  to  reach  a final  consistent  model  of  the  chromosphere. 

The  values  of  N^zn,o  tabulated  in  page  61  to  page  62,  were  used  to 
estimate  ^bsn”l  at  each  height.  Then,  by  subtracting  this  from  Rn>  ob- 
tained from  the  line  and  continuum  intensities,  (cf.  figure  VII ),  we 
obtained  the  values  of  , that  are  plotted  against  n in 

Figure  X.  The  values  corresponding  to  the  blended  lines,  bj^,  b2g,  and 
b30  were  omitted.  The  dif.fei'ent  lines  used  to  distinguish  each  height 
are  the  same  as  vie  used  in  figure  Vll.  For  each  height,  vie  have  included 
in  our  calculations  only  the  lines  for  which  was  less  than  1. 

(For  example,  for  the  lowest  height,  n<15  was  not  included).  To  estimate- 
log  bn,  wc  subtracted  loge^n  computed  for  the  appropriate  temferature. 

Since  e*n  is  negligibly  small  for  Te  20,000°  and  n <.  10,  this  diagram 
itself  gives  log  bn  directly  if  the  chromospheric  temperature  is  higher 
than  20,000°.  At  first  glance,  ve  see  that  bn  is  practically  the  same  for 
all  heights  except  the  highest.  This  implies  that  -the  chromosphere  is 
isothermal  at  least  in  the  region  lower  than  2000  km.  We  will  discuss 
this  point  later. 


If  we  assume  that  the  chromospheric  temperature  is  as  low  as  5000°, 
we  must  subtract  e*n  computed  in  Table  XIV  to  estimate  the  observed 
log  bn  corresponding  to  5000°.  We  see  from  Table  XIV  that  the  log  bn 
lines  will  be  lowered  by  a maximum  amount  of  0.14,  at  n = 10.  If,  however, 
Te  = 5000°,  log  bn  should  be  equal  to  zero,  since  the  radiation  temperature 
of  the  chromosphere  is  nearly  5000°,  giving  a iwltzmann  distribution,  with 
bn  = 1.  Thus  we  estimate  the  real  temperature  to  be  somewhat  higher  than 
the  photospheric  radiation  temperature. 


The  apparent  shapes  of  the  log  bn  versus  n curves  seem  to  be  nearly 
straight,  with  perhaps  a slight  rise  toviurds  the  lower  order  n.  There- 
fore, vie  take  a mean  of  log  bn  obtained  for  four  different  heights,  530, 
1070,  1500,  and  2040  km  to  determine  the  value  of  log  by,  at  each  n.  Then, 
an  analytical  expression  of  log  bn  as  a function  of  n may  be  obtained  by 
fitting  these  mean  values  with  either  of  the  following  two  forms: 

log  bn  = -0.13  (±0.022)  + 6.7  (±0.28)  • £ (4.1) 

or 

log  bn  = 0.73  (±0.049)  - 0.024  (±0.0011)  • n . (4.2) 

The  form  (4.1)  was  used  for  fitting  the  theoretical  computations  for  Te  = 
35,000°  given  in  Chapter  TV,  and  the  linear  form  (4.2)  was  used  for  the 


,I29(±0.022)+ “6.71  (±0.28) 
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Log  bn 

FIGURE  X Log  bn  vs.  n (a) 
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determination  of  in  section  2.  Prom  the  probable  errors  shown,  we 

see  that  log  bn  is  well  represented  by  (4.2).  In  section  2.  we  obtained 
a value  of  B,  corresponding  to  the  coefficient  of  n in  (4.2),  of  0.024  or 
0.022.  This  remarkably  good  agreement  may  be  regarded  as  a check  on  the 
consistency  of  both  computations. 


FIGURE  XE  Log  bn  vs.  n (b) 


The  approximate  theoretical  values  of  log  bn  for  Te  = 35,000°  are 
indicated  byX  in  Figure  X.  This  shows  that  the  observed  bjj  indicate  a 
chromospheric  electron  temperature  somewhat  lower  than  35,000°.  To  get 
further  information  about  the  exact  value  of  the  temperature,  we  have  to 
compute  bn  for  several  ranges  of  Te  lower  than  35,000°. 

Figure  XI  is  a similar  diagram  of  log  bn  derived  from  Figure  IX.  Rn 
in  this  case  was  obtained  by  taking  the  intensity  of  Balmor  lines  relative 
to  the  higher  lines  instead  of  to  the  continuum.  It  was  then  averaged 
for  four  successive  heights  in  the  same  way  that  we  averaged  the  corres- 
ponding results  for  the  continuum  comparison.  If  our  assumption  of  negli- 
gible $ Abs  y,  \ and  byj  for  the  lines  higher  than  n = 27  as  well  as  the 
continuum  is  true,  then  both  figure  X and  XI  should  be  identical  except 
for  a slight  difference  of  height  represented  on  the  diagrams.  We  ob- 
tained the  following  algebraic  expressions  for  log  bn'  in  figure  XI: 
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(U.3) 


log  bn  = -0.22  (+  0.016)  + 7.1  (+  0.20)  * | 
log  bn  = 0.70  (+  0.097)  - 0.027  (+  0.0022)  • n.  (h,k) 


The  values  represented  by  the  above  formulae  are  slightly  smaller  than 
those  in  (U.l)  and  (lu2). 

Finally,  \re  shall  consider  the  following  computations  as  a rough 
estimation  of  the  contribution  of  self-absorption  to  the  Rn  curves.  In 
Chapter  IV,  we  obtained  N2  from  1932  eclipse  data  from  a different  method. 

As  we  have  discussed  previ ously,  various  observed  results  indicate  that 
the  chromosphere  in  1932  was  in  a considerably  different  condition.  There- 
fore one  may  not,  in  a rigorous  sense,  apply  the  results  of  1932  to  1932, 
However,  it  may  be  interesting  to  estimate  log  bn  using  N2  obtained  from 
1932,  and  compare  the  results  to  those  discussed  above.  Figure  XII  shows 
the  diagram  corresponding  to  Figure  X,  thus  obtained.  The  mean  log  bn 
curve  in  Figure  XII  approximately  coincides  with  the  mean  curve  in  Figure  X. 
This  is  a rather  unexpected  result  if  we  take  into  account  the  considerable 
difference  in  N2  values  between  the  two  eclipses,  as  discussed  in  section  2. 
It  may  be  partly  due  to  the  small  contribution  of  self-absorption  to  Rn  for 
higher  order  lines. 


§3.  Comparison  to  the  model  obtained  from  the  Balmer  continuum,  and  the 
determination  of  bp. 

In  the  preceding  sections,  it  was  pointed  out  that  we  can  determine  the 
values  of  Te  and  ne  independently  from  the  emission  in  the  continuum.  The 
analysis  of  the  continuum  data  from  the  1932  eclipse  was  carried  out  by 
Athay,  Pecker,  and  Thomas  (193k).  From  the  intensity  of  the  Balmer  free- 
bound  emission  relative  to  the  photospheric  scattered  light  and  emission 
from  the  negative  hydrogen  ion,  at  the  Balmer  series  limit,  they  determined 
the  values  of  Te  and  n„  at  each  height.  Table  XV  shows  their  results  for 
the  region  lower  than  n = 3000  Ion, 

Table  XV 


Model  of  chromosphere 

obtained 

from  the 

Balmer  continuums 

Height  (km) 

500 

1000 

1300 

2000 

2300 

3000 

Te 

3,3oo 

6,000 

6,300 

7,100 

7,800 

8,600 

log  ne 

11.63 

11. kk 

11.21 

11.02 

10.82 

10.61 

-:f-From  Table  VI  in  the  forthcoming  paper  by  Athay,  Pecker,  and  Thomas. 

We  see  in  the  above  table  that  the  electron  temperature  is  slightly 
higher  than  the  photospheric  radiation  temperature  and  increases  slowly 
upward , According  to  the  method  discussed  in  section  2,  we  can  calculate 
from  equation  (2,6)  the  values  of  ^(h)  at  each  height  corresponding  to  the 
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temperature  given  in  Table  XV.  The  values  of  N2(h)  thus  obtained  are 
shown  in  Table  XVI  in  comparison  witn  the  values  corresponding  to 
Tq  - 6,500°. 

Table  XVI 


Values  of  N2(h) 


Height  (km) 

500 

1000 

i5oo 

2000 

2500 

3000 

log  N2(h) 
for  T0  in  Table 

16.10 

XV 

I5.8ii 

15.56 

15.29 

15.02 

1U.75 

log  l'!o(h) 
for  Tq  - 6,500° 

16.114 

15.06 

15.56 

15.2? 

114.98 

1I+.69 

The  differences 

between  N2(h)  for 

the  two  cases 

are  very 

small. 

Hence  the 

assumption  of  constant  temperature  does  not  necessarily  lead  to  an  appreci- 
able error,  at  least  for  the  determination  of  N2(h). 

In  the  previous  section,  we  calculated  the  observational  values  of 
bn.  These  values  load  to  one  determination  of  electron  temperature  and 
electron  density.  If  we  have  the  theoretical  value  of  bn  for  small  inter- 
vals of  Te  and  Ne,  we  can  obtain  Te  and  Ne  from  these  results  and  compare 
with  those  in  Table  XV.  At  tho  present,  we  have  concluded  only  that  T8 
estimated  from  the  line  intensity  is  lower  than  35,000°  and  may  be 
constant  throughout  tho  lower  region  of  the  chromosphere.  The  values  of 
Ta  in  Table  XV  seem  to  agree  with  this  conclusion. 

We  now  proceed  to  calculate  values  of  bp  from  the  Te,  ne,  and  N2 
values  In  Tables  XV  and  XVI.  From  equation  (2,8)  in  Chapter  III,  the 
number  of  hydrogen  atoms  per  cubic  centimeter  that  are  in  the  second  level 
at  chromospheric  height  h may  be  written 

ItO  . (5.D 


In  the  present  case,  we  may  transform  from  n2(h)  to  N2(h),  the  total 
number  of  hydrogen  atoms  in  the  second  level  within  a column  of  unit 
cross-section  along  the  line  of  sight  by  tlie  following  formulat 


N,d-o 


1 

0 


(5.2) 


;'.o 

Q 

in  which  we  set  =•  1,33  • 10~°cm"’x  according  to  the  result  obtained  In 
equation  (2.lj). 

Hence,  after  substituting  numerical  values,  we  have  from  equations 
(5.1)  and  (5.2)  ' '■  ' 

~/SL 

Hx C<)  = M0G  l^C-O  T^)  (5.3)  I 

* : ; g 

•4 

where  we  introduced  the  usual  approximately,  n^  a rie,  in  the  chromosphere. 

Taking  the  logarithms  of  equation  (5.3),  we  may  write  bg(h)  in  the  following 
formt 


tf 


hj  - S’.o  2.  + 


/,  nnr  ■ / o 1 
Te.Ct) 


By  substituting  the  values  Np(h)  given  in  Table  XVI  and  ne(h)  and 
T0(h)  in  Table  XV,  we  can  obtain  the  empirical  value  of  log  b2  at  each 
height.  The  results  are  shown  in  Table  XVII.  We  see  that  log  b2  increases 
considerably  with  height.  By  comparing  the  values  of  log  b2  thus  obtained 


Table  XVII 


Observed  values  of  log  b2 


Height  (km) 

5oo 

1000  1500  2000 

2500 

3000 

Log  b2 

Of  3.? 

0,79  1,0  24  1,63 

2,03 

2,45 

0.91* 

0.90  0.78  0.80 

0.8U 

with  the  theoretical  values  computed  for  small  interval  of  Te  and  ne,  we 
can  examine  the  model  derived  from  the  continuum  for  internal  consistency. 
For  constant  values  of  n^,  bn  is  larger  for  higher  Te  than  for  lower  Te. 
Since  bn  is  a function  of  both  Te  and  r^,  however,  we  may  not  conclude 
immediately  that  the  increase  of  bn  with  height,  shown  in  Table  XVII, 
indicates  the  increase  of  Te  upwards.  However,  theoretical  calculations 
by  Thomas  (19k9a)  have  shown  that  bn  is  not  as  sensitive  to  a change  of 
Ne  as  it  i3  to  a change  of  Te,  Tho  above  table  shows  that  bn  increases 
by  a factor  of  100  over  a height  range  of  2,500  km.  Hence,  we  may  conclude 
that  a temperature  increasing  upwards  is  implied. 

We  can  prove  the  internal  consistency  of  the  results  given  in  Table 
XVII  in  the  following  manner: 
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Taking  logarithms  of  equation  (3.1)  Chapter  V,  and  differentiating  with 
respect  to  h,  we  have  the  height  gradient  of  the  [turner  free-bound 
emission  at  the  seri.es  limit 

3 ol^Te 

4il  " AX  ^ ^ ( 


Vfe  have  ignored  the  term  - dfin§/dh  from  the  above  equation,  since  we  showed 
earlier  that  p is  constant,  at  least  for  the  region  lower  than  3,000  km 
in  the  chromosphere. 

Similarly,  from  equation  (5.1),  we  can  write 


AlN.'Ne 
~~  ""  * 

&JIm.  j)  2. 

_ Xn 

3 *erf 

cA^.  4 X 

AX 

2.  AH 

in  which  we  note  that 

d ^ ^ 

k £,N;Ne 

/ 

AX 

for  d^np/dh  = 0,  as  we  mentioned  above.  Substituting  equation  (5.6) 
into  (5.5),  we  have 


d-t'E.o  _ Atk  imW'  j% 

~~  A$  <*n  + T<t  a-s. 


Hence,  using  the  above  equation,  we  can  calculate  the  bg-height  gradient 
from  the  emission-height  gradient  in  the  continuum.  Substituting  the 
following  values  obtained  from  the  1952  data: 


d 

A fL 


-S' 

2->  l°l  ‘ l 0 j 


djkNa 

AX 


— 1,3V  10 


and  the  values  of  Ta  shown  in  Table  XV,  we  obtain  the  b2-height  gradients. 
These  values  of  dlogb^dh  are  shown  in  Tablo  XVIII,  Comparing  these 
results  with  those  in  Table  XVII,  we  see  that  the  agreement  of  b2~height 
gradient  is  satisfactory.  This  may  prove  that  the  present  determination 
of  the  values  of  N2  are  consistent  with  the  results  obtained  from  the 
continuum. 
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Table  XVIII 


Value  of  d l0£  b_2 
cIH 


Height  (Ion) 

500 

1000 

1500 

2000 

2500 

3000 

d lofr  b'2  . lo^cm"1 
dh 

0.62 

0.58 

0.60 

0.55 

0.52 

0.62 

§6,  Discussion  of  the  comparison  of  the  density-height  gradient  with  the 

emission  gradient. 

In  Chapter  I it  was  mentioned  that  the  emission-height  gradient  of  the 
flash  spectrum  has  been  considered  to  be  one  of  the  anomalous  features 
that  can  not  be  explained  by  a classical  model  of  the  chromosphere#  All 
the  earlier  eclipse  observations  have  shown  that  the  emission-height 
gradient  of  the  Bnlmer  lines  increase  considerably  towards  the  higher  order 
lines.  From  the  1932  eclipse,  Athay  and  others  havo  found  that  the 
emission  hoight  gradient,  an,  increases  monotonically  from  1.56  x 10“®cm~^- 
for  Hxo  to  the  value  of  2.12  x 10-8cm-l  for  H31.  The  value  of 
2.19  x 10-8cm~l  is  obtained  for  the  height  gradient  of  the  continuum  at 
the  Balmor  series  limit. 

In  the  preceding  sections,  we  have  mentioned  several  pieces  of 
evidence  indicating  that  the  assumptions  of  constant  temperature  and 
exponential  decrease  of  hydrogen  atoms  are  justifiable  as  a first  approxi- 
mation, in  the  lower  region  of  the  chromosphere.  Under  such  assumptions, 
the  integrated  emission  gradient  determined  from  slitless  flash  spectra 
is  equal  to  the  height  gradient  cf  emission  per  unit  volume.  Therefore, 
if  we  ignore  the  effect  of  self-absorption  as  in  the  earlier  analysis, 
the  emission  gradient  of  the  line  En2  must  be  equal  to  the  density  gradient 
of  the  hydrogen  atoms  in  the  n^h  level.  Further,  if  departures  from 
thermodynamic  equilibrium  are  ignored,  the  density  gradients  of  all  the 
energy  states  are  the  same.  Thus,  in  such  a classical  model,  the  emission 
height  gradients  for  all  lines  are  identical  also. 

In  section  5 of  Chapter  IV,  the  value  of  0,86  x 10“®cm“-^  was  obtained 
from  the  1932  eclipse  for  P2,  the  density  gradient  of  hydrogen  2nd  state# 

In  section  2 of  this  chapter,  p?  a 1«33  x 10“®cm~l  was  obtained  from  the 
1952  data.  These  values  are  considerably  lower  tiian  the  emission-height 
gradient  determined  from  the  higher  order  lines.  This  disagreement  is 
considered  to  be  duo  to  the  effects  of  departures  from  thermodynamic 
equilibrium.  As  was  discussed  in  Chapter  II,  departures  from  thermodynamic 
equilibrium  decrease  towards  the  higher  members  of  the  series  and  are 
negligible  at  the  series  limit.  At  the  same  time,  the  self-absorption  for 
the  highest  line  is  negligibly  small  compared  to  that  for  lower  lines. 

Honce,  from  the  above  considerations,  the  emission  gradient,  aw=  2.19,  at 
the  series  limit  may  be  considered  as  the  density  gradient  of  nydrogen 
in  the  highest  level.  Therefore,  we  see  that  there  is  considerable 
difference  in  our  present  results  between  P2  and  p«. 


la  the  present  thesis,  the  density-height  gradients  were  assumed  to  be 
the  same  for  all  levels.  Since  our  results  seem  Lo  show  different  values 
of  j for  different  n,  the  assumption  of  a constant  gradient  is  not 
strictly  justified.  We  have  to  take  the  dependence  of  pn  upon  n into 
account,  From  equation  (2.7)  in  Chapter  III,  we  see  that  the  emission  in 
the  Calmer  lines  is  a function  of  nn  and  n£.  Hence,  if  we  include  the 
depen  tence  of  pn  on  n,  both  pn  and  p n enter  into  the  integration  over  the 
volume.  Therefore,  it  is  impossible  to  perform  the  integration  analytically 
as  was  done  in  Chapter  ITT,  and  it  must  bo  clone  numerically.  Further, 
since  pn  depends  on  both  bn  and  the  self-abosrption,  it  is  difficult  to 
derive  its  value  for  every  line.  It  seems  that  the  only  way  to  carry  out 
the  numerical  integration  is  by  a trial  and  error  method  of  successive 
approxima ti on . 

At  the  present,  it  is  more  difficult  to  estimate  numerically  the 
errors  brought  into  the  results  by  the  assumption  of  p constant  with  n, 
than  to  estimate  the  errors  due  to  the  assumption  regarding  line  profiles 
- namely,  we  assumed  that  the  central  intensity  of  each  line  is  solely 
determined  by  the  amount  of  Doppler  broadening.  In  the  expression  for 
the  self -absorption  derived  in  Chapter  III,  seif -absorption  is  a function 
of  tie  atomic  absorption  coefficient  at  the  line  center.  Hence,  it  might 
be  necessary  to  consider  Stark  broadening,  particularly  since  the  electron 
temperature  was  found  to  be  less  than  30, OCX)0  in  the  lower  region  of  the 
chromosphere.  Further,  in  a rigorous  calculation,  one  has  to  consider 
the  optical  depth  at  the  line  center  as  a transfer  problem.  Thus,  the 
problem  ultimately  becomes  much  more  difficult. 

One  should  note,  however,  that  the  present  analysis  was  based  on  the 
relative  intensities  of  the  lines,  considering,  in  particular,  the  relative 
amount  of  self-absorption  and  bn.  Therefore,  the  aggregate  affect  of  all 
the  small  errors  may  not  be  appreciable  in  the  final  results. 


§7 . Concluding  remarks. 


Considerable  differences  between  the  Balmer  line  intensities  observed 
at  the  1932  and  1952  eclipses  appeared.  In  the  lower  region  of  the  chromo- 
sphere, the  absolute  values  of  line  intensities  of  the  1952  spectra  are 
about  five  times  as  large  as  those  of  1?32.  The  emission-height  gradients 
of  Balmer  lines  are  also  approximately  Uo#  larger  in  1952  than  in  1932 
data.  The  theoretical  investigation  treated  in  Chapter  IV  and  Chapter  V 
indicated  that  these  differences  are  due  mainly  either  to  an  essential 
variation  of  the  physical  conditions  in  tiie  chromo sphere  from  time  to  time 
or  to  inhomogenie ty,  i.e.,  absence  of  spherical  symmetry  in  the  chromosphere. 
Wq  shall  discuss  these  points  by  comparing  the  numerical  results  of  this 
chapter  with  those  of  Chapter  IV. 

In  Chapter  IV,  it  was  concluded  that  the  departure  from  thermodynamic 
equilibrium  was  not  a predominant  effect  in  the  1932  eclipse  observations, 
and  the  difficulties  that  appeared  in  using  the  classical  theory  to  explain 
the  low  Balmer  decrement  could  be  solved  by  including  the  effects  of  self- 
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absorption  only.  On  the  other  hand,  for  the  1952  data  our  analysis 
showed  that  the  two  effects  are  equally  important  and  the  values  of  bn 
obtained  are  nob  at  all  negligible.  The  calculations  in  section  3 of 
Chapter  IV  for  determining  bn  and  the 'amount  of  self-absorption  simul- 
taneously also  agree  with  this  conclusion.  The  successive  computations 
baset.  on  the' least  square  principle  did  not  converge  in  the  analysis  of 
the  1932  data,  whereas  for  the  1952  data  it  gave  consistent  solutions 
for  bn  and  self-ab sorption.  Hence,  the  failure  in  applying  the  method  to 
the  1932  observations  may  be  the  negligibly  small  va  lues  of  bn#  As  we 
have  discussed  in  section  2,  this  difference  is  clearly  seen  in  the 
characteristic  curve  of  Balmer  decrements  shovai  in  Figure  III. 


Considerable  variations  were  also  apparent  in  the  values  of  the 
density-height  gradient,  p.  The  value  of  0.83  x 10“8cm-l  was  given  for 
(1932),  whereas  pp  (1952)  wa3  1.33  x 10"°cm~l.  Since  these  values  of 
Pji  are  independent  of  Te,  and  all  the  other  assumptions  are  the  same  for 
both  analyses,  the  differences  in  Pg  should  be  interpreted  as  an  actual 
variation  of  chromospheric  conditions.  All  the  above  measurements  and 
comparisons  were  based  on  the  assumption  of  a spherically  symmetric  dis- 
tribution of  physical  parameters  in  the  chromosphere. 


As  we  discussed  in  Chapter  I,  Athay,  Evans,  and  Roberts  (1953)  have 
found  a locally  intensified  region  in  the  1952  spectra  taken  at  the  west 
[Limb  of  the  sun.  Their  preliminary  study  shows  that  the  hydrogen  emission 
gradient  in  this  active  region  is  only  one  third  as  high  as  at  the  east 
limb,  although  there  appears  to  be  an  exponential  decrease  of  intensity 
with  height.  The  Balmer  decrements  in  this  region  are  larger  than  in  the 
east  limb  region;  and  the  observed  decrements  from  H<?  to  H08  exceed  the 
thermodynamic  equilibrium  decrements  at  all  heights.  Similar  over- 
intensification of  the  helium  lines  and  coronal  yellow  lino  are  noted  in 
the  same  west  limb  region. 


,:r 

\t£. 


They  suggest  that  these  observations  indicate  a local  increase  of 
electron  density  and  electron  temperature  so  that  severe  departures  from 
thermodynamic  equilibrium  have  taken  place.  They  report  that  two  flares 
were  observed  at  this  region  on  the  day  of  the  eclipse,  and  that  the 
observations  suggest  a moderate  limb  flare  in  some  stage  of  its  develop- 
ment. Furthermore,  the  west  limb  spectra  show  other  irregularities.  All 
the  chromospheric  lines  are  absent  on  one  side  of  the  region  of  intensified 
emission. 


It  seems  that  these  observations  indicate  evidence  of  an  irregular 
structure  of  the  chromosphere  due  to  surface  activity  of  the  sun.  There 
may  be  finer  irregularities,  such  as  spicule  configurations,  which  are  not 
detectable  to  thin  the  limit  of  resolution  of  the  spectra.  Thus,  the 
chromosphere  might  have  to  be  considered  as  a mixture  of  disturbed  regions, 
as  has  been  suggested  by  Hagen,  (1953). 

Another  possible  cause  of  disagreement  of  the  results  from  one  eclipse 
to  the  other,  may  be  some  technical  error  in  the  measurements.  Particu- 
larly, the  height  determination  based  on  the  moon’s  profile  may  admit 
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considerable  errors. 

It  is  interesting  to  notG  that  in  one  quantity  the  data  from  the 
1932  and  1952  eclipses  agree  quite  veil,  compared  to  the  large  dis- 
crepancies between  absolute  values  of  emission-height  gradients  and 
for  the  two  eclipses.  That  quantity  is  the  difference  between  continuum 
height  gradient  and  It  is  an  important  quantity,  since,  as  dis- 

cussed in  section  it  gives  a measure  of  the  b2-height  gradient. 


REFEiLENCES 


Alfv^n,  H. 

19b7 

1950 

M.N.  101,  281. 

Cosmical  Electrodynamics  (Oxford  Clarendon 
Press) . 

Allen,  C.  W. 

1937 

19U1 

19U9 

Ap.  J,  85,  165. 
M.N.  107,  k26. 
M.N.  109,  3U3. 

Athay,  R.  Q.,  Billings, 
D.  E.,-  Evans,  J.  W. 
and  Roberts,  W.  0. 

195*4 

Ap,  J,  (in  press). 

Athay,  R.  0.,  Evans, 

J,  W.  and  Roberta, 

W.  0. 

1953 

The  Observatory,  73,  2144. 

Athay,  R.  0.,  Pecker, 
J.C.  and  Thomas, 

R.  N. 

195b 

Private  conmmni  cation. 

Baumbach,  S, 

1937 

A.N.  263,  121. 

Bell,  B. 

1951 

Special  Report,  No.  35,  Harvard  University. 

Biermann,  L. 

19W 

Zs.  f.  Ap.  25,  161. 

Chamberlain,  J.  W. 

1953 

Ap.  J.  117,  387. 

Cllli6,  G.  0.  and 
Mensel,  D.  H. 

1935 

Harvard  Circular  UlO. 

Davidson,  C.  R.  and 
Stratton,  F.J.M, 

1927 

Mem.  R.  Astr.  Soc.  6I1,  Part  I4. 

Edels,  H.  and 

Craggs,  J.  D. 

1951 

Proc.  Phys.  Soc,  A,  LXIV,  562. 

Giovanelli,  R.  G. 

19U8a 

19U8b 

19U8c 

19U8d 

Australian  J.  Sci.  Res.,  A,  1,  275. 
Ibid.,  p.  289. 

Ibid.,  p,  305. 

Ibid.,  p.  360. 

Goldberg,  L« 

1939 

Ap.  J.  8£,  673. 

Hagen,  J.  P. 

1951 

1953 

Ap.  J.  113,  5U7. 

Paper  presented  at  A.A.S.  Meeting,  August 

1953. 


87 


Houtgast,  J. 

Hoy!  , F. 

Hulst,  H.  c. 

^ . * _ — “•» 

Koelbloed,  D, 
Veltman, 

Inglia,  D.  R, 
Teller,  I 

Matsushima,  ! 

Matsushima,  i 
Miyamoto 

Menzel,  D.  H, 

Menzel,  D.  H, 
Bell,  B. 

Menzel,  D.  H, 
Cilli6,  ( 

Menzel,  B.  H 
others 

Menzel,  D.  H 
Pekeris, 

Minnaert,  M. 


1952 

Convegno  Volta,  p,  33 • 

19U9 

Some  Recent  Researches  in  Solar  Physics, 
(Cambridge:  At  the  University  Press). 

van  de 

19U6 

B.A.N.  10,  79 

19U7 

AP.  J.  105,  U71. 

1953 

The  Sun,  Chapter  5,  The  Chromosphere  and 
~ the  Corona  (Chicago:  The  University 

• - 

of  Chicago  Press), 

i and 

w. 

1951 

Proc,  K.  Nederl.  Akad.  v,  Wetensch.  5U,  U68. 

, and 

2. 

1939 

Ap.  J.  90,  U39. 

3. 

19U8 

Mem,  Ap,  1,  27. 

.1952 

Ap.  J.  Ilf,  5Mu 

3.  and 

• ' * * • .*  - 

, s. 

19U7 

Mem,  Col.  Sci.  Kyoto  Univ.  25,  63. 

» 

1931 

Pub,  Lick  Ob s.  17,  Part  1. 

. and 

19U8 

Tech.  Report,  No,  3,  ONR  Project  M 720-5. 

. and 

3.  a. 

1937 

Ap,  J,  85,  88, 

. and 

1937-U5  A series  of  papers  in  Ap.  J.  (All  references 
are  given  in  last  paper,  Ap,  J,  102,  239) • 

. and 

C.  L. 

1935 

M.II.  96,  77. 

1953 

The  Sun,  Chapter  3,  The  Photosphere. 

(TThicago:  The  University  of  Chicago 
Press). 

19U7 

Mem.  Col.  Sci.  Kyoto  Univ.  25,  31. 

19U9a 

Publ.  Astr.  Soc.  Japan,  1,  10. 

19U9b 

Ibid.,  p.  lU. 

1951a 

Ibid.,  2,  102. 

1951b 

Ibid.,  p.  113. 

1951c 

Ibid.,  3.  61. 

195ld 

Ibid,,  p.  6?. 

1953 

Zs.  f.  Ap.  31,  282. 

Miyamoto,  S 


88 


Miyamoto,  S.  and 
Kawaguchi,  I. 

1950 

Publ.  Astr.  Soc.  Japan,  1,  ilk. 

Pannokoek,  A.  and 
Kinnaert,  M. 

1928 

Verb.  Akad.  Wetensch.  Amsterdam,  13,  No.  5. 

Pecker,  J.  C, 

1951 

Ann.  d'Ap,  lk,  152. 

Pierce,  K»  and 
Goldberg,  L« 

19k8a 

I9k6b 

Tech.  Report  No.  1.  0NR  Project  M 720-5. 
Ibid.  No.  2. 

Redmann,  R.  0. 

19k  2a 
19k2b 
1952 

M.N.  102,  13k. 

Ibid.,  o.  IkO. 

The  Observatory,  7j2,  8?0. 

Roberts,  V/.  0, 

19U5 

Ap.  J.  101,  136. 

Schwarzschild,  M. 

19U8 

Ap.  J,  107,  1. 

Shane,  C.  D, 

19  Ul 

- "i . 

Pub.  A.S.P.  53,  200. 

Thomas,  R.  N. 

19k8a 

19U8b 

19U9a 

19k9b 

1950a 

1950b 

1950c 

1952 

Ap,  J.  108,  130. 
Ibid.,  p.  Ik2. 
Ibid.,  109,  k80. 
Ibid.,  110,  12. 
Ibid.,  TIT,  165. 
Ibid.,  HZ,  337. 
Ibid.,  p.  3k3 
Ibid.,  115,  550. 

Thomas,  R,  N.,  Krook, 
M.,  Menzel,  D.  H., 
and  Bhatnagar,  P. 

1953 

Harvard  University  Scientific  Report  No.  8 
(AROC  Contract  AF  19<60k)~lk6) 

Ueno,  S. 

1950 

Publ,  Astr.  Soc.  Japan,  £,  61. 

Uns8ld,  A. 

1952 

Convegno  Volta,  p.  22. 

V/ildt,  R. 

19k? 

Ap.  J.  105,  36. 

Woolley,  R.v.d.R.  and 
Allen,  G.  W. 

1950 

M.N.  no,  358. 

Woolley,  R.v.d.R.  and 
Stibbs,  D.W.N. 

1953 

The  Outer  Layers  of  a Star,  (Oxford: 
University  Press). 

Woltjer,  L. 

195k 

Private  communication. 

Wurm,  K. 

19i*8a 

Zs.  f.  Ap.  25,  109. 

19ii8b 

Mitt.  Bergodorf,  21,  Nos.  202  and  206. 

Zanstra,  H, 

1950a 

Proc.  K.  Nedorl.  Akad.  v,  Wetensch,  53, 
1289. 

1950b 

M.N.  110,  1*91. 

